用ggplots直接绘制统计模型

Question

包visreg允许直接绘制统计模型，我发现这些模型非常方便检查是否有任何错误，并检查我们是否正确理解了估算的含义。我希望将visreg的功能与ggplot的令人难以置信的灵活性结合起来。我希望能够在ggplot代码行中直接调用模型。这是否可行（最终通过直接修改visreg函数）？

---

例如：

require(visreg)
require(ggplot2)
y = c(rnorm(40,10,1), rnorm(20,11,1), rnorm(5,12,1))
x=factor(c(rep(1,40), rep(2,20), rep(3,5))) # this line has changed!
dt=data.frame(x=x, y=y)
m = lm(y~x, data=dt)

我希望能够直接调用m函数中的对象ggplot来表示我的统计模型m。使用visreg，它只是：

visreg(m)

下面的代码不是我要找的，因为它没有直接调用对象m

 ggplot(dt, aes(x,y)) + 
          geom_boxplot(aes(group=x), alpha=0.5)+geom_jitter()

Answer 1

仅供参考，visreg现在支持ggplot2输出：

visreg(m, gg=TRUE)

产生