通過自舉獲得 glmer 系數置信區間

Question

我第一次使用 R 中的混合模型進行統計分析。 由於我的數據由二進制結果變量組成，我設法使用 lme4 package 的 glmer function 構建了邏輯 model，我認為它可以按我想要的方式工作。

我現在的目標是調查我的 model 系數的統計顯着性。 我讀到一般來說，廣義混合模型的最佳方法是引導置信區間，但我還沒有在 R 中找到關於如何做到這一點的良好、清晰的解釋。

有人會有什么建議嗎？ R 中是否有任何包可以加快此過程，或者人們通常為此構建自己的功能？ 我之前沒有真正做過任何引導，所以我很感激一些更深入的答案。

Answer 1

如果要計算參數自舉置信區間，內置功能

confint(fitted_model, method = "boot")

應該工作（見?confint.merMod ）

另請參閱此答案（它說明了用戶定義數量的參數化和非參數化自舉）。

如果你有多個核心，你可以通過添加parallel = "multicore", ncpus = parallel::detectCores()-1 （或其他一些合適數量的核心來使用）來加快速度：參見?lme4::bootMer了解詳情。