如何在glmer中获得随机效果的p值

时间:2018-12-08 18:31:12

标签: r lme4 mixed-models p-value multilevel-analysis

我想根据行动和国家/地区的特征,使用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

1 个答案:

答案 0 :(得分: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)