R Dataset / Package MASS / beav2
Attachment | Size |
---|---|
dataset-60094.csv | 1.65 KB |
Documentation |
---|
On this Picostat.com statistics page, you will find information about the beav2 data set which pertains to Body Temperature Series of Beaver 2. The beav2 data set is found in the MASS R package. Try to load the beav2 data set in R by issuing the following command at the console data("beav2"). This may load the data into a variable called beav2. If R says the beav2 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("beav2"). 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 beav2 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 beav2 R data set. The size of this file is about 1,693 bytes. Body Temperature Series of Beaver 2DescriptionReynolds (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. Usagebeav2 FormatThe
SourceP. 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. ReferencesVenables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer. See Also
Examplesattach(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))library(nlme) # for gls and corAR1 beav2.gls <- gls(temp ~ activ, data = beav2, corr = corAR1(0.8), method = "ML") summary(beav2.gls) summary(update(beav2.gls, subset = 6:100)) detach("beav2"); rm(temp, activ) -- Dataset imported from https://www.r-project.org. |
Picostat Manual |
---|
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 / azprocedure | March 9, 2018 - 1:06 PM | Dataset |
OpenIntro Statistics Dataset - fheights | August 9, 2020 - 12:25 PM | Dataset |
R Dataset / Package datasets / LakeHuron | March 9, 2018 - 1:06 PM | Dataset |
R Dataset / Package wooldridge / crime1 | March 9, 2018 - 1:06 PM | Dataset |
R Dataset / Package datasets / faithful | March 9, 2018 - 1:06 PM | Dataset |