如何在 R 中將隨機和/或固定效果添加到 cloglog 回歸中

Question

更新關於固定效應的一個變量的處理問題我計划在 R 中運行阻塞回歸，並使用隨機和/或固定效應來檢查未觀察到的異質性。 現在我還沒有找到任何關於如何在 R 中執行此操作的信息，只有一個 function 在 Stata 中。 有人熟悉這個嗎？ 我有一個我想輸入的變量郵政編碼，這個變量由 205 個級別組成，我有點不知道該怎么做。

這是我的堵塞回歸的一個例子：

    model_simple <- as.formula("completion_yesno ~  ac + ov + UCRate + FirstWeek + LastWeek + DayofWeekSu + DayofWeekMo + DayofWeekTu + DayofWeekWe + DayofWeekTh + DayofWeekFr + MonthofYearJan + MonthofYearFeb + MonthofYearMar + MonthofYearApr +MonthofYearMay+ MonthofYearJun + MonthofYearJul + MonthofYearAug + MonthofYearSep + MonthofYearOct + MonthofYearNov")
clog_simple1 = glm(model_simple,data=cllw,family = binomial(link = cloglog))
summary(clog_simple1)

我嘗試使用來自 package glmmML package 的 glmmboost function 輸入固定效果

> model_simple <- as.formula("completion_yesno ~  ac + ov +
> UCRate + FirstWeek + LastWeek + DayofWeekSu + DayofWeekMo
> + DayofWeekTu + DayofWeekWe + DayofWeekTh + DayofWeekFr + MonthofYearJan + MonthofYearFeb + MonthofYearMar + MonthofYearApr
> +MonthofYearMay+ MonthofYearJun + MonthofYearJul + MonthofYearAug + MonthofYearSep + MonthofYearOct + MonthofYearNov") clog_fe <-
> glmmboot(model_simple, data=cllw,family = binomial(link = cloglog),
> cluster = zip, boot = 5000) summary(clog_fe)

function 只是運行，似乎永遠不會產生結果。 我對每一個小費都感到困惑和高興。

Answer 1

您可以將固定效果 model 與標准glm function 配合使用。 您只需要按興趣級別創建一個虛擬對象。 例如，這種風味的東西：

固定效果

model_FE <- as.formula("completion_yesno ~ factor(groupvar) + ac + ov + UCRate + FirstWeek + LastWeek + DayofWeekSu + DayofWeekMo + DayofWeekTu + DayofWeekWe + DayofWeekTh + DayofWeekFr + MonthofYearJan + MonthofYearFeb + MonthofYearMar + MonthofYearApr +MonthofYearMay+ MonthofYearJun + MonthofYearJul + MonthofYearAug + MonthofYearSep + MonthofYearOct + MonthofYearNov")
glm(model_simple,data=cllw,family = binomial(link = cloglog))

factor(group)將創建 K-1（或 K，如果您想在沒有截距的情況下擬合 model）系數。

隨機效應

您必須使用另一個 package lme4為該 glmer 提出glmer 。 你可以在這里和這里找到一些信息