如何取消列出lm fit以制作箱形图？

Question

我需要使用boxplot可视化数据，但无法生成列表对象。我尝试在lm对象上简单地使用unlist，但是它仍然说数据是一个列表。我已经在R文档中了解到这一点，取消列出lm fit是一个列表，其中包含各个残差作为成分。我该怎么办？

new_data.ref_conc <- lm(formula = conc~OD, data=new_data)
unlist(new_data.ref_conc)
boxplot(new_data.ref_conc~control+treat, data=new_data)

Answer 1

这是一种对线性回归的残差（ε）和拟合值（yhat）进行箱图绘制的方法。由于您没有提供数据，因此我创建了自己的数据：

set.seed(1)
x <- rexp(100, 1)
y <- 1 + 2*x + rnorm(100) 
lm_obj <- lm(y~x)

plotdata <- data.frame(type = rep(c("res", "yhat"), each = 100),
                       value = c(residuals(lm_obj), fitted(lm_obj)))

boxplot(value~type, data = plotdata, col = c("dodgerblue", "hotpink2"), pch = 16,
        names = c("Residuals", "Fitted Values"), main = "My Boxplot")