使用循环生成方差估计

Question

我的问题：

• 计算N (1,12)变量上3个观测值的样本。
• 使用该示例计算s^2和variance估算值。
• 重复操作10000次，并使用估算值验证E[s^2-var] < E[ variance estimate- var]

我知道方差方程式：

  

sum（（x-mean（x））^ 2）/（n-1）它是（样本方差）的等式

     

sum（（x-mean（x））^ 2）/（n + 1）是（方差估计）的等式

生成，3个观察结果是问题的一半！ 我认为这就足够了

rnorm(3,1,12)

因此

x <- rnorm(3,1,12)

创建功能：

 s2fun <- function(n){
 s2 <- sum((x-mean(x))^2)/(n-1)
 s2
 }

创建循环：

A <- matrix(0,10000,1)

for (i in 1:10000){
    n <- i
    A[i,] <- s2fun(n)
} 

但每当我检查（A）时，它会生成10000个结果，从脚本

中删除脚本

即使

  

和（A）/ n的   导致inf。