120 likes | 145 Vues
Mixed Linear Models. An Introductory Tutorial. What we have covered!!!. The Linear model:. Mean Structure or Fixed Effects. Errors. What we have covered. We learned to split the errors into within group, and between group components:. Fixed effects. Random effects. Fixed effects: Linear.
E N D
Mixed Linear Models An Introductory Tutorial
What we have covered!!! • The Linear model: Mean Structure or Fixed Effects Errors
What we have covered • We learned to split the errors into within group, and between group components: Fixed effects Random effects
Fixed effects: Linear Including time as a covariate The Data Assuming linear
Fixed effects: Unstructured Time included as a fixed factor The Data Dummy coded
What about the Random effects Treatment by time interaction. (linear, do the treatments progress differently over the study?) Time as covariate (linear) Treat = 0,1 Time = 0,1,2… Fixed effects What about the errors? They can’t be independent.
Error Structures: possible assumptions • Each time point affects only the time point directly following it. • The correlation between consecutive time points is the same for all time points. (Expect this to be violated when the spacing between time points is unequal) • Equal variances across time (Expect this to be violated when there is a ceiling effect)
Which error structure should we use? • Methods of selection Simple: 1. Use Unstructured -Good option when there are many subjects, and few time points 2. Use Auto-regressive -Good option when there are many equally spaced time points and/or few subjects Data based 1. Use AIC to find error structure -This is an elegant data driven model selection tool
Cont. • Methods NOT to use: 1. Do not use tests of the significance of the covariance parameters: -Remember that testing the random effects parameters can not be done with SPSS 2. Do not pick the structure that makes your research findings significant!! -This is fraud
Akaike’s Information Criterion • Akaike’s Information Criterion (AIC) is a measure of how well the model fits the data that penalizes models with lots of parameters. • A good model should fit the data closely (i.e. have little unaccounted for variance), with the fewest parameters possible. The principle is that the simplest model that provides a good accounting of the data is probably the best one. • The smaller the AIC the better. • So, if a model fit with an AR(1) structure has an AIC of 781.9, and when it is fit with TOEP it has an AIC of 780.5, then we should use TOEP. The extra parameters in the TOEP structure are worth it.
The AIC in SPSS This is the AIC for the model The BIC is also a similar measure that can be used
Our Dataset • The Citalopram study (PI Dr. Zisook) • Does Citalopram reduce the depression in schizophrenic patients with subsyndromal depression • Two Groups: Citalopram vs. Placebo • 3 time points: baseline, week 0, 8, and 12 • Outcome measures: PANSS, CDRS, and HAM-17 • There were two sites, but we will only look at the Cincinnati site.