美学必须是长度1或与数据相同（500）

Question

我正在尝试从线性回归模型输出中重现此boxplots。我的模型和我使用的代码如下。我收到这个错误“美学必须是长度1或与数据相同（500）。”我究竟做错了什么？请帮忙！

out <- lm(L_SHUCK ~ L_VOLUME+CLASS, data = mydata)

r <- residuals(out)

ggplot(out , aes (x = CLASS , y = r, group = CLASS)) + geom_boxplot()

Answer 1

您无法真正放入lm对象的输出来执行此操作。你需要按照这些方针做点什么：

out = lm(mpg ~ wt + as.factor(am), mtcars)
r = residuals(out)
mtcars$residuals = r

ggplot(mtcars) + geom_boxplot(aes(x = as.factor(am), y = residuals))

我认为一个好的ggplot2教程可以帮助你理解究竟发生了什么，例如： this one I wrote