Значения квантилей построения прямоугольника в R внутри цикла for

Предположим, у меня есть данные качества воздуха. Я сделал цикл for, чтобы построить всю коробочную диаграмму набора данных о качестве воздуха.

name <- names(airquality)

classes<-sapply(airquality,class)

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

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

Теперь я хочу отобразить все значения квантилей, то есть первый квантиль, медиана, третий квантиль и среднее значение, как показано на рисунке ниже. Я много искал в интернете, но не нашел ничего подходящего мне. Ниже приведен желаемый график, который я хочу построить:

2 ответа

Вот пример, просто используя атрибут "Ветер".

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)

Boxplot of wind

С вашим кодом:

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)
}

Сюжет здесь, по этой ссылке: Боксплот с правилом IQR

Другие вопросы по тегам