2016-03-28 11 views
1

I eine mehrstufiges logistische Regression leitet, mit der Funktion von glmer Paket lme4 in R. My binomisches Ergebnis (oder response-) Variable wird als cg und codiert.glmer Referenz Ergebnis Zeichenfolge

Meine Frage ist: Wie kann ich wissen, welche Ebene der Ergebnisvariablen als Referenzergebnis genommen wird?

Es scheint die Standardreferenz ist die erste alphabetisch (so c in meinem Fall), ist das wahr?

Ich vermute, ich dies nach der levels und relevel Funktionen:

levels(data$Outcome) 
# [1] "c" "g" 

test <- relevel(data$Outcome, ref = "g") 
# levels(test) 
# [1] "g" "c" 

Dies scheint zu sagen, dass c die Referenz vorher war, aber das ist ziemlich umständlich.

Meine Hausaufgaben:

Für den Fall eines 0/1 codierten Antwort:
ich die Antwort nicht mit ?glmer oder die Online-PDF-Handbuch von lme4 oder diese verwandten Beiträge gefunden
Using glmer for logistic regression, how to verify response reference

Haben Sie die Spitze mit dem ‚recode‘ Funktion aus:
Logistic regression - defining reference level in R

Antwort

0

Sie benötigen zu ändern Referenz-Level-Reihenfolge. This post zeigt, wie es geht. In Ihrem Fall schreiben Sie es auf diese Weise:

data$Outcome <- factor(data$Outcome, levels = c("g", "c")) 

bearbeiten, basierend auf OPs comment- Um Ihre Frage zu beantworten: Ja, sind Faktorstufen standardmäßig alphabetisch. Dieser R-Blogger post diskutiert es mehr.

+0

Nach Ihrem Vorschlag umgekehrt Zeichen für alle Log-Odds Schätzungen (z. B. mein Intercept -1,19 wurde 1,19). Ich erkannte auch, dass "g" mein seltenstes Ergebnis ist, was bedeutet, dass "c" tatsächlich als "0" (was ich als Standard bezeichnete) und "g" als "1" codiert wurde. Mit anderen Worten: Das Ergebnis, das zuerst alphabetisch erscheint, ist in der Tat der "Default"; Die Antwort auf meine Frage lautet "JA". – UncleVinnie