R Dataset / Package MASS / beav1
On this Picostat.com statistics page, you will find information about the beav1 data set which pertains to Body Temperature Series of Beaver 1. The beav1 data set is found in the MASS R package. Try to load the beav1 data set in R by issuing the following command at the console data("beav1"). This may load the data into a variable called beav1. If R says the beav1 data set is not found, you can try installing the package by issuing this command install.packages("MASS") and then attempt to reload the data with library("MASS") followed by data("beav1"). Perhaps strangley, if R gives you no output after entering a command, it means the command succeeded. If it succeeded you can see the data by typing beav1 at the command-line which should display the entire dataset.
If you need to download R, you can go to the R project website. You can download a CSV (comma separated values) version of the beav1 R data set. The size of this file is about 1,909 bytes.
Body Temperature Series of Beaver 1
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.
The observation at 22:20 is missing.
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.
Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.
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)
Dataset imported from https://www.r-project.org.
How To Register With a Username
How To Register With Google Single Sign On (SSO)
How To Login With a Username and Password
How To Login With Google Single Sign On (SSO)
How To Import a Dataset
How To Perform Statistical Analysis with Picostat
How To Use Educational Applications with Picostat
|Recent Queries For This Dataset|
No queries made on this dataset yet.
|Title||Authored on||Content type|
|R Dataset / Package COUNT / medpar||March 9, 2018 - 1:06 PM||Dataset|
|R Dataset / Package gap / crohn||March 9, 2018 - 1:06 PM||Dataset|
|R Dataset / Package DAAG / SP500W90||March 9, 2018 - 1:06 PM||Dataset|
|R Dataset / Package datasets / AirPassengers||March 9, 2018 - 1:06 PM||Dataset|
|R Dataset / Package psych / Holzinger||March 9, 2018 - 1:06 PM||Dataset|