我想有一个很好的关于我从lm()模型得到的残差的情节.目前我使用plot(model$residuals),但我希望有更好的东西.如果我尝试使用ggplot绘制它,我会收到错误消息:
ggplot2不知道如何处理类numeric的数据
Gop*_*ala 23
不再推荐使用Fortify,根据Hadley的说法,可能会弃用Fortify.
你可以使用扫帚包做类似的事情(更好):
library(broom)
y <-rnorm(10)
x <-1:10
mod <- lm(y ~ x)
df <- augment(mod)
ggplot(df, aes(x = .fitted, y = .resid)) + geom_point()
Run Code Online (Sandbox Code Playgroud)
使用ggfortify::autoplot()的gg回归诊断图的版本.看到这个小插图.
fit <- lm(mpg ~ hp, data = mtcars)
library(ggfortify)
autoplot(fit)
Run Code Online (Sandbox Code Playgroud)