EVOLUTION-MANAGER
Edit File: beav1.html
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"><html xmlns="http://www.w3.org/1999/xhtml"><head><title>R: Body Temperature Series of Beaver 1</title> <meta http-equiv="Content-Type" content="text/html; charset=utf-8" /> <link rel="stylesheet" type="text/css" href="R.css" /> </head><body> <table width="100%" summary="page for beav1 {MASS}"><tr><td>beav1 {MASS}</td><td style="text-align: right;">R Documentation</td></tr></table> <h2> Body Temperature Series of Beaver 1 </h2> <h3>Description</h3> <p>Reynolds (1994) describes a small part of a study of the long-term temperature dynamics of beaver <em>Castor canadensis</em> 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. </p> <h3>Usage</h3> <pre> beav1 </pre> <h3>Format</h3> <p>The <code>beav1</code> data frame has 114 rows and 4 columns. This data frame contains the following columns: </p> <dl> <dt><code>day</code></dt><dd> <p>Day of observation (in days since the beginning of 1990), December 12–13. </p> </dd> <dt><code>time</code></dt><dd> <p>Time of observation, in the form <code>0330</code> for 3.30am. </p> </dd> <dt><code>temp</code></dt><dd> <p>Measured body temperature in degrees Celsius. </p> </dd> <dt><code>activ</code></dt><dd> <p>Indicator of activity outside the retreat. </p> </dd> </dl> <h3>Note</h3> <p>The observation at 22:20 is missing. </p> <h3>Source</h3> <p>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) <em>Case Studies in Biometry.</em> New York: John Wiley and Sons. </p> <h3>References</h3> <p>Venables, W. N. and Ripley, B. D. (2002) <em>Modern Applied Statistics with S.</em> Fourth edition. Springer. </p> <h3>See Also</h3> <p><code><a href="beav2.html">beav2</a></code> </p> <h3>Examples</h3> <pre> beav1 <- within(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) rm(temp, activ) </pre> <hr /><div style="text-align: center;">[Package <em>MASS</em> version 7.3-51.4 <a href="00Index.html">Index</a>]</div> </body></html>