R中的箱形图到错误图的转换
box plot to error plot conversion in R
我们如何转换从
生成的箱形图
# 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 分位数)?
@Keniajin 的 link 使您走上了正确的轨道,但您将需要用于分位数的函数。这是一个 ggplot 解决方案:
require(ggplot2)
首先,我们创建函数来计算分位数25 en 75
Q25 <- function(x) {quantile(x, .25)}
Q75 <- function(x) {quantile(x, .75)}
然后我们使用 stat_summary
和这些函数制作绘图。请注意,如果您愿意,可以将函数替换为标准错误的中值、最小值、最大值。
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")
我们如何转换从
# 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 分位数)?
@Keniajin 的 link 使您走上了正确的轨道,但您将需要用于分位数的函数。这是一个 ggplot 解决方案:
require(ggplot2)
首先,我们创建函数来计算分位数25 en 75
Q25 <- function(x) {quantile(x, .25)}
Q75 <- function(x) {quantile(x, .75)}
然后我们使用 stat_summary
和这些函数制作绘图。请注意,如果您愿意,可以将函数替换为标准错误的中值、最小值、最大值。
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")