The analysis of designed experiments and longitudinal data by using smoothing splines - Discussion

In designed experiments and in particular longitudinal studies, the aim may be to assess the effect of a quantitative variable such as time on treatment effects. Modelling treatment effects can be complex in the presence of other sources of variation. Three examples are presented to illustrate an approach to analysis in such cases. The first example is a longitudinal experiment on the growth of cows under a factorial treatment structure where serial correlation and variance heterogeneity complicate the analysis. The second example involves the calibration of optical density and the concentration of a protein DNase in the presence of sampling variation and variance heterogeneity. The final example is a multienvironment agricultural field experiment in which a yield-seeding rate relationship is required for several varieties of lupins. Spatial variation within environments, heterogeneity between environments and variation between varieties all need to be incorporated in the analysis. In this paper, the cubic smoothing spline is used in conjunction with fixed and random effects, random coefficients and variance modelling to provide simultaneous modelling of trends and covariance structure. The key result that allows coherent and flexible empirical model building in complex situations is the linear mixed model representation of the cubic smoothing spline. An extension is proposed in which trend is partitioned into smooth and non-smooth components. Estimation and inference, the analysis of the three examples and a discussion of extensions and unresolved issues are also presented.

[1]  C. Ansley,et al.  The Signal Extraction Approach to Nonlinear Regression and Spline Smoothing , 1983 .

[2]  G. Robinson.  That BLUP is a Good Thing: The Estimation of Random Effects , 1991 .

[3]  D. Bates,et al.  Nonlinear mixed effects models for repeated measures data. , 1990, Biometrics.

[4]  G. Wahba.  Bayesian "Confidence Intervals" for the Cross-validated Smoothing Spline , 1983 .

[5]  Robert Kohn,et al.  The Performance of Cross-Validation and Maximum Likelihood Estimators of Spline Smoothing Parameters , 1991 .

[6]  G. Wahba.  Improper Priors, Spline Smoothing and the Problem of Guarding Against Model Errors in Regression , 1978 .

[7]  B. Yandell,et al.  Semi-Parametric Generalized Linear Models. , 1985 .

[8]  Robin Thompson,et al.  Average information REML: An efficient algorithm for variance parameter estimation in linear mixed models , 1995 .

[9]  A. Verbyla,et al.  Modelling Variance Heterogeneity: Residual Maximum Likelihood and Diagnostics , 1993 .

[10]  Marie Davidian,et al.  Some Simple Methods for Estimating Intraindividual Variability in Nonlinear Mixed Effects Models , 1993 .

[11]  P. Diggle.  Analysis of Longitudinal Data , 1995 .

[12]  C. Mcgilchrist.  Estimation in Generalized Mixed Models , 1994 .

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

[14]  Brian R. Cullis,et al.  Genotype by Environment Variance Heterogeneity in a Two-Stage Analysis , 1997 .

[15]  A. Harvey.  Estimating Regression Models with Multiplicative Heteroscedasticity , 1976 .

[16]  G. Wahba.  A Comparison of GCV and GML for Choosing the Smoothing Parameter in the Generalized Spline Smoothing Problem , 1985 .

[17]  L. Skovgaard.  NONLINEAR MODELS FOR REPEATED MEASUREMENT DATA. , 1996 .

[18]  M. Kenward,et al.  Small sample inference for fixed effects from restricted maximum likelihood. , 1997, Biometrics.

[19]  P. Diggle,et al.  Semiparametric models for longitudinal data with application to CD4 cell numbers in HIV seroconverters. , 1994, Biometrics.

[20]  P. A. P. Moran,et al.  Maximum-likelihood estimation in non-standard conditions , 1971, Mathematical Proceedings of the Cambridge Philosophical Society.

[21]  R. Schall.  Estimation in generalized linear models with random effects , 1991 .

[22]  H. Chernoff.  On the Distribution of the Likelihood Ratio , 1954 .

[23]  J. G. Rowell,et al.  Analysing data with repeated observations on each experimental unit , 1976, The Journal of Agricultural Science.

[24]  Michael G. Kenward,et al.  A Method for Comparing Profiles of Repeated Measurements , 1987 .

[25]  Brian R. Cullis,et al.  Accounting for natural and extraneous variation in the analysis of field experiments , 1997 .

[26]  P J Diggle,et al.  Nonparametric estimation of covariance structure in longitudinal data. , 1998, Biometrics.

[27]  M. Aitkin.  Modelling variance heterogeneity in normal regression using GLIM , 1987 .

[28]  James E. Grizzle,et al.  Applications of multivariate analysis of variance to repeated measurements experiments. , 1966 .

[29]  D. Stram,et al.  Variance components testing in the longitudinal mixed effects model. , 1994, Biometrics.

[30]  D. M. Allen,et al.  Analysis of growth and dose response curves. , 1969, Biometrics.

[31]  G. Wahba,et al.  A Correspondence Between Bayesian Estimation on Stochastic Processes and Smoothing by Splines , 1970 .

[32]  N. Breslow,et al.  Approximate inference in generalized linear mixed models , 1993 .

[33]  C. E. Rogers,et al.  Symbolic Description of Factorial Models for Analysis of Variance , 1973 .

[34]  G. Wahba.  Spline models for observational data , 1990 .

[35]  B. Silverman,et al.  Nonparametric regression and generalized linear models , 1994 .

[36]  P. Green.  Penalized Likelihood for General Semi-Parametric Regression Models. , 1987 .

[37]  John Wishart,et al.  GROWTH-RATE DETERMINATIONS IN NUTRITION STUDIES WITH THE BACON PIG, AND THEIR ANALYSIS , 1938 .

[38]  C. Mcgilchrist,et al.  A model for the analysis of growth data from designed experiments. , 1990, Biometrics.

[39]  B. Cullis,et al.  The Analysis of Multistratum and Spatially Correlated Repeated Measures Data , 1992 .

[40]  P. Diggle.  An approach to the analysis of repeated measurements. , 1988, Biometrics.

[41]  J. Ware,et al.  Random-effects models for longitudinal data. , 1982, Biometrics.

[42]  K. Liang,et al.  Asymptotic Properties of Maximum Likelihood Estimators and Likelihood Ratio Tests under Nonstandard Conditions , 1987 .

[43]  Rotha.  GENSTAT 5 Release 3 Reference Manual , 1993 .