我有一个数据框df1:
number=c(4,3,2,3,4,1)
year=c("2000","2000","2000", "2015", "2015", "2015")
items=c(12, 10, 15, 5, 10, 7)
df1=data.frame(number, year, items)
setDT(df1)[, Prop := number/sum(number), by = year]
Run Code Online (Sandbox Code Playgroud)
这样看起来像这样:
number year items Prop
1: 4 2000 12 0.4444444
2: 3 2000 10 0.3333333
3: 2 2000 15 0.2222222
4: 3 2015 5 0.3750000
5: 4 2015 10 0.5000000
6: 1 2015 7 0.1250000
Run Code Online (Sandbox Code Playgroud)
我想获取每年项目数的平均值,因此我尝试使用此功能:
mean.df1=aggregate((df1$number*df1$Prop),list(df1$year), mean)
Run Code Online (Sandbox Code Playgroud)
但返回的均值错误值。我希望它返回:
Group.1 x
1 2000 2.918918
2 2015 2.296296
Run Code Online (Sandbox Code Playgroud)
其中Group.1是年份,x是正确的平均值。
谢谢!
为了aggregate意味着/年项目数
aggregate(number ~ year, data=df1, mean)
# year number
# 1 2000 3.000000
# 2 2015 2.666667
Run Code Online (Sandbox Code Playgroud)
对于基数R中的加权平均值,您可以进行标准的拆分应用合并
sapply(split(df1, df1$year), function(x) weighted.mean(x$number, w=x$items))
Run Code Online (Sandbox Code Playgroud)
要么
sapply(split(df1, df1$year), function(x) sum(x$number*x$items)/sum(x$items))
# 2000 2015
# 2.918919 2.818182
Run Code Online (Sandbox Code Playgroud)