0

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

对于我最初的几个箱形图，我使用了 ggplot2，效果很好。然而，最近我发现stat_summary我在暑假期间研究的一个网站的平均值绘制在错误的位置。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 个百分位数？

4

0 回答 0