Learn R Programming

brms (version 2.9.0)

epilepsy: Epileptic seizure counts

Description

Breslow and Clayton (1993) analyze data initially provided by Thall and Vail (1990) concerning seizure counts in a randomized trial of anti-convulsant therapy in epilepsy. Covariates are treatment, 8-week baseline seizure counts, and age of the patients in years.

Usage

epilepsy

Arguments

Format

A data frame of 236 observations containing information on the following 9 variables.

Age

The age of the patients in years

Base

The seizure count at 8-weeks baseline

Trt

Either 0 or 1 indicating if the patient received anti-convulsant therapy

patient

The patient number

visit

The session number from 1 (first visit) to 4 (last visit)

count

The seizure count between two visits

obs

The observation number, that is a unique identifier for each observation

zAge

Standardized Age

zBase

Standardized Base

Examples

Run this code
# NOT RUN {
## poisson regression without random effects. 
fit1 <- brm(count ~ zAge + zBase * Trt, 
            data = epilepsy, family = poisson())
summary(fit1) 
plot(fit1)             
    
## poisson regression with varying intercepts of patients
## as well as normal priors for overall effects parameters.    
fit2 <- brm(count ~ zAge + zBase * Trt + (1|patient), 
            data = epilepsy, family = poisson(), 
            prior = set_prior("normal(0,5)"))
summary(fit2) 
plot(fit2)
# }
# NOT RUN {
 
# }

Run the code above in your browser using DataLab