ich rechne derzeit eine Mehrebenenanalyse mit 32 Ländern (Ländervariable "CNTRY3"). Die abhängige Variable ist die Zahlungsbereitschaft für den Umweltschutz "WTP" (Skala 1 - 5, zentriert). Ich habe im dritten Schritt (Random-Intercept-Random-Slope) der Mehrebenenanalyse vier Random Slopes aufgenommen (diese möchte ich auf Grund meiner Theorie variieren lassen):
- RINC_ALL_z -> Einkommen pro Person (z-standardisiert)
- social_trust (Skala 1 - 5, zentriert) -> Vertrauen in andere Menschen
- political_trust (Skala 1 - 5, zentriert) -> Vertrauen in die Politik
- EC_cen (Skala 1 - 5, zentriert) -> Umweltbewusstsein
Code: Alles auswählen
modell.3a <- lmer(WTP ~ RINC_ALL_z + social_trust_cen + political_trust_cen +
Men + lowest_degree + middle_degree + requirement_university +
uncompleted_university + university_degree + AGE_cen + urban +
EC_cen + (RINC_ALL_z + social_trust_cen + political_trust_cen +
EC_cen|CNTRY3), data=ISSP2010_1)
convergence code: 0
Model failed to converge with max|grad| = 0.00527884 (tol = 0.002, component 1)
Nun konnte ich recherchieren, dass wenn man das Modell z. B. mit dem Optimizer "BOBYQA" rechnet, diese "convergence warnings" umgangen werden können. Und tatsächlich, wenn ich das Modell so rechne, dann erscheint kein convergence warning mehr:
Code: Alles auswählen
modell.3b <- lmer(WTP ~ RINC_ALL_z + social_trust_cen + political_trust_cen +
Men + lowest_degree + middle_degree + requirement_university +
uncompleted_university + university_degree + AGE_cen + urban +
EC_cen + (RINC_ALL_z + social_trust_cen + political_trust_cen +
EC_cen|CNTRY3), data=ISSP2010_1, control = lmerControl(optimizer = "bobyqa",
optCtrl=list(maxfun=1e5)))
Bild "loglikmhjf2.png" anzeigen.
Bild "t-valuecomparisont7kr7.png" anzeigen.
Die erste Frage wäre, was so ein Optimizer überhaupt macht (man liest nur immer das man Optimizer benutzen soll, um convergence issues zu vermeiden)?
Die zweite Frage die sich mir stellt, würdet Ihr diese Interpretation der beiden Schaubilder ebenfalls so sehen?
Ich würde mich sehr über eine Antwort freuen, ich zerbreche mir schon mehrere Tage darüber den Kopf...
Viele Grüße
Jörn