Wd.var未给出R中预期的结果

Wd.var未给出R中预期的结果,r,R,我正在使用SDMTools包中的wt.var函数来计算加权方差。我试着在网上查找文档,在网上找到了cran.r-project.org/web/packages/SDMTools/SDMTools.pdf。有一个指向的链接。但是,该页面不存在。为什么加权方差计算中存在差异?第三行4.75是加权平均数。但是wt.var也没有使用这个公式。参见维基百科的公式: > wt.var(c(1,8,9),c(0.5,0.25,0.25)) [1] 22.7 > 0.5*(1-6)^2+0.25*

我正在使用SDMTools包中的wt.var函数来计算加权方差。我试着在网上查找文档,在网上找到了cran.r-project.org/web/packages/SDMTools/SDMTools.pdf。有一个指向的链接。但是,该页面不存在。为什么加权方差计算中存在差异?第三行4.75是加权平均数。但是wt.var也没有使用这个公式。

参见维基百科的公式:

> wt.var(c(1,8,9),c(0.5,0.25,0.25))
[1] 22.7
> 0.5*(1-6)^2+0.25*(8-6)^2+0.25*(9-6)^2
[1] 15.75
> 0.5*(1-4.75)^2+0.25*(8-4.5)^2+0.25*(9-4.75)^2 
[1] 14.60938
x1
x1<-c(1,8,9)
w1<-c(0.5,0.25,0.25)
wt.mean<-sum(x1*w1)/sum(w1)
wt.var<-(sum(w1)/(sum(w1)^2-sum(w1^2)))*(sum(w1*(x1-wt.mean)^2))
wt.var
[1] 22.7