Как получить p-значения для случайных эффектов в glmer

Question

Я хочу проанализировать, когда требования протеста направлены на государство, основываясь на действии и характеристиках уровня страны, используя glmer. Итак, я хотел бы получить p-значения как фиксированных, так и случайных эффектов. Моя модель выглядит так:

targets <- glmer(state ~ ENV + HLH + HRI + LAB + SMO + Capital + 
(1 + rile + parties + rep + rep2 + gdppc + election| Country), 
data = df, family = binomial)

Выходные данные дает мне только Variance & Std.Dev. случайных эффектов, а также корреляции между ними, что имеет смысл для большинства многоуровневых анализов, но не для моих целей. Есть ли способ, которым я могу получить что-то вроде оценок и p-значений для случайных эффектов?

Если это невозможно сделать с помощью R, существует ли какое-либо другое статистическое программное обеспечение, которое дало бы такой вывод?

ОБНОВЛЕНИЕ: Следуя приведенным здесь предложениям, я переместил этот вопрос на перекрестную проверку: https://stats.stackexchange.com/questions/381208/r-how-to-get-estimates-and-p-values-for-random-effects-in-glmer

Answer 1

library(lme4)

library(lattice) 

xyplot(incidence/size ~ period|herd, cbpp, type=c('g','p','l'),
   layout=c(3,5), index.cond = function(x,y)max(y))

gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
          data = cbpp, family = binomial)

summary(gm1)