Comparing linear and nonlinear mixed model approaches to cosinor analysis

The cosinor model, used for variables governed by circadian and other biological rhythms, is a nonlinear model in the amplitude and acrophase parameters that has a linear representation upon transformation. With linear cosinor analysis, amplitude and acrophase for each harmonic can be computed as nonlinear functions of the estimated linear regression coefficients. Here a flexible mixed model approach to cosinor analysis is considered, where the fixed effect parameters may enter nonlinearly as acrophase and amplitude for each harmonic or linearly after transformation to regression coefficients. In addition, the random effects may enter nonlinearly as subject-specific deviations from the acrophases and amplitudes or linearly as subject-specific deviations from the regression coefficients. It is also possible for the fixed effects to enter nonlinearly while the random effects enter linearly. Additionally, we evaluate whether including higher order linear harmonic terms as random effects, that is, Rao-Khatri 'covariance adjustment', improves precision. Applying the delta method to nonlinear functions of the parameters from linear mixed cosinor models to obtain approximate variances produces results that are often identical to results from nonlinear mixed models. Consequently, traditional linear cosinor analysis can often be used to estimate and compare the nonlinear parameters of interest, that is, amplitudes and acrophases, via the delta method. This is advantageous since the nonlinear mixed model may have convergence difficulties for more complex models. However, for some multiple-group analyses, the linear cosinor transformation should not be used and we clarify when the two methods are equivalent and when they differ.