R Dataset / Package MASS / epil
Attachment | Size |
---|---|
dataset-19894.csv | 14.54 KB |
Documentation |
---|
On this Picostat.com statistics page, you will find information about the epil data set which pertains to Seizure Counts for Epileptics. The epil data set is found in the MASS R package. Try to load the epil data set in R by issuing the following command at the console data("epil"). This may load the data into a variable called epil. If R says the epil 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("epil"). 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 epil 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 epil R data set. The size of this file is about 14,894 bytes. Seizure Counts for EpilepticsDescriptionThall and Vail (1990) give a data set on two-week seizure counts for 59 epileptics. The number of seizures was recorded for a baseline period of 8 weeks, and then patients were randomly assigned to a treatment group or a control group. Counts were then recorded for four successive two-week periods. The subject's age is the only covariate. Usageepil FormatThis data frame has 236 rows and the following 9 columns:
SourceThall, P. F. and Vail, S. C. (1990) Some covariance models for longitudinal count data with over-dispersion. Biometrics 46, 657–671. ReferencesVenables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth Edition. Springer. Examplessummary(glm(y ~ lbase*trt + lage + V4, family = poisson, data = epil), cor = FALSE) epil2 <- epil[epil$period == 1, ] epil2["period"] <- rep(0, 59); epil2["y"] <- epil2["base"] epil["time"] <- 1; epil2["time"] <- 4 epil2 <- rbind(epil, epil2) epil2$pred <- unclass(epil2$trt) * (epil2$period > 0) epil2$subject <- factor(epil2$subject) epil3 <- aggregate(epil2, list(epil2$subject, epil2$period > 0), function(x) if(is.numeric(x)) sum(x) else x[1]) epil3$pred <- factor(epil3$pred, labels = c("base", "placebo", "drug"))contrasts(epil3$pred) <- structure(contr.sdif(3), dimnames = list(NULL, c("placebo-base", "drug-placebo"))) summary(glm(y ~ pred + factor(subject) + offset(log(time)), family = poisson, data = epil3), cor = FALSE)summary(glmmPQL(y ~ lbase*trt + lage + V4, random = ~ 1 | subject, family = poisson, data = epil)) summary(glmmPQL(y ~ pred, random = ~1 | subject, family = poisson, data = epil3)) -- 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
|
R Output | Date |
---|---|
Picostat Output - Correlation Coefficients | Dec 7, 2021 |
Picostat Output - Correlation Coefficients | Dec 7, 2021 |
Picostat Output - Correlation Coefficients | Dec 7, 2021 |
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 |