Значения квантилей построения прямоугольника в 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)
С вашим кодом:
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