小编Gyv*_*yve的帖子

我是R的新手,最近用它做了一些Boxplots.我还在箱线图中添加了均值和标准差.我想知道我是否可以在不同的百分位数中添加某种刻度线或圆圈.假设我想在每个HOUR箱图中标记第85个,第90个百分位数,有没有办法做到这一点？我的数据包括每小时MW的一年负载量和我的输出包括每月每小时24个箱图.我每个月都在做,因为我不确定是否有办法同时运行所有96个(每个月,工作日/周末,4个不同区域)的箱形图.在此先感谢您的帮助.

JANWD <-read.csv("C:\\My Directory\\MWBox2.csv")
JANWD.df<-data.frame(JANWD)
JANWD.sub <-subset(JANWD.df, MONTH < 2 & weekend == "NO")

KeepCols <-c("Hour" , "Houston_Load")
HWD <- JANWD.sub[ ,KeepCols]

sd <-tapply(HWD$Houston_Load, HWD$Hour, sd)
means <-tapply(HWD$Houston_Load, HWD$Hour, mean)

boxplot(Houston_Load ~ Hour, data=HWD, xlab="WEEKDAY HOURS", ylab="MW Differnce", ylim= c(-10, 20), smooth=TRUE ,col ="bisque", range=0)

points(sd, pch = 22, col= "blue")
points(means, pch=23, col ="red")

#Output of the subset of data used to run boxplot for month january in Houston 
str(HWD)
'data.frame':   504 obs. of  2 variables:
 `$ Hour        : …