2016-03-28 80 views
1

我使用lme4中的函數glmer運行多級邏輯迴歸。我的二項結果(或響應)變量編碼爲cgglmer參考結果字符串

我的問題是:我怎麼知道結果變量的哪個級別作爲參考結果?

它似乎默認的引用是第一個字母(所以c在我的情況),這是真的嗎?

我使用的levelsrelevel功能後猜測這:

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

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

這似乎是說,c是之前的參考,但是這是比較間接的。

我的家庭作業:
我沒有找到使用?glmer答案,或lme4在線PDF手冊,或者這些相關文章:

對於0/1編碼的響應的情況下:
Using glmer for logistic regression, how to verify response reference

得到使用「重新編碼」功能,從尖:
Logistic regression - defining reference level in R

回答

0

您需要更改參考級別的順序。 This post演示如何做。在你的情況寫這樣說:

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

編輯基於OP的comment-要回答你的問題:是的,因子水平在默認情況下按字母順序排列。此R-Blogger的post對此進行了更多討論。

+0

根據您的建議,所有對數賠率估計值都會完全顛倒符號(例如,我的攔截-1.19變爲1.19)。我也意識到'g'是我最罕見的結果,這意味着'c'確實被編碼爲'0'(我將其稱爲默認值),'g'編碼爲'1'。換句話說,首先按字母順序排列的結果確實是'默認';我的問題的答案是'是'。 – UncleVinnie