使用R進行迴歸模型時變更對照組的做法(兩個,目前其中之一故障)

騰堯同學發現老師教科書7.2(p.213)所使用的方式sjmisc::to_factor()目前是失效的,經老師確認,的確如此。

# 重設對照組
kao06$partyID <- sjmisc::to_factor(kao06$partyID, ref.lvl="民進黨")
kao06.mod.2 <- glm(turnout ~ gender + age + partyID,
family=binomial,
data=kao06)
summary(kao06.mod.2)

老師已向作者反應此問題,
https://github.com/strengejacke/sjmisc/issues/133
在2021年套件作者尚未解決這問題。

目前解決之道是使用較傳統、不用套件、不用標籤只用數值的作法:
kao06$partyID <- ref_lvl(kao06$partyID, lvl=2)

另一替代做法是relevel(),但不推薦這方法,因為它會把既有的標籤給刪去!
fruit <- factor(c("apple", "banana", "orange", "banana", "apple"))
fruit <- relevel(fruit, ref = "banana")



October 19, 2022 更新:目前仍在這個串中討論著
https://github.com/strengejacke/sjlabelled/issues/56