我有3个时间序列,缺少值(T1,T2和T3)。我将它们合并为一个名为x的动物园对象(您可以从此处下载:https : //www.dropbox.com/s/1xk4qn9hthibqpu/x.rda?dl=0)。
我想使用2组权重(w1和w2)来计算加权平均值。这是我的工作:
library(zoo)
load("x.rda")
w1 <- c(0.38,0.04,0.58)
w2 <- c(0.53,0,0.47)
average1 <- zoo(apply(x, 1, function(rw,w) weighted.mean(rw, weights=w, na.rm=TRUE), w1 ), order.by=index(x) )
average2 <- zoo(apply(x, 1, function(rw,w) weighted.mean(rw, weights=w, na.rm=TRUE), w2 ), order.by=index(x) )
plot(average1-average2)
但是average1和average2是相同的。我究竟做错了什么?
请参阅?weighted.mean
。参数的名称是w
,不是weights
:
avg1 <- avg2 <- x[, 1]
avg1[] <- apply(x, 1, weighted.mean, w = w1, na.rm = TRUE)
avg2[] <- apply(x, 1, weighted.mean, w = w2, na.rm = TRUE)
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句