让X1,…,Xn是独立同分布的随机变量,并定义X¯=X1+X2⋯+Xnn.
Pr{X¯≠0}=1Xii=1,…nXi/X¯X1X¯∼X2X¯∼⋯∼XnX¯.
E[Xi/X¯]E[X1X¯]=E[X2X¯]=⋯=E[XnX¯],
i=1,…,nE[XiX¯]=1n(E[X1X¯]+E[X2X¯]+⋯+E[XnX¯])=1nE[X1X¯+X2X¯+⋯+XnX¯]=1nE[X1+X2+⋯+XnX¯]=1nE[nX¯X¯]=nnE[X¯X¯]=1.
让我们看看是否可以通过简单的蒙特卡洛检查。
x <- matrix(rgamma(10^6, 1, 1), nrow = 10^5)
mean(x[, 3] / rowMeans(x))
[1] 1.00511
很好,结果在重复下不会有太大变化。