Следуя примеру с этого веб-сайта, мы можем построить график плотности для нескольких факторов в пределах data.frame
library(sm)
attach(mtcars)
sm.density.compare(mpg, cyl, xlab="Miles Per Gallon")
У меня очень простой вопрос: как добавить вертикальную линию для каждого фактора, который представляет медианное или среднее?
Это пример для среднего. Чтобы вычислить медианное значение, просто замените «FUN = mean» на «FUN = median» в агрегатной функции.
library(sm)
attach(mtcars)
sm.density.compare(mpg, cyl, xlab="Miles Per Gallon")
means <- aggregate(mpg ~ cyl, FUN = mean)
abline(v = means[1,2], col = 2)
abline(v = means[2,2], col = 3, lty = 2)
abline(v = means[3,2], col = 4, lty = 3)
