如何计算列中每个类别的份额?

时间:2016-11-10 08:36:57

标签: r ggplot2 bar-chart geom-bar

df = data.frame(week = as.factor(rep(c(1, 2), times = 5)),
                name = as.factor(rep(LETTERS[1:5], times = 2)),
                count = rpois(n = 10, lambda = 20))

    > df
     week   name count
1       1      A    16
2       2      B    14
3       1      C    23
4       2      D    15
5       1      E    12
6       2      A    15
7       1      B    23
8       2      C    22
9       1      D    22
10      2      E    26

我想计算每个姓名的每周份数。 起初我打算使用以下方法:

transform(df, week1_share = ifelse(week == "1", round((df$count / sum(df$count)  * 100),2), NA))
transform(df, week2_share = ifelse(week == "2", round((df$count / sum(df$count)  * 100),2), NA))

然后让每一列合并,最终把它作为条形图上的标签,似乎效率太低。必须有一些类型的快速解决方案,我还不知道。

基本上我想做的是如下所示,但添加可能已按上述方式计算的份额%,以便在每个方框内匹配。

ggplot(df, aes(reorder(week, -count),count, color = "white", group = name, fill = name))+
        geom_bar(position = "stack", stat = "identity") +
        scale_y_continuous(labels=comma)+
        ggthemes::scale_color_tableau()

enter image description here

我不知道为什么重新排序功能经常会失败。如果您有任何关于desc订单排序的提示,请分享。

2 个答案:

答案 0 :(得分:2)

使用splitunsplitprop.table的基础R解决方案将是:

df2 <- unsplit(lapply(split(df, df$week), 
                  function(x){
                    x$prop <- prop.table(x$count)
                    x}
                  ), df$week)

简而言之,split会返回一个data.frames分配列表,其中包含第二个参数,unsplit会返回由分割产生的列表。

使用data.table包,这更短:

library(data.table)
dt <- data.table(df)
dt[, prop := prop.table(count), by=week]

我对dplyr不是很精通,但我确信这也是一个非常简短直接的解决方案。

编辑:这是我在dplyr / magrittr中提出的:

library(dplyr)
df3 <- df %>%
   group_by(week) %>%
   mutate(freq = prop.table(count))

答案 1 :(得分:1)

您使用的数据已被使用:

# Loading the required data
df = data.frame(week = as.factor(rep(c(1, 2), times = 5)),
                name = as.factor(rep(LETTERS[1:5], times = 2)),
                count = rpois(n = 10, lambda = 20))

使用plyr包裹功能,计算了百分比和标签的相对位置。

#Loading the required packages    
library(plyr)
library(ggplot2)

# Calculating the percentages
df = ddply(df, .(week), transform, percent = round(count/sum(count) * 100))

# Calculating the position for plotting
df = ddply(df, .(week), transform, pos = cumsum(percent) - (0.5 * percent))

使用上面计算的信息,绘图已经完成。

# Basic graph
p10 <- ggplot() + geom_bar(aes(y = percent, x = week, fill = name), 
                       data = df, stat="identity")

# Adding data labels
p10 <- p10 + geom_text(data=df, aes(x = week, y = pos, 
                                label = paste0(percent,"%")), size=4)
p10

这是你一直在寻找的吗?

enter image description here