r - R 中箱形图到误差图的转换

Question

我们如何转换从 生成的箱线图

# Boxplot of MPG by Car Cylinders 
boxplot(mpg~cyl,data=mtcars, main="Car Milage Data", 
    xlab="Number of Cylinders", ylab="Miles Per Gallon")

以最小的努力绘制误差图(第 25 和第 75 分位数)？

Answer 1

@Keniajin 的链接使您走上了正确的道路，但您将需要分位数的函数。这是一个 ggplot 解决方案:

require(ggplot2)

首先，我们创建函数来计算分位数25 en 75

Q25 <- function(x) {quantile(x, .25)}
Q75 <- function(x) {quantile(x, .75)}

然后我们使用 stat_summary 和这些函数来绘制图表。请注意，如果您愿意，您可以将函数替换为 stderr 的中值、最小值、最大值。

ggplot(data=mtcars, aes(x=cyl,y=mpg)) +
  stat_summary(fun.y=mean,fun.ymin=Q25, fun.ymax=Q75) +
  ggtitle("Car Milage Data") +
  xlab("Number of Cylinders") +
  ylab("Miles Per Gallon")