Как добавить случайные и / или фиксированные эффекты в регрессию клогов в R

Question

Я планирую запустить регрессию Cloglog в R со случайными и / или фиксированными эффектами для проверки ненаблюдаемой гетерогенности. Теперь я не нашел ничего о том, как сделать это в R только одна функция в Stata. Кто-нибудь знаком с этим?

Это пример моей регрессии в клоглоге:

    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)

Answer 1

Вы можете приспособить модель с фиксированными эффектами к стандартной функции glm. Вам просто нужно создать пустышку по уровню интереса. Например, что-то в этом аромате:

Фиксированные эффекты

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, если вы хотите соответствовать модели без перехвата).

Случайные эффекты

Вы должны использовать другой пакет lme4, который предлагает функцию для этого glmer. Вы можете найти некоторую информацию здесь и здесь