r-如何在glmer中获得随机效应的p值

Question

我想使用glmer，根据行动和国家层面的特征，分析抗议的主张何时针对国家。所以，我想得到固定和随机效应的p值。我的模型是这样的：

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

输出只给我方差&随机效应的标准差，以及它们之间的相关性，这对大多数多层次分析来说是有意义的，但对我的目的来说不是。有没有什么方法可以得到随机效应的估计值和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)