beav1(MASS)
beav1()所属R语言包:MASS
Body Temperature Series of Beaver 1
1海狸体内温度系列
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Reynolds (1994) describes a small part of a study of the long-term temperature dynamics of beaver Castor canadensis in north-central Wisconsin. Body temperature was measured by telemetry every 10 minutes for four females, but data from a one period of less than a day for each of two animals is used there.
雷诺兹(1994年)描述了一个海狸在北部和中部威斯康星州的蓖麻黄花长期的温度变化研究的一小部分。测量体温遥测四女,每隔10分钟,但有用于数据从一个周期不到一天每两个动物。
用法----------Usage----------
beav1
格式----------Format----------
The beav1 data frame has 114 rows and 4 columns. This data frame contains the following columns:
beav1数据框有114行和4列。这个数据框包含下列资料:
Day of observation (in days since the beginning of 1990), December 12–13.
观察日(自1990年初以来的天),12月12-13日。
Time of observation, in the form 0330 for 3.30am.
时间的观察,在形式0330为上午03时半。
Measured body temperature in degrees Celsius.
测量体温摄氏度。
Indicator of activity outside the retreat.
活动指示灯以外的撤退。
注意----------Note----------
The observation at 22:20 is missing.
22:20观察失踪。
源----------Source----------
P. S. Reynolds (1994) Time-series analyses of beaver body temperatures. Chapter 11 of Lange, N., Ryan, L., Billard, L., Brillinger, D., Conquest, L. and Greenhouse, J. eds (1994) Case Studies in Biometry. New York: John Wiley and Sons.
雷诺兹的PS(1994)的时间序列分析,海狸体温。第11章,朗格,北,瑞安,研究,台球,研究,Brillinger,征服,研究和温室,J. EDS(1994)生物统计学中的案例研究。纽约:John Wiley和Sons。
参考文献----------References----------
Modern Applied Statistics with S. Fourth edition. Springer.
参见----------See Also----------
beav2
beav2
举例----------Examples----------
attach(beav1)
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, cor = FALSE)
detach("beav1"); rm(temp, activ)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|