Niech X1,…,Xn będą niezależnymi i identycznie rozmieszczonymi losowymi zmiennymi i zdefiniuj X¯=X1+X2⋯+Xnn.
Załóżmy, że Pr{X¯≠0}=1 . Ponieważ Xi są identycznie rozmieszczone, symetria mówi nam, że dla i=1,…n , (zależne) zmienne losowe Xi/X¯ mają ten sam rozkład:
X1X¯∼X2X¯∼⋯∼XnX¯.
Jeśli istnieją oczekiwania E[Xi/X¯] (jest to kluczowy punkt), to
E[X1X¯]=E[X2X¯]=⋯=E[XnX¯],
a dla i=1,…,n mamy
E[XiX¯]=1n(E[X1X¯]+E[X2X¯]+⋯+E[XnX¯])=1nE[X1X¯+X2X¯+⋯+XnX¯]=1nE[X1+X2+⋯+XnX¯]=1nE[nX¯X¯]=nnE[X¯X¯]=1.
Zobaczmy, czy możemy to sprawdzić za pomocą prostego Monte Carlo.
x <- matrix(rgamma(10^6, 1, 1), nrow = 10^5)
mean(x[, 3] / rowMeans(x))
[1] 1.00511
Dobrze, a wyniki nie zmieniają się zbytnio przy powtarzaniu.