THE SAMPLING VARIANCE OF THE RELATIVE EFFICIENCY OF INDIRECT TO DIRECT SELECTION WHEN USING VARIANCE-COVARIANCE COMPONENTS

A formula to compute the sampling variance of the relative efficiency of indirect to direct selection is developed for the case when the genetic and phenotypic components are estimated from the analyses of variance and covariance. This formula holds provided that the genetic and phenotypic variances and covariances are bounded away from zero. It is also assumed that the number of offspring per mating group, k , is constant. The effect of varying k values and the distribution of the ratio estimator of indirect to direct selection, û, are briefly discussed. An example from an expriment using Tribolium castaneum is also given.     

ESTIMATING THE NUMBER OF VISITS TO THE DOCTOR

The frequency of doctor consultations has direct consequences for health care budgets, yet little statistical analysis of the determinants of doctor visits has been reported. We consider the distribution of the number of visits to the doctor and, in particular, we model its dependence on a number of demographic factors. Examination of the Australian 1995 National Health Survey data reveals that generalized linear Poisson or negative binomial models are inadequate for modelling the mean as a function of covariates, because of excessive zero counts, and a mean‐variance relationship that varies enormously over covariate values. A negative binomial model is used, with parameter values estimated in subgroups according to the discrete combinations of the covariate values. Smoothing splines are then used to smooth and interpolate the parameter values. In effect the mean and the shape parameters are each modelled as (different) functions of gender, age and geographical factors. The estimated regressions for the mean have simple and intuitive interpretations. However, the dependence of the (negative binomial) shape parameter on the covariates is more difficult to interpret and is subject to influence by extreme observations. We illustrate the use of the model by estimating the distribution of the number of doctor consultations in the Statistical Local Area of Ryde, based on population numbers from the 1996 census.     

Bias of ml and reml estimators in regression models with arma errors

Data are simulated for a regression model in which the errors have an autoregressive, moving average structure. The parameters of this structure together with the error variance are estimated using both MLand REML techniques. Average biases of estimators from each technique are reported for a range of true parameter values.     

Least significant spacing for ‘one versus the rest’ normal populations

Consider sample means from k(≥2) normal populations where the variances and sample sizes are equal. The problem is to find the ‘least significant difference’ or ‘spacing’ (LSS) between the two largest means, so that if an observed spacing is larger we have confidence 1 - α that the population with largest sample mean also has the largest population mean.

When the variance is known it is shown that the maximum LSS occurs when k = 2, provided a < .2723. In other words, for any value of k we may use the usual (one-tailed) least significant difference to demonstrate that one population has a population mean greater than (or equal to) the rest.

When the variance is estimated bounds are obtained for the confidence which indicate that this last result is approximately correct.     

Inference based on taguchi's saturated designs

The present paper studies the validity of inferential procedures which follow the Taguchi method, under saturated designs. The distribution of the signal to noise (S/N) ratio Y [ILM0001] is investigated,for normal parent distributions. We further investigate the distribution of orthonormal contrasts of such S/N variables. Finally, we discuss and provide critical values for mod-F tests of significance of parameters, when the k smallest SS values are pooled to serve as error variance estimate     

NONPARAMETRIC ESTIMATION OF GENEWISE VARIANCE FOR MICROARRAY DATA

Estimation of genewise variance arises from two important applications in microarray data analysis: selecting significantly differentially expressed genes and validation tests for normalization of microarray data. We approach the problem by introducing a two-way nonparametric model, which is an extension of the famous Neyman-Scott model and is applicable beyond microarray data. The problem itself poses interesting challenges because the number of nuisance parameters is proportional to the sample size and it is not obvious how the variance function can be estimated when measurements are correlated. In such a high-dimensional nonparametric problem, we proposed two novel nonparametric estimators for genewise variance function and semiparametric estimators for measurement correlation, via solving a system of nonlinear equations. Their asymptotic normality is established. The finite sample property is demonstrated by simulation studies. The estimators also improve the power of the tests for detecting statistically differentially expressed genes. The methodology is illustrated by the data from MicroArray Quality Control (MAQC) project.     

A GROUP OF SPLIT-PLOT DESIGNS WITH A HETEROSCEDASTIC MODEL

