Smoothing spline estimation of generalised varying-coefficient mixed model

The generalised varying-coefficient model with longitudinal data faces a challenge that the data are correlated, as multiple observations are measured from each individual. In this article we consider the generalised varying-coefficient mixed model (GVCMM) which uses a varying-coefficient model to fit mean functions, while accounting for overdispersion and correlation by adding random effects. Smoothing splines are used to estimate the smooth but arbitrary nonparametric coefficient functions. The usually intractable integration involved in evaluating the quasi-likelihood function is approximated by the Laplace method. This suggests that the GVCMM can be approximately represented by a generalised linear mixed model. Hence, the smoothing parameters and the variance components can be estimated by using the restricted maximum log-likelihood (REML) approach, where the smoothing parameters are treated as an extra variance component vector. We illustrate the performance of the proposed method through some simulation and an application to a real data set.

[1]  Jianqing Fan,et al.  Two‐step estimation of functional linear models with applications to longitudinal data , 1999 .

[2]  L. Tierney,et al.  Accurate Approximations for Posterior Moments and Marginal Densities , 1986 .

[3]  Jianhua Z. Huang,et al.  Varying‐coefficient models and basis function approximations for the analysis of repeated measurements , 2002 .

[4]  Tom L. Burr,et al.  Modeling Longitudinal Data , 2006, Technometrics.

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

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

[7]  N. Breslow,et al.  Bias Correction in Generalized Linear Mixed Models with Multiple Components of Dispersion , 1996 .

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

[9]  D. Harville.  Maximum Likelihood Approaches to Variance Component Estimation and to Related Problems , 1977 .

[10]  Jeffrey S. Morris,et al.  Wavelet‐based functional mixed models , 2006, Journal of the Royal Statistical Society. Series B, Statistical methodology.

[11]  Chong Gu.  Smoothing Spline Anova Models , 2002 .

[12]  Local Asymptotics for B-Spline Estimators of the Varying Coefficient Model , 2004 .

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

[14]  Chin-Tsang Chiang,et al.  Smoothing Spline Estimation for Varying Coefficient Models With Repeatedly Measured Dependent Variables , 2001 .

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

[16]  Li Ping Yang,et al.  Nonparametric smoothing estimates of time-varying coefficient models with longitudinal data , 1998 .

[17]  R. Kohn,et al.  Nonparametric spline regression with prior information , 1993 .

[18]  X. Lin,et al.  Inference in generalized additive mixed modelsby using smoothing splines , 1999 .

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

[20]  R. Tibshirani,et al.  Varying‐Coefficient Models , 1993 .

[21]  Jianqing Fan,et al.  Functional-Coefficient Regression Models for Nonlinear Time Series , 2000 .