6

我想使用一个函数，它可以快速给我一个矢量广告的标准偏差，允许我在矢量中包含元素的权重。IE

sd(c(1,2,3))     #weights all equal 1
#[1] 1
sd(c(1,2,3,3,3))  #weights equal 1,1,3 respectively
#[1] 0.8944272

对于加权平均值，我可以wt.mean()使用library(SDMTools)例如

>  mean(c(1,2,3))
[1] 2
>     wt.mean(c(1,2,3),c(1,1,1))
[1] 2
> 
>     mean(c(1,2,3,3,3))
[1] 2.4
>     wt.mean(c(1,2,3),c(1,1,3))
[1] 2.4

但该wt.sd功能似乎并没有提供我想要的东西：

>   sd(c(1,2,3))
[1] 1
>     wt.sd(c(1,2,3),c(1,1,1))
[1] 1
>     sd(c(1,2,3,3,3))
[1] 0.8944272
>     wt.sd(c(1,2,3),c(1,1,3))
[1] 1.069045

我期待一个0.8944272从我返回的函数 weighted sd。最好我会在 data.frame 上使用它，例如：

data.frame(x=c(1,2,3),w=c(1,1,3))
4

2 回答 2

9 
library(Hmisc)
sqrt(wtd.var(1:3,c(1,1,3)))
#[1] 0.8944272
于 2013-08-09T10:57:49.510 回答
4

您可以使用rep根据权重复制值。然后，sd可以计算得到的向量。

x <- c(1, 2, 3) # values
w <- c(1, 1, 3) # weights

sd(rep(x, w))
[1] 0.8944272
于 2013-08-09T10:58:59.907 回答