Problem hierarchische Regression

[tniemann]

Liebe alle,

ich habe ein Problem:

Wenn ich folgende Syntax durchlaufen lasse:

disc1 <- lm(disc_T1 ~ stud, data = diss)
disc2 <- lm(disc_T1 ~ sex, data = diss)
anova(disc1,disc2)

erhalte ich diese Meldung:

Fehler in anova.lmlist(object, ...) :
Modelle sind nicht alle mit der gleichen Datensatzgröße angepasst worden

Als Lösung recherchierte ich:
disc1 <- lm(disc_T1 ~ stud, data = diss)
disc2 <- lm(disc_T1 ~ sex, data = disc1$model)

anova(disc1,disc2)

dann erhalte ich folgende Fehlermeldung:

Fehler in eval(predvars, data, env) : Objekt 'sex' nicht gefunden

Ich bin leider eine komplette Anfängerin und hoffe auf Euren Rat

[bigben]

Hallo,

disc1 <- lm(disc_T1 ~ stud, data = diss)
disc2 <- lm(disc_T1 ~ sex, data = diss)
anova(disc1,disc2)
[...]Fehler in anova.lmlist(object, ...) :
Modelle sind nicht alle mit der gleichen Datensatzgröße angepasst worden

Das dürfte am ehesten daran liegen, dass stud und sex nicht gleich lang sind, was am ehesten an missing data liegt. Was sagt denn

Code: Alles auswählen

length(na.omit(diss$stud))
length(na.omit(diss$sex))

Als Lösung recherchierte ich:
disc1 <- lm(disc_T1 ~ stud, data = diss)
disc2 <- lm(disc_T1 ~ sex, data = disc1$model)

Du könntest auch einfach eine Arbeitskopie von diss anlegen und darin alle NA weglassen:

Code: Alles auswählen

diss.1 <- na.omit(diss[,c("disc_T1", "sex", "stud")])
disc1 <- lm(disc_T1 ~ stud, data = diss.1)
disc2 <- lm(disc_T1 ~ sex, data = diss.1)

Wenn das das Problem nicht löst, empfehle ich Dir, ein funktionierendes Minimalbeispiel zu posten, wie es z. B. hier erklärt wird: viewtopic.php?f=20&t=11

HTH,
Bernhard

[tniemann]

Dankeschön Bernhard