我有合适的模型,想通过"效果"来绘制效果。 R中的包(第3.3.1节)。到目前为止,一切正常。
dipt.moi <- lmer(log_dipt ~ stemp + shum + swind + scc + ct + mt + (1|plot)
#Plotting:
library(effects)
plot(allEffects(dipt.moi)[1])
#or:
plot(Effect(focal.predictors = "stemp", dipt.moi)
但是当我看到我的情节时,我注意到x轴并没有覆盖我数据的所有范围。它只是&#34;切断&#34;边缘的一些数据。 奇怪的是,这并不适用于我的所有情节,仅适用于某些情况,我无法在此发现特定模式。
请注意我在stemp上的实际数据范围是从-2.1到2.9
我试图调整:
plot(Effect(focal.predictors = "stemp", dipt.moi),xlevels=c(-2.1,2.9))
#and
dlist<-list(L$stemp)
plot(Effect(focal.predictors = "stemp", dipt.moi),xlevels=list(dlist))
绘制了地块但不会改变一点。
其他论坛也提出了这个问题(https://stat.ethz.ch/pipermail/r-help/2014-February/366237.html),但我尝试使用xlevels命令(如上所述),xlim命令,但它没有用。
有谁知道如何在这里调整x轴?也许它甚至不可能?!
问候,卢卡斯。
答案 0 :(得分:0)
它不是plot()
问题,但由于allEffects()
计算从-2到2.您可以使用参数xlevels
来定义范围,但是您不能给它正确。
library(effects); library(lme4)
dipt.moi <- lmer(log_dipt ~ stemp + shum + swind + scc + ct + mt + (1|plot))
dlist <- list(stemp = seq(-2.1, 2.9, 0.4)) # the example values
plot(allEffects(dipt.moi, xlevels = dlist)[1])