我有一个计数的数据框。我想计算加权比例，绘制比例，并绘制这些加权比例的标准误差线。
我的数据框示例：
head(df[1:4,])
badge year total b_1 b_2 b_3 b_4 b_5 b_6 b_7 b_8 b_9 b_10
1 15 2014 14 3 2 1 1 1 1 1 1 1 1
2 15 2015 157 13 12 11 8 6 6 6 5 5 5
3 15 2016 15 5 3 1 1 1 1 1 1 1 0
4 2581 2014 13 1 1 1 1 1 1 1 1 1 1
[该数据包含给定年份中，警察在十个不同的节拍（b_1，b_2 ...）中响应的911电话呼叫计数。因此，第15名警务人员在2014年总共回应了14个电话，其中第3拍属于第1拍，第2拍属于第2拍，依此类推。
本质上，我想要获得每个节拍中发生的通话的总体比例。但我希望这些比例可以按通话总数加权。
到目前为止，我已经能够通过将每个b_列和total列中的值相加并计算比例来计算此值。我已经用简单的条形图绘制了这些。我还无法弄清楚如何计算按总数加权的标准误差。
我不喜欢如何绘制数据。我主要是想获得正确的标准错误。
这里是我到目前为止的代码：
sums_by_beat <- apply(df[, grep('b_', colnames(df2))], 2, sum) props_by_beat <- sums_by_beat / sum(df$total) # Bar plot of proportions by beat barplot(props_by_beat, main='Distribution of Calls by Beat', xlab="Nth Most Common Division", ylim=c(0,1), names.arg=1:length(props_by_beat), ylab="Percent of Total Calls")以及我的数据的30行样本：
df <- structure(list(badge = c(15, 15, 15, 2581, 2581, 2745, 2745,
3162, 3162, 3162, 3396, 3650, 3650, 3688, 3688, 3688, 3698, 3698,
3698, 3717, 3717, 3717, 3740, 3740, 3740, 3813, 3873, 3907, 3930,
4007), year = c(2014, 2015, 2016, 2014, 2015, 2015, 2016, 2014,
2015, 2016, 2016, 2014, 2015, 2014, 2015, 2016, 2014, 2015, 2016,
2014, 2015, 2016, 2014, 2015, 2016, 2016, 2015, 2014, 2014, 2014
), total = c(14, 157, 15, 13, 29, 1, 1, 754, 1172, 1039, 14,
1, 2, 34, 57, 146, 3, 7, 28, 593, 1036, 1303, 461, 952, 1370,
1, 4, 41, 5, 451), b_1 = c(3, 13, 5, 1, 3, 1, 1, 33, 84, 83,
2, 1, 2, 5, 10, 14, 2, 7, 7, 39, 72, 75, 42, 69, 81, 1, 1, 7,
1, 36), b_2 = c(2, 12, 3, 1, 2, 0, 0, 33, 61, 52, 2, 0, 0, 3,
6, 8, 1, 0, 2, 37, 65, 70, 29, 65, 75, 0, 1, 5, 1, 23), b_3 = c(1,
11, 1, 1, 2, 0, 0, 32, 57, 45, 2, 0, 0, 3, 5, 8, 0, 0, 2, 34,
62, 67, 28, 50, 73, 0, 1, 3, 1, 22), b_4 = c(1, 8, 1, 1, 2, 0,
0, 31, 44, 39, 2, 0, 0, 3, 3, 7, 0, 0, 2, 34, 61, 67, 26, 42,
72, 0, 1, 3, 1, 21), b_5 = c(1, 6, 1, 1, 1, 0, 0, 30, 42, 37,
1, 0, 0, 3, 3, 7, 0, 0, 1, 33, 53, 61, 23, 42, 67, 0, 0, 2, 1,
21), b_6 = c(1, 6, 1, 1, 1, 0, 0, 30, 40, 36, 1, 0, 0, 2, 2,
6, 0, 0, 1, 32, 53, 61, 22, 41, 63, 0, 0, 2, 0, 21), b_7 = c(1,
6, 1, 1, 1, 0, 0, 26, 39, 35, 1, 0, 0, 2, 2, 6, 0, 0, 1, 30,
47, 58, 22, 39, 62, 0, 0, 2, 0, 21), b_8 = c(1, 5, 1, 1, 1, 0,
0, 26, 39, 33, 1, 0, 0, 2, 2, 6, 0, 0, 1, 30, 47, 58, 21, 38,
59, 0, 0, 2, 0, 19), b_9 = c(1, 5, 1, 1, 1, 0, 0, 24, 34, 33,
1, 0, 0, 2, 2, 5, 0, 0, 1, 30, 43, 57, 20, 37, 57, 0, 0, 2, 0,
15), b_10 = c(1, 5, 0, 1, 1, 0, 0, 23, 34, 32, 1, 0, 0, 1, 2,
5, 0, 0, 1, 27, 40, 56, 18, 36, 55, 0, 0, 2, 0, 14)), row.names = c(NA,
30L), class = "data.frame")
我有一个计数的数据框。我想计算加权比例，绘制比例，并绘制这些加权比例的标准误差线。我的数据框示例：head（...
您的“比例”实际上不是加权（按total）观测值的平均值吗？然后，我们可以简单地据此计算加权的colMeans。
据我所知，没有内置的R函数来计算加权平均值的标准误，但计算起来相当简单-某些假设在您描述的情况下可能是有效的。例如，请参阅：https://en.wikipedia.org/wiki/Weighted_arithmetic_mean#Standard_error