我绘制了每年 PM2.5 水平的箱线图。
Boxplot(PM2.5~year, data=subset(dat, hour==12), las=1)
如何从箱线图中提取中位数等值?
默认boxplot
函数不可见地返回摘要,您只需将其分配给一个变量:
res <- boxplot(Sepal.Length ~ Species, data=iris)
其中res
存在一个元素stats
:
> res$stats
[,1] [,2] [,3]
[1,] 4.3 4.9 5.6
[2,] 4.8 5.6 6.2
[3,] 5.0 5.9 6.5
[4,] 5.2 6.3 6.9
[5,] 5.8 7.0 7.9
这些是方框的四分位数摘要。中位数是中间的,所以:
> res$stats[3,]
[1] 5.0 5.9 6.5