我如何 plot 只是 plotly 或 ggplot2 中的残差？

Question

如果我想 plot 只是残差，我可以

plot(model$residuals)

我会得到一个很好的分散 plot。 我如何在 ggplot2 或 plotly 中做同样的事情？ 我不想 plot 残差与拟合值。

谢谢，阿迪

Answer 1

也许像这样？

model <- lm(Sepal.Width ~ Petal.Length, data = iris)

ggplot(data.frame(x = seq(model$residuals), y = model$residuals)) +
  geom_point(aes(x, y)) +
  labs(x = "Index", y = "Residuals", 
       title = paste("Residuals of", format(model$call)))

这比：

plot(model$residuals)

Answer 2

与@AllanCameron 类似，您可以使用broom package ，它还提供有关df （中的变量）结果的其他选项：

library(ggplot2)
library(broom)
#Data
data("iris")
#Model
m1 <- lm(Sepal.Length~Sepal.Width,data=iris)
df <- augment(m1)
ggplot(df, aes(x = 1:nrow(df), y = .resid)) + geom_point() + xlab('x')