美學的長度必須為1或與數據相同（500）

Question

我正在嘗試從線性回歸模型輸出中重現此箱型圖 。 我的模型和使用的代碼如下。 我收到此錯誤“美學必須為長度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教程可以幫助您准確地了解正在發生的事情，例如我寫的這篇文章 。