For a group of split-plot designs it is assumed that the error variance is constant for a particular experiment but it varies from experiment to experiment. Assuming error variances to be known, estimators of the treatment parameters are obtained by weighted (generalised) least squares method and the corresponding analysis is given. For unknown error variances, an adjustment of the statistics using estimated weights is proposed for removing much of the resulting bias. The adjustment stems from a theorem due to Meier (1953).     

Estimation in one–way classification with heteroscedastic error

For an unbalanced one way calssification under the random effect model the problem of estimation of the fixed effect parameter and the variance is considered. Tje error variance which are funtionally related to the above set of parameters are assumed to fall into k classes with constant error varaince for a class. The asymptotic properties of the proposed estimate is established for increasing number of classes k, assuming the number of observations in the classes form a fixed sequence     

ESTIMATION, PREDICTION AND INFERENCE FOR THE LASSO RANDOM EFFECTS MODEL

The least absolute shrinkage and selection operator (LASSO) can be formulated as a random effects model with an associated variance parameter that can be estimated with other components of variance. In this paper, estimation of the variance parameters is performed by means of an approximation to the marginal likelihood of the observed outcomes. The approximation is based on an alternative but equivalent formulation of the LASSO random effects model. Predictions can be made using point summaries of the predictive distribution of the random effects given the data with the parameters set to their estimated values. The standard LASSO method uses the mode of this distribution as the predictor. It is not the only choice, and a number of other possibilities are defined and empirically assessed in this article. The predictive mode is competitive with the predictive mean (best predictor), but no single predictor performs best across in all situations. Inference for the LASSO random effects is performed using predictive probability statements, which are more appropriate under the random effects formulation than tests of hypothesis.     

A CUSUM control chart for monitoring the variance when parameters are estimated

CUSUM control chart has been widely used for monitoring the process variance. It is usually used assuming that the nominal process variance is known. However, several researchers have shown that the ability of control charts to signal when a process is out of control is seriously affected unless process parameters are estimated from a large in-control Phase I data set. In this paper we derive the run length properties of a CUSUM chart for monitoring dispersion with estimated process variance and we evaluate the performance of this chart by comparing it with the same chart but with assumed known process parameters.     

Product-limit survival functions with correlated survival times

A simple variance estimator for product-limit survival functions is demonstrated for survival times with nested errors. Such data arise whenever survival times are observed within clusters of related observations. Greenwood's formula, which assumes independent observations, is not appropriate in this situation. A robust variance estimator is developed using Taylor series linearized values and the between-cluster variance estimator commonly used in multi-stage sample surveys. A simulation study shows that the between-cluster variance estimator is approximately unbiased and yields confidence intervals that maintain the nominal level for several patterns of correlated survival times. The simulation study also shows that Greenwood's formula underestimates the variance when the survival times are positively correlated within a cluster and yields confidence intervals that are too narrow. Extension to life table methods is also discussed.     

AN EXAMINATION OF ESTIMATED RESIDUALS IN A REGRESSION WITH AN INFINITE ORDER PARAMETRIC MODEL

We consider a linear regression with the error term that obeys an autoregressive model of infinite order and estimate parameters of the models. The parameters of the autoregressive model should be estimated based on estimated residuals obtained by means of the method of ordinary least squares, because the errors are unobservable. The consistency of the coefficients, variance and spectral density of the model obeyed by the error term is shown. Further, we estimate the coefficients of the linear regression by means of the method of estimated generalized least squares. We also show the consistency of the estimator.

     

On testing homogeneity of t normal means against ordered alternatives in r groups

Given k=rt normal populations with a common but unknown variance consisting of t treatments applied to r different groups of units, and supposing that in each group the means are monoto-nically non-decreasing (or non-increasing), then the likelihood ratio test of homogeneity of the means in each group against the simple order alternative is considered. Critical values are provided when one observation is drawn from each of the k populations.     

Prediction in a class of mixed models with two variance components

A best unbiased predictor (BUP) of an arbitrary linear combination of fixed and random effects in mixed linear models is available when the true values of the variance ratios are known. When the true values are unknown, a two-stage predictor, obtained from the BUP by replacing the true values by estimated values, can be used. In this article, exact mean squared errors of two-stage predictors are obtained for a class of mixed models with two variance components that includes the balanced one-way random model and other analysis-of-variance models with proportional frequencies and one balanced random factor.     

