我有一个使用R版本3.6.2和lmerTEST 3.1-2的带二进制变量的广义线性预测模型,我想得到置信区间。所以我使用了confint()函数。
这里:
model71<-glmer(fetale_Schwangerschaftskomplikationen~Infektion_in_Schwangerschaft+(1|Pat_ID),data=CED,family="binomial")使用summary()函数,我得到了没有警告的结果,但对于置信区间,我得到了:
confint(model71)
Computing profile confidence intervals ...
Error in zetafun(np, ns) : profiling detected new, lower deviance我试图通过使用confint.merMod提升devtol参数来解决上述here所描述的问题。但这种提振对我不起作用。将devtol提升到1e-8、1e-7或1e-6会得到相同的输出:
confint(model71)
Computing profile confidence intervals ...
Error in zetafun(np, ns) : profiling detected new, lower deviance有没有人对我有别的想法,我该如何解决这个问题?(我为我糟糕的英语感到抱歉)谢谢
发布于 2020-05-02 19:11:05
正如推荐中所建议的,您可以使用confint.merMod()和method参数(如confint.merMod(model, method = "Wald") )指定用于生成置信区间的方法。
选项包括引导(boot)、Wald (Wald)和配置文件(profile)。一如既往,你的实验设置将决定什么最适合你的用例。
https://stackoverflow.com/questions/61485322
复制相似问题