当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R beav2 海狸体温系列2


R语言 beav2 位于 MASS 包(package)。

说明

Reynolds (1994) 说明了威斯康星州中北部海狸加拿大蓖麻 (Castor canadensis) 长期温度动态研究的一小部分。每 10 分钟通过遥测技术测量四只雌性动物的体温,但这里使用的是两只动物每只不到一天的一段时间的数据。

用法

beav2

格式

beav2 数据帧有 100 行和 4 列。该 DataFrame 包含以下列:

day

观察日(自 1990 年初以来的天数),11 月 3 日至 4 日。

time

观察时间,格式为 0330,表示凌晨 3.30。

temp

测量的体温以摄氏度为单位。

activ

静修处外活动的指标。

例子

attach(beav2)
beav2$hours <- 24*(day-307) + trunc(time/100) + (time%%100)/60
plot(beav2$hours, beav2$temp, type = "l", xlab = "time",
   ylab = "temperature", main = "Beaver 2")
usr <- par("usr"); usr[3:4] <- c(-0.2, 8); par(usr = usr)
lines(beav2$hours, beav2$activ, type = "s", lty = 2)

temp <- ts(temp, start = 8+2/3, frequency = 6)
activ <- ts(activ, start = 8+2/3, frequency = 6)
acf(temp[activ == 0]); acf(temp[activ == 1]) # also look at PACFs
ar(temp[activ == 0]); ar(temp[activ == 1])

arima(temp, order = c(1,0,0), xreg = activ)
dreg <- cbind(sin = sin(2*pi*beav2$hours/24), cos = cos(2*pi*beav2$hours/24))
arima(temp, order = c(1,0,0), xreg = cbind(active=activ, dreg))

## IGNORE_RDIFF_BEGIN
library(nlme) # for gls and corAR1
beav2.gls <- gls(temp ~ activ, data = beav2, correlation = corAR1(0.8),
                 method = "ML")
summary(beav2.gls)
summary(update(beav2.gls, subset = 6:100))
detach("beav2"); rm(temp, activ)
## IGNORE_RDIFF_END

来源

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.

也可以看看

beav1

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Body Temperature Series of Beaver 2。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。