0

假设我有一个数据框空气质量。我做了一个 for 循环来绘制空气质量数据集的所有箱线图。

name <- names(airquality)

classes<-sapply(airquality,class)

airquality[is.na(airquality)] <- 0

for (name in name[classes == 'numeric']) {
  boxplot(airquality[,name])
}

现在我想显示所有分位数值，即第一分位数、中位数、第三分位数和平均值，如下图所示。我在网上搜索了很多，但没有找到任何适合我需要的东西。下面是我想要绘制的所需图表：

在此处输入图像描述

4

2 回答 2

3

这是一个示例，仅使用“风”属性。

B = boxplot(airquality[,"Wind"])
text(1.3, B$stats, B$stats)

IQR = B$stats[4] - B$stats[2]
segments(0.5, c(B$stats[2], B$stats[4]), 0.7, c(B$stats[2], B$stats[4]))
text(0.6, B$stats[3], IQR)
arrows(0.6, B$stats[3]+0.5, 0.6, B$stats[4]-0.1, 0.1)
arrows(0.6, B$stats[3]-0.5, 0.6, B$stats[2]+0.1, 0.1)

风箱图

于 2016-12-18T18:02:41.507 回答
1

使用您的代码：

name <- names(airquality)
classes<-sapply(airquality,class)
airquality[is.na(airquality)] <- 0
for (name in name[classes == 'numeric']) {
  boxplot(airquality[,name])
  text(x=1.25,y=fivenum(airquality[,name]), labels =fivenum(airquality[,name]))
  text(x=0.75,y=median(airquality[,name]), labels=IQR(airquality[,name]))
  arrows(0.77, fivenum(airquality[,name])[2], 0.77, fivenum(airquality[,name])[4], angle= 90 ,length=0.07,code=3)
}

情节在这里，在这个链接中：Boxplot with IQR rule

于 2016-12-19T06:28:12.307 回答