A note on maximum Z-scores for control charts for individuals

Control chart limits are often constructed retrospectively based on a sequence of individual measurements. It is shown that the usual control chart limits cannot be crossed for small numbers of measurements.     

On the choice of precedence tests

The best precedence test (BPT) is derived for testing the hypothesis that the lifetimes of two types of items on test have the same distribution. The test has maximum power in the class of the Lehmann type of alternatives F - 1 - (1-G) , A > 1, where F and G are probability distributions of the lifetimes of two types of items on test. This class includes exponential distributions, the Weibull distribution differing only in scale and distributions with proportional hazard rates. Exact power of the BPT is compared with other nonparametrie and parametric tests. The test may terminate before all the lifetimes of the items on test are recorded. In comparing with competing tests of equal size, the power functions are similar but a considerable number of items can be saved and the time on test can be reduced by using the BPT     

Approximating discrete multivariate distributions prom known moments

To approximate the joint distribution of the two-colony stepping-stone model, a finite mixture approach is proposed for constructing discrete multi-variate distributions. This approach generahzes the classic method of linear combinations of independent variables. The stepping-stone model is approximated through matching known moments. Numerical examples from entomology are given. Comparisons are made with the work by Wehrly et al (1993).     

Discrimination of AR,MA and ARMA time series models

The problem of discrimination between two stationary ARMA time series models is considered, and in particular AR(p), MA(p), ARMA(1,1) models. The discriminant based on the likelihood ration leads to a quadratic form that is generally too complicated to evaluated explicitly. The discriminant can be expressed approximately as a linear combination of independent chi–squared random varianles each with one degree of freedom, the coefficients, of which are eigenvalues of cumbersome matrices. An analytical solution which gives the coefficients approximately is suggested.     

Extension to higher dimensions of the jaeschke-eicker result on the standardized empirical process

The asymptotic distribution of the sup-norm of the heavily weighted empirical process is established in the multidimensional case. This theorem extends in particular the famous result in Jaeschke (1975, 1979) to higher dimensions. There is a striking difference between the behaviour for higher dimensions and that for dimension one, especially the limiting distribution is now a simple transformation of a standard exponential random variable.     

An augmented data scoring algorithm for maximum likelihood

The expectation-maximization (EM) method facilitates computation of max¬imum likelihood (ML) and maximum penalized likelihood (MPL) solutions. The procedure requires specification of unobservabie complete data which augment the measured or incomplete data. This specification defines a conditional expectation of the complete data log-likelihood function which is computed in the E-stcp. The EM algorithm is most effective when maximizing the iunction Q{0) denned in the F-stnp is easier than maximizing the likelihood function.

The Monte Carlo EM (MCEM) algorithm of Wei & Tanner (1990) was introduced for problems where computation of Q is difficult or intractable. However Monte Carlo can he computationally expensive, e.g. in signal processing applications involving large numbers of parameters. We provide another approach: a modification of thc standard EM algorithm avoiding computation of conditional expectations.     

Slope-rotatable designs with equal maximum directional variance for second order response surface models

The concept of sloperotaiability with equal maximum directional vari ance for second order response surface models is introduced as a new design property. This requires that the maximum variance of the estimated slope over all possible directions be only a function of p, which is the distance from the design originif is shown that a rotatable design satisfies this property Also, minimization of tiie maximum variance of the estimated slope over all possible directions is proposed as a new design optirnality criterion, and op¬timal designs are called slope-directional minirnax designs. For the class of cquiradial designs, the slope-directional minirnax designs are compared with D— optimal designs.     

On the coefficient of variation for residence time distributions of some stochastic compartmental models

The study of residence time distributions is motivated by the desire to develop new practical tools for the statistical analysis of compartmental systems. In particular, Gibaldi and Perrier (1982) describe three alternative models for a two-compartment system, which were noted to be "indistinguishable based solely on plasma or urinary excretion data," defining a residence time distribution. In this paper, properties of the coefficient of variation of the residence time distributions are developed for these three models. In addition to the standard Markovian model with exponential retention times, properties are also derived for a non-Markovian model with gamma retention times. A coefficient of variation may be estimated from commonly available elimination data, and may be used in principle to discriminate between these three models.     

Umvu estimators for the population mean and variance based on random effects models for lognormal data

Cross-classified data are often obtained in controlled experimental situations and in epidemiologic studies. As an example of the latter, occupational health studies sometimes require personal exposure measurements on a random sample of workers from one or more job groups, in one or more plant locations, on several different sampling dates. Because the marginal distributions of exposure data from such studies are generally right-skewed and well-approximated as lognormal, researchers in this area often consider the use of ANOVA models after a logarithmic transformation. While it is then of interest to estimate original-scale population parameters (e.g., the overall mean and variance), standard candidates such as maximum likelihood estimators (MLEs) can be unstable and highly biased. Uniformly minimum variance unbiased (UMVU) cstiniators offer a viable alternative, and are adaptable to sampling schemes that are typiral of experimental or epidemiologic studies. In this paper, we provide UMVU estimators for the mean and variance under two random effects ANOVA models for logtransformed data. We illustrate substantial mean squared error gains relative to the MLE when estimating the mean under a one-way classification. We illustrate that the results can readily be extended to encompass a useful class of purely random effects models, provided that the study data are balanced.     

Maximum likelihood estimation of parameter structures for the Wishart distribution using constraints

Maximum likelihood estimation under constraints for estimation in the Wishart class of distributions, is considered. It provides a unified approach to estimation in a variety of problems concerning covariance matrices. Virtually all covariance structures can be translated to constraints on the covariances. This includes covariance matrices with given structure such as linearly patterned covariance matrices, covariance matrices with zeros, independent covariance matrices and structurally dependent covariance matrices. The methodology followed in this paper provides a useful and simple approach to directly obtain the exact maximum likelihood estimates. These maximum likelihood estimates are obtained via an estimation procedure for the exponential class using constraints.