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 章。纽约：约翰·威利父子公司。