Unbiasedness of two-stage estimation and prediction procedures for mixed linear models

The traditional method for estimating or predicting linear combinations of the fixed effects and realized values of the random effects in mixed linear models is first to estimate the variance components and then to proceed as if the estimated values of the variance components were the true values. This two-stage procedure gives unbiased estimators or predictors of the linear combinations provided the data vector is symmetrically distributed about its expected value and provided the variance component estimators are translation-invariant and are even functions of the data vector. The standard procedures for estimating the variance components yield even, translation-invariant estimators.     

The Effect of Estimation of Unknown Degrees of Freedom on Estimation of Remaining Parameters in Spanos’ Conditional t Heteroskedastic Model

This article investigates the effect of estimation of unknown degrees of freedom on efficient estimation of remaining parameters in Spanos’ conditional t heteroskedastic model. We compare by simulation three maximum likelihood estimators (MLEs) of the remaining parameters in the model: the MLE of the remaining parameters when all the parameters are estimated by the MLE, when the degrees of freedom is estimated by a method of moments estimator, and when the degrees of freedom is erroneously specified. The latter two methods are found to perform poorly compared to the former method for the inference of variance parameters in the model. Thus, efficient estimation of degrees of freedom by the MLE is important to estimate efficiently the remaining variance parameters.     

INDEX TO VOLUME 4 (1986)

This article provides an empirical investigation of the risk-neutral variance process and the market price of variance risk implied in the foreign-currency options market. There are three principal contributions. First, the parameters of Heston's mean-reverting square-root stochastic volatility model are estimated using dollar/mark option prices from 1987 to 1992. Second, it is shown that these implied parameters can be combined with historical moments of the dollar/mark exchange rate to deduce an estimate of the market price of variance risk. These estimates are found to be nonzero, time varying, and of sufficient magnitude to imply that the compensation for variance risk is a significant component of the risk premia in the currency market. Finally, the out-of-sample test suggests that the historical variance and the Hull and White implied variance contain no more information than that imbedded in the Heston implied variance.     

Comparison of disease prevalence in two populations under double-sampling scheme with two fallible classifiers

A disease prevalence can be estimated by classifying subjects according to whether they have the disease. When gold-standard tests are too expensive to be applied to all subjects, partially validated data can be obtained by double-sampling in which all individuals are classified by a fallible classifier, and some of individuals are validated by the gold-standard classifier. However, it could happen in practice that such infallible classifier does not available. In this article, we consider two models in which both classifiers are fallible and propose four asymptotic test procedures for comparing disease prevalence in two groups. Corresponding sample size formulae and validated ratio given the total sample sizes are also derived and evaluated. Simulation results show that (i) Score test performs well and the corresponding sample size formula is also accurate in terms of the empirical power and size in two models; (ii) the Wald test based on the variance estimator with parameters estimated under the null hypothesis outperforms the others even under small sample sizes in Model II, and the sample size estimated by this test is also accurate; (iii) the estimated validated ratios based on all tests are accurate. The malarial data are used to illustrate the proposed methodologies.     

Comparing Several Variances with a Control in Direction-mixed Families of Hypotheses

This article deals with simultaneous comparisons of k(k ? 1) populations with a control population with respect to variance in direction-mixed families of hypotheses. Computation of critical constants required for the implementation of the proposed procedures is discussed and selected values of the critical constants are tabulated. Power comparison of the proposed procedures with their existing competitors is carried out using Monte Carlo simulation. An example is given to illustrate the implementation of proposed procedures.     

The mixed model for survey regression estimation

Estimation of the population mean under the regression model with random components is considered. Conditions under which the random components regression estimator is design consistent are given. It is shown that consistency holds when incorrect values are used for the variance components. The regression estimator constructed with model parameters that differ considerably from the true parameters performed well in a Monte Carlo study. Variance estimators for the regression predictor are suggested. A variance estimator appropriate for estimators constructed with a biased estimator for the between-group variance component performed well in the Monte Carlo study.