如何在效果图中调整x轴范围

时间:2016-09-16 13:49:37

标签: r plot statistics

我有合适的模型,想通过"效果"来绘制效果。 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;边缘的一些数据。 奇怪的是,这并不适用于我的所有情节,仅适用于某些情况,我无法在此发现特定模式。

effect plot

请注意我在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轴?也许它甚至不可能?!

问候,卢卡斯。

1 个答案:

答案 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])