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()
我不知道为什么重新排序功能经常会失败。如果您有任何关于desc订单排序的提示,请分享。
答案 0 :(得分:2)
使用split
,unsplit
和prop.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
这是你一直在寻找的吗?