Mixed Model mit zwei Messzeitpunkten

Varianzanalyse, Diskriminanzanalyse, Kontingenzanalyse, Faktorenanalyse, Clusteranalyse, MDS, ....

Moderator: EDi

Antworten
juppi
Beiträge: 19
Registriert: Di Mär 13, 2018 3:55 pm

Mixed Model mit zwei Messzeitpunkten

Beitrag von juppi »

Liebes Forum,

ich habe folgenden Datensatz:

Code: Alles auswählen

head(dat_long)
  patnr Gruppe T0_SOZ_GESCHLECHT_1     T0_SOZ_TEXT_50 Messzeitpunkt Madrs
1     1   2                      2 dt                             1    19
2     1   2                      2 dt                             2     0
3     5   2                      1 dt                             1    16
4     5   2                      1 dt                             2     4
5     7   2                      2 dt                             1    19
6     7   2                      2 dt                             2    16
Ich möcht gerne diese drei verschiedenen Mixed Model berechnen:

Code: Alles auswählen

library(lattice)
library(lme4)
lmm.1 <- lmer(Madrs ~ Messzeitpunkt  + (1 | patnr), dat_long, REML=FALSE)
#Modell mit zufällligen y-Achsenabschnitten
lmm.2 <- lmer(Madrs ~ Messzeitpunkt  + (Messzeitpunkt | patnr), dat_long, REML=FALSE)
#Modell mit zufälligen Steigungen
lmm.3 <- lmer(Madrs ~ Messzeitpunkt*Gruppe  + (Messzeitpunkt | patnr), dat_long, REML=FALSE)
#Modell mit unterschiedlichen Verläufen der Behandlungsgruppen
Leider kommt ab dem zweiten Modell folgende Fehlermeldung:

Code: Alles auswählen

lmm.2 <- lmer(Madrs ~ Messzeitpunkt  + (Messzeitpunkt | patnr), dat_long, REML=FALSE)
Error: number of observations (=266) <= number of random effects (=266) for term (Messzeitpunkt | patnr); the random-effects parameters and the residual variance (or scale parameter) are probably unidentifiable
Ich verstehe die Fehlermeldung leider nicht :( Kann mir damit jemand weiterhelfen?

Vielen Dank!
Benutzeravatar
EDi
Beiträge: 1599
Registriert: Sa Okt 08, 2016 3:39 pm

Re: Mixed Model mit zwei Messzeitpunkten

Beitrag von EDi »

Ich verstehe die Fehlermeldung leider nicht :( Kann mir damit jemand weiterhelfen?
Wie die Fehlermeldung schon sagt, dad Model ist zu komplex für die Daten. Du versuchst 266 random effects zu schätzen (patnr & messzeitpunkt identifizieren genau eine beobachtung) und hast nur 266 beobachtungen. Somit sind keine freiheitsgrade mehr übrig und das model ist nicht schätzbar.

Mehr Daten oder einfacheres Model...
Bitte immer ein reproduzierbares Minimalbeispiel angeben. Meinungen gehören mir und geben nicht die meines Brötchengebers wieder.

Dieser Beitrag ist lizensiert unter einer CC BY 4.0 Lizenz
Bild.
Antworten