r - 带有砝码的快速标准偏差

Question

我想使用一个函数，该函数可以快速为我提供矢量广告的标准偏差，并允许我在矢量中包含元素的权重。 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返回 sd的函数。最好我将在data.frame上使用它，例如:

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

Answer 1

library(Hmisc)
sqrt(wtd.var(1:3,c(1,1,3)))
#[1] 0.8944272