Smooth Random Effects Distribution in a Linear Mixed Model
- 6 December 2004
- journal article
- Published by Oxford University Press (OUP) in Biometrics
- Vol. 60 (4), 945-953
- https://doi.org/10.1111/j.0006-341x.2004.00250.x
Abstract
Summary A linear mixed model with a smooth random effects density is proposed. A similar approach to P-spline smoothing of Eilers and Marx (1996, Statistical Science11, 89–121) is applied to yield a more flexible estimate of the random effects density. Our approach differs from theirs in that the B-spline basis functions are replaced by approximating Gaussian densities. Fitting the model involves maximizing a penalized marginal likelihood. The best penalty parameters minimize Akaike's Information Criterion employing Gray's (1992, Journal of the American Statistical Association87, 942–951) results. Although our method is applicable to any dimensions of the random effects structure, in this article the two-dimensional case is explored. Our methodology is conceptually simple, and it is relatively easy to fit in practice and is applied to the cholesterol data first analyzed by Zhang and Davidian (2001, Biometrics57, 795–802). A simulation study shows that our approach yields almost unbiased estimates of the regression and the smoothing parameters in small sample settings. Consistency of the estimates is shown in a particular case.Keywords
This publication has 22 references indexed in Scilit:
- Selecting the Number of Knots for Penalized SplinesJournal of Computational and Graphical Statistics, 2002
- Selecting the Number of Knots for Penalized SplinesJournal of Computational and Graphical Statistics, 2002
- Empirical Bayes Estimation via the Smoothing by Roughening ApproachJournal of Computational and Graphical Statistics, 1999
- A Smooth Nonparametric Estimate of a Mixing Distribution Using Mixtures of GaussiansJournal of the American Statistical Association, 1996
- A Smooth Nonparametric Estimate of a Mixing Distribution Using Mixtures of GaussiansJournal of the American Statistical Association, 1996
- A Linear Mixed-Effects Model with Heterogeneity in the Random-Effects PopulationJournal of the American Statistical Association, 1996
- Flexible Methods for Analyzing Survival Data Using Splines, with Applications to Breast Cancer PrognosisJournal of the American Statistical Association, 1992
- Nonparametric Maximum Likelihood Estimation of a Mixing DistributionJournal of the American Statistical Association, 1978
- Nonparametric Maximum Likelihood Estimation of a Mixing DistributionJournal of the American Statistical Association, 1978