Hi!
I've been working on thses stats for a long moment and I still can't figure out what is the best way to deal with this problem.
Participants in a clinical study have written in a diary emotions and other variables trough a therapy for OCD. I am using LMM in SPSS to analyze a set of emotions (anxiety, sadness, joy) in relation with a treatment variable (outcome) accross therapy sessions (20 weeks) and I am using a very small N (8 persons).
Now, I know I can perform LMM even if there are few missing data, like someone skipped a day here and there, but the problem is that some people haven't rated some emotions. For example, 8 out of 8 rated anxiety and sadness but only 5 rate joy.
If I perform an analysis with all emotions as covariates, SPSS only calculates for the 5 persons who rated all emotions. But I'm wondering if I can analyze all emotions in different models. Is it possible? Or is there a way I can perform an anlysis with everything in one model and not losing participants?
Examples of syntax:
**For everything in one model**
MIXED
IP BY Sujet WITH Semaines Anxiete Tristesse Joie Colère
/PRINT = SOLUTION TESTCOV G
/SAVE = FIXPRED SEFIXP DFFIXP PRED SEPRED DFPRED RESID
/CRITERIA = CIN(95) MXITER(100) MXSTEP(5) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, ABSOLUTE) LCONVERGE(0, ABSOLUTE)
/METHOD = ML
/FIXED = Anxiete Tristesse Joie Colère Semaines Semaines*Anxiete Semaines*Tristesse Semaines*Joie
Semaines*Colère | SSTYPE(3) NOINT
/REPEATED=Semaines | SUBJECT(Sujet) COVTYPE(AR1)
/RANDOM = intercept Semaines | COVTYPE(UN) SUBJECT(Sujet).
**One emotion per model**
MIXED
IP BY Sujet WITH Anxiete Semaines
/PRINT = SOLUTION TESTCOV G
/SAVE = FIXPRED SEFIXP DFFIXP PRED SEPRED DFPRED RESID
/CRITERIA = CIN(95) MXITER(100) MXSTEP(5) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, ABSOLUTE) LCONVERGE(0, ABSOLUTE)
/METHOD = ML
/REPEATED=Semaines | SUBJECT(Sujet) COVTYPE(AR1)
/FIXED = Semaines Anxiete Semaines*Anxiete | SSTYPE(3) NOINT
/RANDOM = intercept Semaines | COVTYPE(UN) SUBJECT(Sujet).
And if I do them seperatly, how can I interpret the effect of time, since I get different results with each model?
Help MUCH appreciated
I've been working on thses stats for a long moment and I still can't figure out what is the best way to deal with this problem.
Participants in a clinical study have written in a diary emotions and other variables trough a therapy for OCD. I am using LMM in SPSS to analyze a set of emotions (anxiety, sadness, joy) in relation with a treatment variable (outcome) accross therapy sessions (20 weeks) and I am using a very small N (8 persons).
Now, I know I can perform LMM even if there are few missing data, like someone skipped a day here and there, but the problem is that some people haven't rated some emotions. For example, 8 out of 8 rated anxiety and sadness but only 5 rate joy.
If I perform an analysis with all emotions as covariates, SPSS only calculates for the 5 persons who rated all emotions. But I'm wondering if I can analyze all emotions in different models. Is it possible? Or is there a way I can perform an anlysis with everything in one model and not losing participants?
Examples of syntax:
**For everything in one model**
MIXED
IP BY Sujet WITH Semaines Anxiete Tristesse Joie Colère
/PRINT = SOLUTION TESTCOV G
/SAVE = FIXPRED SEFIXP DFFIXP PRED SEPRED DFPRED RESID
/CRITERIA = CIN(95) MXITER(100) MXSTEP(5) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, ABSOLUTE) LCONVERGE(0, ABSOLUTE)
/METHOD = ML
/FIXED = Anxiete Tristesse Joie Colère Semaines Semaines*Anxiete Semaines*Tristesse Semaines*Joie
Semaines*Colère | SSTYPE(3) NOINT
/REPEATED=Semaines | SUBJECT(Sujet) COVTYPE(AR1)
/RANDOM = intercept Semaines | COVTYPE(UN) SUBJECT(Sujet).
**One emotion per model**
MIXED
IP BY Sujet WITH Anxiete Semaines
/PRINT = SOLUTION TESTCOV G
/SAVE = FIXPRED SEFIXP DFFIXP PRED SEPRED DFPRED RESID
/CRITERIA = CIN(95) MXITER(100) MXSTEP(5) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, ABSOLUTE) LCONVERGE(0, ABSOLUTE)
/METHOD = ML
/REPEATED=Semaines | SUBJECT(Sujet) COVTYPE(AR1)
/FIXED = Semaines Anxiete Semaines*Anxiete | SSTYPE(3) NOINT
/RANDOM = intercept Semaines | COVTYPE(UN) SUBJECT(Sujet).
And if I do them seperatly, how can I interpret the effect of time, since I get different results with each model?
Help MUCH appreciated