Problem bei Dummycodierung in multivariatem Regressionsmodell
Verfasst: Mo Mär 20, 2017 3:08 pm
Hallo zusammen,
für eine Ausarbeitung wollte ich den Einfluss verschiedener Faktoren auf die psychische Gesundheit betrachten. Dabei wollte ich ein multivariates lineares Regressionsmodell schätzen. Für mein Modell habe ich unter anderem eine dichotome Dummyvariable für den Erwerbsstatus erstellt, welche die beiden Ausprägungen erwerbstätig (1) und nicht erwerbstätig (0) umfasst.
Wenn ich nun das folgende Modell schätzen möchte erhalte ich jedoch immer die angegebene Fehlermeldung.
Hier findet ihr auch nochmal meine Syntax, mit der ich die Dummyvariable codiert habe.
Vielen Dank für eure Hilfe und liebe Grüße
Joshua
für eine Ausarbeitung wollte ich den Einfluss verschiedener Faktoren auf die psychische Gesundheit betrachten. Dabei wollte ich ein multivariates lineares Regressionsmodell schätzen. Für mein Modell habe ich unter anderem eine dichotome Dummyvariable für den Erwerbsstatus erstellt, welche die beiden Ausprägungen erwerbstätig (1) und nicht erwerbstätig (0) umfasst.
Wenn ich nun das folgende Modell schätzen möchte erhalte ich jedoch immer die angegebene Fehlermeldung.
Sobald ich die Dummyvariable für den Erwerbsstatus aus dem Modell rausnehme, funktioniert die Berechnung. Jedoch verstehe ich die Fehlermeldung in diesem Zusammenhang nicht, da meine Dummyvariable ja aus zwei Faktoren (0;1) besteht.> multi1 <- lm(psychegesund ~ alter + b.income + factor(frau) + factor(dummy_erwerbsstatus)
+ factor(sammelunterkunft) + factor(verheiratet) + factor(dummy_europa)
+ factor(bmi) + factor(mediensprache_deutsch) + factor(g.humanitaer)
, data = data_c, na.action = na.exclude)
Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) :
Kontraste können nur auf Faktoren mit 2 oder mehr Stufen angewendet werden
Hier findet ihr auch nochmal meine Syntax, mit der ich die Dummyvariable codiert habe.
Ich würde mich freuen, wenn ihr mir bei meinem Problem weiterhelfen könntetdummy_erwerbsstatus <- data_t$bep12
dummy_erwerbsstatus [data_t$bep12 < 0] <- NA
dummy_erwerbsstatus [data_t$bep12 == 9 ] <- 0 #Nicht Erwerbstätig
dummy_erwerbsstatus [data_t$bep12 == 1 |
data_t$bep12 == 2 |
data_t$bep12 == 3 |
data_t$bep12 == 4] <- 1 #Erwerbstätig
Vielen Dank für eure Hilfe und liebe Grüße
Joshua