如题,R语言中如何计算样本方差(不用Var函数)?
- 方法一:(不推荐,因为我用SAS用的多,编程习惯很多都是SAS上的习惯,R里面有更简单的方式。比如下面的for循环获取每个值就很SAS)
test<-rnorm(100)
#循环计算的两种方法
a<-0
b<-0
c<-0
for (i in 1:length(test)) {
a=a+test[i]^2
b=b+test[i]
c=c+(test[i]-mean(test))^2
}
#第一种
(a-b^2/length(test))/(length(test)-1)
#第二种
c/(length(test)-1)
这里其实就是通过for循环获取向量中每个值进行计算。这里面的第一种计算公式是:
第二种计算公式是:
- 方法二: