Как отменить вывод lm fit для того, чтобы сделать коробочный сюжет?

Question

Мне нужно визуализировать данные с помощью boxplot, но он не может сгенерировать список объектов. Я попытался просто использовать unlist для lm-объекта, но он все еще говорит, что данные являются списком. Я читал об этом в документации на 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")