glmer和比值比

Question

有没有办法计算glmer模型的优势比？

我的模型定义如下：

results_reduced <- glmer(R0A1 ~ MPHW_Perc + AG_Perc + Shrub_Perc + Dist_PrimaryRoads 
                         + Dist_SecondaryRoads + (1 | ID), data = secondorder_st,
                         family = binomial)

我想计算这个模型的优势比和相关的置信区间;但是，我将在GLM中使用的以下语法继续运行而不停止：

#odds ratios and 95% CI
exp(cbind(OR = coef(results_reduced), confint(results_reduced)))

因此，我认为随机效应会导致一些挂起问题。有没有办法输出混合效应模型的优势比？

Answer 1

随机效应可能不是问题所在。 1，coef(glmm_model)拉出固定和随机效果的列表。你想要fixef(glmm_model)。你的数据集有多大？计算glmm模型的置信区间需要时间（有时是几分钟），因此您可能不会等待足够长的时间。 3，对于glmer，计算置信区间的默认方式是＆＃34; profile。＆＃34;根据您的数据，它会经常发出警告，当您嵌入cbind和exp的来电时，这些警告可能不合作。你可以尝试confint(..., method ="Wald")，或者一件一件地做事。