R beavers 兩隻海狸的體溫係列

說明

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

用法

beaver1
beaver2

格式

beaver1 DataFrame 有 114 行和 4 列，以 10 分鍾為間隔進行體溫測量。

beaver2 DataFrame 有 100 行和 4 列，以 10 分鍾為間隔進行體溫測量。

變量如下：

天

觀察日(自 1990 年初以來的天數)：12 月 12 日至 13 日 ( beaver1 ) 和 11 月 3 日至 4 日 ( beaver2 )。

時間

觀察時間，格式為0330，淩晨 3:30

溫度

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

活躍的

靜修處外活動的指標。

注意

beaver1 中缺少 22:20 的觀察結果。

例子

require(graphics)
(yl <- range(beaver1$temp, beaver2$temp))

beaver.plot <- function(bdat, ...) {
  nam <- deparse(substitute(bdat))
  with(bdat, {
    # Hours since start of day:
    hours <- time %/% 100 + 24*(day - day[1]) + (time %% 100)/60
    plot (hours, temp, type = "l", ...,
          main = paste(nam, "body temperature"))
    abline(h = 37.5, col = "gray", lty = 2)
    is.act <- activ == 1
    points(hours[is.act], temp[is.act], col = 2, cex = .8)
  })
}
op <- par(mfrow = c(2, 1), mar = c(3, 3, 4, 2), mgp = 0.9 * 2:0)
 beaver.plot(beaver1, ylim = yl)
 beaver.plot(beaver2, ylim = yl)
par(op)

來源

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