Question

我正在尝试使用R中的boxplot()函数来显示第10和第95百分位数，而不是默认情况下的1.5 * IQR。人们已经在这个论坛上提供了解决方案，但它们似乎都不适合我。

对于我最初的几个盒子图，我使用了ggplot2，效果很好。然而，最近我发现stat_summary在夏天我正在研究的其中一个网站错误地绘制了平均值。Image generated by ggplot2 我只能使用Boxplot()使用以下内容解决此问题;

 boxplot(mydata$tgm~interaction(mydata$site, mydata$season), 
         las = 1,
         names = c("BFa", "MBa", "STa", "BFsp", "MBsp", "STsp", "BFsu", "MBsu", "STsu", "BFw", "MBw", "STw"), 
         col = c("royalblue", "red", "green", "royalblue", "red", "green", "royalblue", "red", "green", "royalblue", "red", "green"), 
         outline = FALSE, 
         xlab = "Site per season", 
         ylab = "TGM concentration (ng/m3)")

然后：

means <- aggregate(mydata$tgm~interaction(mydata$site, mydata$season), data = mydata, mean)

points(means, col = "black", pch = 16)

我是否还有改变我的论点以便显示第10和第90百分位？

使用boxplot（）

0 个答案: