geom_errorbar 上限为正均值,下限为负均值
geom_errorbar upper for positive mean and lower for negative mean
我想在 geom_bar 绘图上绘制误差线。我知道您可以调整 aes(ymax, ymin) 以仅显示上部或下部条形图。
但是,我有一些负平均值和正平均值,我想分别显示它们的误差线。
这里有一些示例代码和数据可以玩
library(ggplot)
myData <- aggregate(mtcars$mpg,
by = list(cyl = mtcars$cyl, gears = mtcars$gear),
FUN = function(x) c(mean = mean(x), sd = sd(x),
n = length(x)))
myData <- do.call(data.frame, myData)
myData$se <- myData$x.sd / sqrt(myData$x.n)
colnames(myData) <- c("cyl", "gears", "mean", "sd", "n", "se")
myData$names <- c(paste(myData$cyl, "cyl /",
myData$gears, " gear"))
myData$sign <- c(1, -1, 1, 1, -1, 1, -1, 1)
myDataN <- mutate(myData, mean_new = mean*sign, se_new=se+30)
dodge <- position_dodge(width = 0.9)
limits <- aes(
ymax = myDataN$mean_new + myDataN$se_new,
ymin = myDataN$mean_new - myDataN$se_new)
#ymax = myDataN$mean_new, #for upper
#ymin = myDataN$mean_new) #for lower
p <- ggplot(data = myDataN, aes(x = names, y = mean_new, fill = names))
p + geom_bar(stat = "identity", position = dodge) +
geom_errorbar(limits, position = dodge, width = 0.25) +
theme(axis.text.x=element_blank(), axis.ticks.x=element_blank(),
axis.title.x=element_blank())
谢谢
-深
您可以重新定义限制,使它们以 mean_new 为条件(无论是正数还是负数):
limits <- aes(
ymax = myDataN$mean_new + (myDataN$mean_new > 0)*myDataN$se_new,
ymin = myDataN$mean_new - (myDataN$mean_new < 0)*myDataN$se_new)
根据上面的定义,你得到如下图:
最好去除接触条的胡须。但是,我目前无法弄清楚该怎么做。
我想在 geom_bar 绘图上绘制误差线。我知道您可以调整 aes(ymax, ymin) 以仅显示上部或下部条形图。
但是,我有一些负平均值和正平均值,我想分别显示它们的误差线。
这里有一些示例代码和数据可以玩
library(ggplot)
myData <- aggregate(mtcars$mpg,
by = list(cyl = mtcars$cyl, gears = mtcars$gear),
FUN = function(x) c(mean = mean(x), sd = sd(x),
n = length(x)))
myData <- do.call(data.frame, myData)
myData$se <- myData$x.sd / sqrt(myData$x.n)
colnames(myData) <- c("cyl", "gears", "mean", "sd", "n", "se")
myData$names <- c(paste(myData$cyl, "cyl /",
myData$gears, " gear"))
myData$sign <- c(1, -1, 1, 1, -1, 1, -1, 1)
myDataN <- mutate(myData, mean_new = mean*sign, se_new=se+30)
dodge <- position_dodge(width = 0.9)
limits <- aes(
ymax = myDataN$mean_new + myDataN$se_new,
ymin = myDataN$mean_new - myDataN$se_new)
#ymax = myDataN$mean_new, #for upper
#ymin = myDataN$mean_new) #for lower
p <- ggplot(data = myDataN, aes(x = names, y = mean_new, fill = names))
p + geom_bar(stat = "identity", position = dodge) +
geom_errorbar(limits, position = dodge, width = 0.25) +
theme(axis.text.x=element_blank(), axis.ticks.x=element_blank(),
axis.title.x=element_blank())
谢谢 -深
您可以重新定义限制,使它们以 mean_new 为条件(无论是正数还是负数):
limits <- aes(
ymax = myDataN$mean_new + (myDataN$mean_new > 0)*myDataN$se_new,
ymin = myDataN$mean_new - (myDataN$mean_new < 0)*myDataN$se_new)
根据上面的定义,你得到如下图:
最好去除接触条的胡须。但是,我目前无法弄清楚该怎么做。