library(reshape2) source('../../Data/get_tourette_data.R') DF = get_tourette_data() DF$period_I = DF$period_I_motor + DF$period_I_phonic ## total tics DF$period_II = DF$period_II_motor + DF$period_II_phonic DF$period_III = DF$period_III_motor + DF$period_III_phonic DF$period_I_motor = NULL DF$period_I_phonic = NULL DF$period_II_motor = NULL DF$period_II_phonic = NULL DF$period_III_motor = NULL DF$period_III_phonic = NULL DF$patient = NULL DF_melt = melt(DF, id=c('naloxone_dose'), variable.name='period', value.name='n_tics') DF_melt$period = as.factor(DF_melt$period) m = glm(n_tics ~ naloxone_dose + period, family=poisson, data=DF_melt) print(summary(m))