R beav1 海狸體溫係列1

說明

Reynolds (1994) 說明了威斯康星州中北部海狸加拿大蓖麻 (Castor canadensis) 長期溫度動態研究的一小部分。每 10 分鍾通過遙測技術測量四隻雌性動物的體溫，但這裏使用的是兩隻動物每隻不到一天的一段時間的數據。

用法

beav1

格式

beav1 數據幀有 114 行和 4 列。該 DataFrame 包含以下列：

day

觀察日(自 1990 年初以來的天數)，12 月 12 日至 13 日。

time

觀察時間，格式為 0330，表示淩晨 3.30。

temp

測量的體溫以攝氏度為單位。

activ

靜修處外活動的指標。

注意

22:20 的觀測值缺失。

例子

beav1 <- within(beav1,
               hours <- 24*(day-346) + trunc(time/100) + (time%%100)/60)
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 <- ts(c(beav1$temp[1:82], NA, beav1$temp[83:114]),
           start = 9.5, frequency = 6)
activ <- ts(c(beav1$activ[1:82], NA, beav1$activ[83:114]),
            start = 9.5, frequency = 6)

acf(temp[1:53])
acf(temp[1:53], 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])
alpha <- 0.80
stemp <- as.vector(temp - alpha*lag(temp, -1))
sX <- X[-1, ] - alpha * X[-115,]
beav1.ls <- lm(stemp ~ -1 + sX, na.action = na.omit)
summary(beav1.ls, correlation = FALSE)
rm(temp, activ)

來源

P. S. Reynolds (1994) 海狸體溫的時間序列分析。 Lange, N.、Ryan, L.、Billard, L.、Brilllinger, D.、Conquest, L. 和 Greenhouse, J. eds (1994) 生物測定案例研究第 11 章。紐約：約翰·威利父子公司。

參考

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.

也可以看看

beav2