Standards of Probability Sampling for Legal Evidence

Brief recognition is given to bicentennials of the death of James Stirling (1692–1770) and birth of Ferdinand Rudolph Hassler (1770–1843); to sesquicentennials of the publication of the 3rd definitive edition of Laplace's Théorie analytique des probabilités, and of the births of William Chauvenet (1820–1870) and Isaac Todhunter (1820–1884); to centennials of the births of Louis Bachelier (1870–1946) and Jean Perrin (1870–1942); and to the semicentennial of the initial publication (1920) of Metron, an International Review of Statistics.     

Optimal designs for asymmetrical factorial paired comparison experiments

Three forms of a general null hypothesis Ho on the factorial parameters of a general asymmetrical factorial paired comparison experiment are considered. A class of partially balanced designscorresponding to each form of H₀ is constructed and the A,D and ioptimal design, minimizing the trace, determinant and largest eigenvalue of a defined covariance matrix of related maximumlikelihoodestimators, in that class is determined. Moreover, the optimal design in each class maximizes the noncentrality parameter λ² of the asymptotic noncentral chi-square distribution of the likelihood ratiostatistic -2 log λ for testing H_o under defined local alternatives. These results apply directly to symmetrical factorial paired comparison experiments as special casesExamples are given forillustrating applications of the developed results     

Being Mapuche today in Chile: judgement criteria for ethnic identification among the new generation

This study seeks to identify the underlying factors that sustain ethnic identification judgements – in this case regarding the Mapuche ethnic group – elaborated by a sample of Chilean young people who live in the Santiago metropolitan region. To this end a factorial survey method was used, incorporating quasi-experimental vignettes, which demonstrate social and individual determinants of human judgements. The results show that both self-identification as Mapuche and last name – particularly the paternal surname – are the criteria considered most relevant to inform judgements about ethnic identity. It is clear that there is a complex idea about what it means to be Mapuche today among young people in Chile, combining both subjective and objective factors. This conjunction shows the narrowness of the constructivist and essentialist proposals for understanding ethnicity.     

Finite Sample Comparison of Parametric, Semiparametric, and Wavelet Estimators of Fractional Integration

In this paper we compare through Monte Carlo simulations the finite sample properties of estimators of the fractional differencing parameter, d. This involves frequency domain, time domain, and wavelet based approaches, and we consider both parametric and semiparametric estimation methods. The estimators are briefly introduced and compared, and the criteria adopted for measuring finite sample performance are bias and root mean squared error. Most importantly, the simulations reveal that (1) the frequency domain maximum likelihood procedure is superior to the time domain parametric methods, (2) all the estimators are fairly robust to conditionally heteroscedastic errors, (3) the local polynomial Whittle and bias-reduced log-periodogram regression estimators are shown to be more robust to short-run dynamics than other semiparametric (frequency domain and wavelet) estimators and in some cases even outperform the time domain parametric methods, and (4) without sufficient trimming of scales the wavelet-based estimators are heavily biased.     

A fractional order statistic towards defining a smooth quantile function for discrete data

This work is motivated in part by a recent publication by Ma et al. (2011) who resolved the asymptotic non-normality problem of the classical sample quantiles for discrete data through defining a new mid-distribution based quantile function. This work is the motivation for defining a new and improved smooth population quantile function given discrete data. Our definition is based on the theory of fractional order statistics. The main advantage of our definition as compared to its competitors is the capability to distinguish the uth quantile across different discrete distributions over the whole interval, u∈(0,1). In addition, we define the corresponding estimator of the smooth population quantiles and demonstrate the convergence and asymptotic normal distribution of the corresponding sample quantiles. We verify our theoretical results through a Monte Carlo simulation, and illustrate the utilization of our quantile function in a Q-Q plot for discrete data.     

Cointegration in Fractional Systems with Unknown Integration Orders

