用
xa=aggregate(x$avg,by=list(x$value),FUN=weighted.mean,w=x$weight)
给我一个错误
weighted.mean.default(X [[1L]],...)出错:'x'和'w'必须 长度相同
但是
weighted.mean(x$avg,w=x$weight);
工作正常。
答案 0 :(得分:7)
根据old R thread的建议,您可以改为使用by
:
wt <- c(5, 5, 4, 1)/15
x <- c(3.7,3.3,3.5,2.8)
xx <- data.frame(avg=x, value=gl(2,2), weight=wt)
by(xx, xx$value, function(x) weighted.mean(x$avg, x$weight))
答案 1 :(得分:5)
这是一种“百万种方式来涂抹猫”问题,这是一个plyr
解决方案(使用@ chl的示例数据):
ddply(xx,.(value),summarise, wm = weighted.mean(avg,weight))