我试图从使用包glmer
生成的一组二项lme4
个对象中绘制预测概率。它们的形式如下:
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默认使用未包含在调用中的预测变量的中位数,但我已为breaks
提供了Mound
。我如何指定或至少确认正在使用的木质的价值?
答案 0 :(得分:1)
我混淆了breaks
和cond
个参数。文档显示了
cond = list(woody = Z)
我可以控制woody的值,其中Z
是woody
的值。也许这会对某人有所帮助。