Мне интересно, могу ли я рассчитать устойчивые стандартные ошибки и как я могу это сделать (если это возможно) после многоуровневого квазипуассона, используя glmmpql.
Примером может быть следующий код:
glmmPQL(nabs ~ gen + sco + pedu, random = ~ 1 | region/school, family = quasipoisson, data = data))
После этого, как я могу (или можно) рассчитать устойчивые стандартные ошибки?