I would like to detect group effects (if any) along with statistical confidences. I have a hierarchical data set structured as follows:
Drug Groupsconsist of Control, Drug1, and Drug2.
- There are 300
Patientsdistributed amongst the
Drug Groupswith unequal sizes (150, 90, 60 respectively).
Meter Readingsare continuous variables recorded on a daily basis where
Patientswill have at least 150 such readings, though some have more records than others. There's good reason to expect a reading is correlated with the previous day's.
Pooled all together,
Meter Readings aren't normally distributed. Within
Drug Group pools, the
Meter Readings don't have equal variances and are again not normal.
I've decided to pursue a Bayesian route, since the various frequentist comparison of means tests have constraints on normality and/or equal variance. I'm a Pythonist (and a total ignoramus wrt statistics), so I want to use PyMC3 for the modeling. I've checked out Thomas Wiecki's post on Hierarchical Linear Regression as well as Dan Saber's on Multilevel Logistic Regression, but I'm having some trouble getting started.
- Is a multilevel model appropriate for my problem set, considering I only have 3 top level groups?
- How do I go about modeling the nested structure of my data, specifically between
Patientdifferences and the time-dependency in
- How would I interpret the model results to determine if there are measurable differences due to