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.
也可以看看
相关用法
- R beav1 海狸体温系列1
- R biopsy 乳腺癌患者的活检数据
- R boxcox 线性模型的 Box-Cox 变换
- R bandwidth.nrd 通过正态参考分布得到的密度()带宽
- R birthwt 与婴儿出生体重低相关的风险因子
- R bcv 带宽选择的有偏交叉验证
- R bacteria 药物治疗后细菌的存在
- R summary.rlm 鲁棒线性模型的总结方法
- R ginv 广义逆矩阵
- R housing 哥本哈根住房条件调查的频率表
- R predict.qda 根据二次判别分析进行分类
- R contr.sdif 连续差异对比编码
- R Melanoma 恶性黑色素瘤的生存率
- R predict.glmmPQL glmmPQL 拟合的预测方法
- R ucv 带宽选择的无偏交叉验证
- R theta.md 估计负二项式的 theta
- R parcoord 平行坐标图
- R rlm 线性模型的稳健拟合
- R npk 经典 N、P、K 阶乘实验
- R Cars93 1993 年美国销售的 93 辆汽车的数据
- R predict.lda 通过线性判别对多变量观测值进行分类
- R geyser 老忠实间歇泉数据
- R summary.negbin 类“negbin”对象的摘要方法函数
- R Aids2 澳大利亚艾滋病生存数据
- R truehist 绘制直方图
注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Body Temperature Series of Beaver 2。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。