Dear R users, I'm trying to add a regression line on my boxplots (something like:boxplot(c(1:3),c(4:6),c(5:8))) But I can't see it. Please help !!! It's not a April fool's joke!!! Yan [[alternative HTML version deleted]]
On Apr 1, 2011, at 11:24 AM, Yan Jiao wrote:> boxplot(c(1:3),c(4:6),c(5:8))df <- structure(list(grp = c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3), val = c(1L, 2L, 3L, 4L, 5L, 6L, 5L, 6L, 7L, 8L)), .Names = c("grp", "val"), row.names = c(NA, -10L), class = "data.frame") boxplot(val~grp,data=df) lm(val~grp, data=df) #Call: #lm(formula = val ~ grp, data = df) #Coefficients: #(Intercept) grp # 0.04348 2.21739 ?abline abline(reg=lm(val~grp, data=df) ) -- David Winsemius, MD West Hartford, CT
Am 01.04.2011 11:24, schrieb Yan Jiao:> Dear R users, > > I'm trying to add a regression line on my boxplots (something > like:boxplot(c(1:3),c(4:6),c(5:8))) But I can't see it. Please help > !!! It's not a April fool's joke!!!Your sample data does not make any sense (at least to me). I would do it as follows: set.seed(8) df <- data.frame(y = rnorm(100), x = factor(rep(1:4))) boxplot(y ~ x, data = df) regline <- lm(y ~ as.numeric(x), data = df) abline(regline) HTH, Bernd