Reynolds (1994) describes a small part of a study of the long-term temperature dynamics of beaver Castor canadensis in north-central Wisconsin. Body temperature was measured by telemetry every 10 minutes for four females, but data from a one period of less than a day for each of two animals is used there.
The
beav1
data frame has 114 rows and 4 columns.
This data frame contains the following columns:
0330
for 3.30am
The observation at 22:20 is missing.
P. S. Reynolds (1994) Time-series analyses of beaver body temperatures.
Chapter 11 of
Lange, N., Ryan, L., Billard, L., Brillinger, D., Conquest, L.
and Greenhouse, J. eds (1994)
Case Studies in Biometry.
New York: John Wiley and Sons.
beav1 <- beav1 # make local copy attach(beav1) beav1$hours <- 24*(day-346) + trunc(time/100) + (time%%100)/60 detach() plot(beav1$hours, beav1$temp, type = "l", xlab = "time", ylab = "temperature", main = "Beaver 1") usr <- par("usr"); usr[3:4] <- c(-0.2, 8); par(usr=usr) lines(beav1$hours, beav1$activ, type = "s", lty = 2) temp <- rts(c(temp[1:82], NA, temp[83:114]), start = 9.5, frequency = 6, units = "hours") activ <- rts(c(activ[1:82], NA, activ[83:114]), start = 9.5, frequency = 6, units = "hours") acf(temp[1:53]) # and also type="partial" ar(temp[1:53]) act <- c(rep(0, 10), activ) X <- cbind(1, act = act[11:125], act1 = act[10:124], act2 = act[9:123], act3 = act[8:122]) arima.(temp, c(1,0,0), xreg = X)