对GLMM使用confint()估计CI会导致zetafun(np,ns)的误差。

Question

我有一个二元变量的GLMM，使用R版本3.6.2和lmerTEST 3.1-2，我想得到置信区间。所以我使用了 confint() 功能。

在这里。

model71<-glmer(fetale_Schwangerschaftskomplikationen~Infektion_in_Schwangerschaft+(1|Pat_ID),data=CED,family="binomial")

有了这个... summary() 函数，我得到的resulst没有一个警告，但对于置信区间我得到。

confint(model71)
Computing profile confidence intervals ...
Error in zetafun(np, ns) : profiling detected new, lower deviance

我试着解决这个问题，就像描述的那样 此处，通过提高 devtol 参数，使用 confint.merMod.但是这个提升对我来说没有用。提升 devtol 到1e-8,1e-7或1e-6给我同样的输出。

confint(model71)
Computing profile confidence intervals ...
Error in zetafun(np, ns) : profiling detected new, lower deviance

有没有人有其他的想法，我，我怎么能解决这个问题？(我很抱歉我的英语不好)谢谢

Answer 1

正如在表扬中所提出的，您可以在下面指定用于生成置信区间的方法。confint.merMod() 随着 method 参数，如 confint.merMod(model, method = "Wald").

选项包括引导(boot)、瓦尔德(Wald)，和型材(profile).一如既往，你的实验设置将决定什么最适合你的usecase。