平均的垂直線增加密度圖

Question

繼從這個website的例子中，我們可以生成內data.frame

library(sm) 
attach(mtcars) 
sm.density.compare(mpg, cyl, xlab="Miles Per Gallon") 

我的問題很簡單的幾個因素密度圖，怎麼可以添加一個垂直對於代表的中位數或的每個因子的平均值是？

Answer 1

這是的意思是的一個例子。要計算中位數，只需在集合函數中將「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)