Fitting nonlinear models with ARMA errors to biological rhythm data.

Many behavioural and physiological processes are periodic with a period of approximately 24 hours. For descriptive purposes, linear regression using a simple sinusoid with a fixed 24 hour period is sometimes an adequate tool for analysing data from such processes. Inference based on regression models under the assumption of independent and identically distributed errors, however, can often mislead seriously. In this paper we present a general class of models for fitting biological rhythms, with use of higher-order harmonic terms of one or more unknown fundamentals and ARMA processes for the errors. We describe the procedures for model specification and estimation and give the theoretical justification for these procedures. Analysis of a series of human core body temperature illustrates the methodology.

[1]  G. Box,et al.  Distribution of Residual Autocorrelations in Autoregressive-Integrated Moving Average Time Series Models , 1970 .

[2]  Gwilym M. Jenkins,et al.  Time series analysis, forecasting and control , 1972 .

[3]  Chris Chatfield,et al.  Introduction to Statistical Time Series. , 1976 .

[4]  T. W. Anderson,et al.  Statistical analysis of time series , 1972 .

[5]  A. M. Walker,et al.  A Survey of statistical work on the Mackenzie River series of annual Canadian lynx trappings for the , 1977 .

[6]  A. Ronald Gallant,et al.  Nonlinear Regression with Autocorrelated Errors , 1976 .

[7]  M. G. Bulmer,et al.  A Statistical Analysis of the 10-Year Cycle in Canada , 1974 .

[8]  Gail Eugene. Bachman.  Analysis of stationary time series , 1963 .

[9]  E. J. Hannan,et al.  Non-linear time series regression , 1971, Journal of Applied Probability.

[10]  A. Winfree.  Book reviewThe clocks that time us : Martin Moore-Ede, Frank Sulzman, and Charles Fuller Harvard University Press, Cambridge, 1982, 448 pp., $25.00 , 1983 .

[11]  H. White,et al.  Nonlinear Regression with Dependent Observations , 1984 .

[12]  B. Efron,et al.  Assessing the accuracy of the maximum likelihood estimator: Observed versus expected Fisher information , 1978 .

[13]  R. Fisher.  Tests of significance in harmonic analysis , 1929 .

[14]  G. C. Tiao,et al.  Consistent Estimates of Autoregressive Parameters and Extended Sample Autocorrelation Function for Stationary and Nonstationary ARMA Models , 1984 .

[15]  R. Marks,et al.  Predicting hospital census using time series regression methods , 1977 .

[16]  Ruey S. Tsay.  Regression Models with Time Series Errors , 1984 .

[17]  David A. Pierce,et al.  Least squares estimation in the regression model with autoregressive-moving average errors , 1971 .

[18]  G. C. Tiao,et al.  Estimation of time series parameters in the presence of outliers , 1988 .

[19]  D J Kupfer,et al.  REM latency: a psychobiologic marker for primary depressive disease. , 1976, Biological psychiatry.

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

[21]  D. A. Pierce.  Distribution of Residual Autocorrelations in the Regression Model with Autoregressive-Moving Average Errors , 1971 .

[22]  Tong Yl,et al.  Parameter estimation in studying circadian rhythms. , 1976 .

[23]  G. Box,et al.  On a measure of lack of fit in time series models , 1978 .

[24]  D. G. Watts,et al.  Relative Curvature Measures of Nonlinearity , 1980 .

[25]  D. Kupfer,et al.  Reduced cortisol latency in depressive illness. , 1983, Archives of general psychiatry.

[26]  A. I. McLeod,et al.  On the Distribution of Residual Autocorrelations in Box–Jenkins Models , 1978 .

[27]  D. B. Preston.  Spectral Analysis and Time Series , 1983 .