堆叠条形图中误差线位置错误（r）

Question

堆叠条形图中误差线位置错误（r）

4

我尝试根据先前回答的问题（在R中使用指定的误差条值制作堆叠条形图）创建一个带有误差棒的堆叠条形图。

然而，我的误差条位置不正确。我尝试更改SD的顺序，这会使误差棒相对于彼此移动，但它们仍然与堆叠条不对齐。

x <- data.frame(Period = c("Mon","Tue","Wed","Thu","Fri","Sat","Sun"), Sample = c("Day","Day","Day","Day","Day","Day","Day"), "12AM-6AM" = c(23.14,16.43,16,22.71,36.86,87.14,110.71), "6AM-12PM" = c(16.14,20.86,18.43,16.71,15.14,14.29,16), "12PM-6PM" = c(26.86,23.71,25.57,23.43,35.29,38,30), "6PM-12AM" =c(35.86,34.14,31.71,36.43,45.57,44,27.14))

library(dplyr)
library(reshape2)

 mx <- melt(x,  id.vars=1:2)
 mx <- mx %>% group_by(Period) %>%
      mutate(pos = cumsum(value)) %>%
      ungroup() %>%
      mutate(sd = c( 5.4, 2.7, 4.7, 4.4, 8.2, 13.2, 20.7,
                    5.6, 2.3, 5.4, 5.5, 1.6, 4.1, 3.1,
                   5.2, 5.6, 5.9, 3.5, 6.3, 4.5, 6.9,
                   3.5, 6.0, 5.9, 6.2, 8.2, 9.0, 2.4
                   ),
             upper = pos + sd/2,
             lower = pos - sd/2)

    days<-c("Mon", "Tue", "Wed","Thu","Fri", "Sat", "Sun")
    ggplot(mx, aes(x=Period, y=value, fill=variable), xLabels=NA) +
      geom_bar(stat="identity") +
      facet_grid(~Sample) + geom_errorbar(aes(ymin = lower, ymax = upper),
 width = .2, col = "red") +
      theme_bw() + scale_x_discrete(limits = days) +
      xlab(" Day of the Week") + scale_fill_grey() +
      ylab ("Number of calls")

此外，我对R非常陌生，对stackoverflow更是初次接触（这是我的第一个问题！）-因此，任何有关如何更清晰地格式化问题或更有效地搜索答案的建议，都将不胜感激 :)。

- Miranda

1个回答

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- Clemens Hug · Accepted Answer

我认为这是因为在y轴上绘制变量时，cumsum(value)中的总和是假设它们是从下往上绘制的。我在进行求和之前重新排列了数据帧，然后它就可以工作了。

x <- data.frame(Period = c("Mon","Tue","Wed","Thu","Fri","Sat","Sun"), Sample = c("Day","Day","Day","Day","Day","Day","Day"), "12AM-6AM" = c(23.14,16.43,16,22.71,36.86,87.14,110.71), "6AM-12PM" = c(16.14,20.86,18.43,16.71,15.14,14.29,16), "12PM-6PM" = c(26.86,23.71,25.57,23.43,35.29,38,30), "6PM-12AM" =c(35.86,34.14,31.71,36.43,45.57,44,27.14))

library(tidyverse)
library(reshape2)
#> 
#> Attaching package: 'reshape2'
#> The following object is masked from 'package:tidyr':
#> 
#>     smiths

mx <- melt(x,  id.vars=1:2)
mx <- mx %>% 
  mutate(sd = c( 5.4, 2.7, 4.7, 4.4, 8.2, 13.2, 20.7,
                 5.6, 2.3, 5.4, 5.5, 1.6, 4.1, 3.1,
                 5.2, 5.6, 5.9, 3.5, 6.3, 4.5, 6.9,
                 3.5, 6.0, 5.9, 6.2, 8.2, 9.0, 2.4)
  ) %>%
  group_by(Period) %>%
  arrange(desc(variable)) %>%
  mutate(
    pos = cumsum(value),
    upper = pos + sd/2,
    lower = pos - sd/2
  ) %>%
  ungroup()

days<-c("Mon", "Tue", "Wed","Thu","Fri", "Sat", "Sun")
ggplot(mx, aes(x=Period, y=value, fill=variable), xLabels=NA) +
  geom_bar(stat="identity") +
  facet_grid(~Sample) + geom_errorbar(aes(ymin = lower, ymax = upper),
                                      width = .2, col = "red") +
  theme_bw() + scale_x_discrete(limits = days) +
  xlab(" Day of the Week") + scale_fill_grey() +
  ylab ("Number of calls")

^{此文档由 reprex package (v0.2.1) 于2018年11月10日创建}