A method for fitting regression splines with varying polynomial order in the linear mixed model
- 12 September 2005
- journal article
- research article
- Published by Wiley in Statistics in Medicine
- Vol. 25 (3) , 513-527
- https://doi.org/10.1002/sim.2232
Abstract
The linear mixed model has become a widely used tool for longitudinal analysis of continuous variables. The use of regression splines in these models offers the analyst additional flexibility in the formulation of descriptive analyses, exploratory analyses and hypothesis‐driven confirmatory analyses. We propose a method for fitting piecewise polynomial regression splines with varying polynomial order in the fixed effects and/or random effects of the linear mixed model. The polynomial segments are explicity constrained by side conditions for continuity and some smoothness at the points where they join. By using a reparameterization of this explicitly constrained linear mixed model, an implicitly constrained linear mixed model is constructed that simplifies implementation of fixed‐knot regression splines. The proposed approach is relatively simple, handles splines in one variable or multiple variables, and can be easily programmed using existing commercial software such as SAS or S‐plus. The method is illustrated using two examples: an analysis of longitudinal viral load data from a study of subjects with acute HIV‐1 infection and an analysis of 24‐hour ambulatory blood pressure profiles. Copyright © 2005 John Wiley & Sons, Ltd.Keywords
This publication has 28 references indexed in Scilit:
- Ambulatory Blood Pressure MonitoringSouthern Medical Journal, 2003
- Effect of Dietary Patterns on Ambulatory Blood PressureHypertension, 1999
- Prediction-interval procedures and (fixed-effects) confidence-interval procedures for mixed linear modelsCommunications in Statistics - Theory and Methods, 1988
- Maximum Likelihood Computations with Repeated Measures: Application of the EM AlgorithmJournal of the American Statistical Association, 1987
- [Generalized Additive Models]: CommentStatistical Science, 1986
- ML estimation for the multivariate normal distribution with general linear model mea1 and linear-structure covariance matrix; one-population, complete-data caseCommunications in Statistics - Theory and Methods, 1986
- Splines in StatisticsJournal of the American Statistical Association, 1983
- Maximum Likelihood Approaches to Variance Component Estimation and to Related ProblemsJournal of the American Statistical Association, 1977
- Cubic Splines as a Special Case of Restricted Least SquaresJournal of the American Statistical Association, 1977
- Spline Functions in Data AnalysisTechnometrics, 1974