On this Picostat.com statistics page, you will find information about the epilepsy data set which pertains to Epilepsy Data . The epilepsy data set is found in the HSAUR R package. You can load the epilepsy data set in R by issuing the following command at the console data("epilepsy"). This will load the data into a variable called epilepsy. If R says the epilepsy data set is not found, you can try installing the package by issuing this command install.packages("HSAUR") and then attempt to reload the data. 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 epilepsy R data set. The size of this file is about 6,683 bytes.

Epilepsy Data

Description

A randomised clinical trial investigating the effect of an anti-epileptic drug.

Usage

data("epilepsy")

Format

A data frame with 236 observations on the following 6 variables.

treatment

the treatment group, a factor with levels placebo and Progabide.

base

the number of seizures before the trial.

age

the age of the patient.

seizure.rate

the number of seizures (response variable).

period

treatment period, an ordered factor with levels 1 to 4.

subject

the patient ID, a factor with levels 1 to 59.

Details

In this clinical trial, 59 patients suffering from epilepsy were randomized to groups receiving either the anti-epileptic drug Progabide or a placebo in addition to standard chemotherapy. The numbers of seizures suffered in each of four, two-week periods were recorded for each patient along with a baseline seizure count for the 8 weeks prior to being randomized to treatment and age. The main question of interest is whether taking progabide reduced the number of epileptic seizures compared with placebo.

Source

P. F. Thall and S. C. Vail (1990), Some covariance models for longitudinal count data with overdispersion. Biometrics, 46, 657–671.

Examples

  data("epilepsy", package = "HSAUR")
  library(lattice)
  dotplot(I(seizure.rate / base) ~ period | subject, data = epilepsy, 
          subset = treatment == "Progabide")
  dotplot(I(seizure.rate / base) ~ period | subject, data = epilepsy, 
          subset = treatment == "Progabide")

--

Dataset imported from https://www.r-project.org.