Cointegrated bivariate nonstationary time series are considered in a fractional context, without allowance for deterministic trends. Both the observable series and the cointegrating error can be fractional processes. The familiar situation in which the respective integration orders are 1 and 0 is nested, but these values have typically been assumed known. We allow one or more of them to be unknown real values, in which case Robinson and Marinucci (2001, 2003) have justified least squares estimates of the cointegrating vector, as well as narrow‐band frequency‐domain estimates, which may be less biased. While consistent, these estimates do not always have optimal convergence rates, and they have nonstandard limit distributional behavior. We consider estimates formulated in the frequency domain, that consequently allow for a wide variety of (parametric) autocorrelation in the short memory input series, as well as time‐domain estimates based on autoregressive transformation. Both can be interpreted as approximating generalized least squares and Gaussian maximum likelihood estimates. The estimates share the same limiting distribution, having mixed normal asymptotics (yielding Wald test statistics with χ² null limit distributions), irrespective of whether the integration orders are known or unknown, subject in the latter case to their estimation with adequate rates of convergence. The parameters describing the short memory stationary input series are √n‐consistently estimable, but the assumptions imposed on these series are much more general than ones of autoregressive moving average type. A Monte Carlo study of finite‐sample performance is included.     

A Mathematical Model for Pathogen Cross‐Contamination Dynamics during the Postharvest Processing of Leafy Greens

We developed a probabilistic mathematical model for the postharvest processing of leafy greens focusing on Escherichia coli O157:H7 contamination of fresh‐cut romaine lettuce as the case study. Our model can (i) support the investigation of cross‐contamination scenarios, and (ii) evaluate and compare different risk mitigation options. We used an agent‐based modeling framework to predict the pathogen prevalence and levels in bags of fresh‐cut lettuce and quantify spread of E. coli O157:H7 from contaminated lettuce to surface areas of processing equipment. Using an unbalanced factorial design, we were able to propagate combinations of random values assigned to model inputs through different processing steps and ranked statistically significant inputs with respect to their impacts on selected model outputs. Results indicated that whether contamination originated on incoming lettuce heads or on the surface areas of processing equipment, pathogen prevalence among bags of fresh‐cut lettuce and batches was most significantly impacted by the level of free chlorine in the flume tank and frequency of replacing the wash water inside the tank. Pathogen levels in bags of fresh‐cut lettuce were most significantly influenced by the initial levels of contamination on incoming lettuce heads or surface areas of processing equipment. The influence of surface contamination on pathogen prevalence or levels in fresh‐cut bags depended on the location of that surface relative to the flume tank. This study demonstrates that developing a flexible yet mathematically rigorous modeling tool, a “virtual laboratory,” can provide valuable insights into the effectiveness of individual and combined risk mitigation options.     

On parameter estimation for locally stationary long-memory processes

We consider parameter estimation for time-dependent locally stationary long-memory processes. The asymptotic distribution of an estimator based on the local infinite autoregressive representation is derived, and asymptotic formulas for the mean squared error of the estimator, and the asymptotically optimal bandwidth are obtained. In spite of long memory, the optimal bandwidth turns out to be of the order n^-1/5$n^{-1/5}$ and inversely proportional to the square of the second derivative of d. In this sense, local estimation of d is comparable to regression smoothing with iid residuals.     

Statistical inference on the drift parameter in fractional Brownian motion with a deterministic drift

In statistical inference on the drift parameter a in the fractional Brownian motion W^H_t with the Hurst parameter H ∈ (0, 1) with a constant drift Y^H_t = at + W^H_t, there is a large number of options how to do it. We may, for example, base this inference on the properties of the standard normal distribution applied to the differences between the observed values of the process at discrete times. Although such methods are very simple, it turns out that more appropriate is to use inverse methods. Such methods can be generalized to non constant drift. For the hypotheses testing about the drift parameter a, it is more proper to standardize the observed process, and to use inverse methods based on the first exit time of the observed process of a pre-specified interval until some given time. These procedures are illustrated, and their times of decision are compared against the direct approach. Other generalizations are possible when the random part is a symmetric stochastic integral of a known, deterministic function with respect to fractional Brownian motion.