我正在尝试绘制使用lme4包生成的一组二项式glmer对象的预测概率。它们的形式如下:
library(lme4)
model = glmer(
LogFlag ~ FRI + Mound + woody + (1 | ID) +
(0 + FRI | ID) + (0 + Mound | ID) +
(0 + woody | ID) + Mound:FRI + Mound:woody,
family = "binomial",
data = BBmod
)为了绘制图,我使用对visreg的调用,如下所示:
visreg(
model, xvar = "FRI", by = "Mound", scale = "response",
rug = 0, overlay = T, partial = F, gg = T,
breaks = c(2, 4, 6, 8, 10, 15, 20, 40, 80),
band = F, line = list(col = rainbow(9)), legend = F
)这通常会生成非常漂亮的图,x轴上有FRI,y轴上有概率,在breaks参数中指定的Mound值有9行单独的行。
我的问题是,正在使用的第三个预测因子woody的价值是什么?在没有任何breaks参数的情况下,visreg默认使用调用中未包含的预测值的中位数,但我为Mound指定了breaks。如何指定或至少确认正在使用的Woody值?
发布于 2017-05-07 03:04:06
我混淆了breaks和cond参数。文档显示,使用
cond = list(woody = Z)我可以控制woody的值,其中Z是woody的值。也许这会对某人有帮助。
https://stackoverflow.com/questions/43823739
复制相似问题