Two-stage analysis based on a mixed model: large-sample asymptotic theory and small-sample simulation results

SUMMARY A two-stage analysis for the mixed model in which variance components due to the random effects are estimated and used to compute generalized least squares estimates of fixed effects is developed. Large-sample theory is used to establish asymptotic properties. An approximate t test that can be used to test linear contrasts among fixed effects is discussed. Two modest simulations, based on a model for a grazing trial (Burns, Harvey, and Giesbrecht, 1981, Proceedings of 14th International Grassland Conference, J. A. Smith and V. W. Hays (eds), 497-500, Boulder, Colorado: Westview Press; Burns et al., 1983, Agronomy Journal 75, 865-871) are used to show that the asymptotic results are reasonable for small samples.

[1]  Raghu N. Kackar,et al.  Unbiasedness of two-stage estimation and prediction procedures for mixed linear models , 1981 .

[2]  N. Kakwani,et al.  Note on the Unbiasedness of a Mixed Regression Estimator , 1968 .

[3]  C. Radhakrishna Rao,et al.  Minimum variance quadratic unbiased estimation of variance components , 1971 .

[4]  Robert I. Jennrich,et al.  Newton-Raphson and Related Algorithms for Maximum Likelihood Variance Component Estimation , 1976 .

[5]  S. R. Searle,et al.  Restricted Maximum Likelihood (REML) Estimation of Variance Components in the Mixed Model , 1976 .

[6]  S. R. Searle Linear Models , 1971 .

[7]  P. A. W. Lewis,et al.  A Pseudo-Random Number Generator for the System/360 , 1969, IBM Syst. J..

[8]  H. Hartley,et al.  Computing Maximum Likelihood Estimates for the Mixed A.O.V. Model Using the W Transformation , 1973 .

[9]  J. Rao,et al.  Combining Independent Estimators and Estimation in Linear Regression with Unequal Variances , 1971 .

[10]  Francis G. Giesbrecht,et al.  An efficient procedure for computing minque of variance components and generalized least squares estimates of fixed effects , 1983 .

[11]  J. Scherpenisse,et al.  A pseudo random number generator , 1974 .

[12]  John F. Monahan,et al.  Monte Carlo Comparison of ANOVA, MIVQUE, REML, and ML Estimators of Variance Components , 1984 .

[13]  M. Kendall,et al.  The advanced theory of statistics , 1945 .

[14]  Joseph C. Burns,et al.  Central Appalachian Hill Land Pasture Evaluation Using Cows and Calves. I. Ordinary and Generalized Least Squares Analysis for an Unbalanced Grazing Experiment1 , 1983 .

[15]  J. N. K. Rao,et al.  On the Estimation of Heteroscedastic Variances , 1973 .

[16]  Calyampudi R. Rao Estimation of Heteroscedastic Variances in Linear Models , 1970 .

[17]  W. Fuller,et al.  Transformations for Estimation of Linear Models with Nested-Error Structure , 1973 .

[18]  F. E. Satterthwaite Synthesis of variance , 1941 .

[19]  H. D. Patterson,et al.  Recovery of inter-block information when block sizes are unequal , 1971 .