主題
Search

樣本方差計算

DOWNLOAD Mathematica Notebook下載 Wolfram 筆記本

當數值計算樣本方差 s 時，必須先計算均值，然後才能確定 s^2。這需要儲存樣本值集。但是，可以使用僅涉及最後一個樣本的遞迴關係來計算 s^2，如下所示。這意味著無需預先計算 mu 本身，並且每一步只需儲存一組執行值。

在下文中，使用略微欠佳的符號 mu_j 來表示從前 j 個樣本計算出的 mu（即，不是第 j 個 矩）

mu_j=(sum_(i=1)^(j)x_i)/j,
(1)

並令 s_j^2 表示從前 j 個樣本計算出的偏差校正樣本方差 s_(N-1)^2 的值。為均值計算出的前幾個值是

mu_1=x_1
(2)
mu_2=(1·mu_1+x_2)/2
(3)
mu_3=(2mu_2+x_3)/3.
(4)

因此，對於 j=2、3，以下等式成立

mu_j=((j-1)mu_(j-1)+x_j)/j.
(5)

因此，透過歸納法，

mu_(j+1)=([(j+1)-1]mu_((j+1)-1)+x_(j+1))/(j+1)
(6)
=(jmu_j+x_(j+1))/(j+1)
(7)
mu_(j+1)(j+1)=(j+1)mu_j+(x_(j+1)-mu_j)
(8)
mu_(j+1)=mu_j+(x_(j+1)-mu_j)/(j+1).
(9)

根據樣本方差的定義，

s_j^2=(sum_(i=1)^(j)(x_i-mu_j)^2)/(j-1)
(10)

對於 j>=2。定義 s_1=0，然後可以使用遞迴方程計算 s_j

js_(j+1)^2=j(sum_(i=1)^(j+1)(x_i-mu_(j+1))^2)/j
(11)
=sum_(i=1)^(j+1)(x_i-mu_(j+1))^2
(12)
=sum_(i=1)^(j+1)[(x_i-mu_j)+(mu_j-mu_(j+1))]^2
(13)
=sum_(i=1)^(j+1)(x_i-mu_j)^2+sum_(i=1)^(j+1)(mu_j-mu_(j+1))^2+2sum_(i=1)^(j+1)(x_i-mu_j)(mu_j-mu_(j+1)).
(14)

處理第一項，

sum_(i=1)^(j+1)(x_i-mu_j)^2=sum_(i=1)^(j)(x_i-mu_j)^2+(x_(j+1)-mu_j)^2
(15)
=(j-1)s_j^2+(x_(j+1)-mu_j)^2.
(16)

使用 (◇) 寫入

x_(j+1)-mu_j=(j+1)(mu_(j+1)-mu_j),
(17)

所以

sum_(i=1)^(j+1)(x_i-mu_j)^2=(j-1)s_j^2+(j+1)^2(mu_(j+1)-mu_j)^2.
(18)

現在處理 (◇) 中的第二項，

sum_(i=1)^(j+1)(mu_j-mu_(j+1))^2=(j+1)(mu_j-mu_(j+1))^2.
(19)

考慮 (◇) 中的第三項，

sum_(i=1)^(j+1)(x_i-mu_j)(mu_j-mu_(j+1))=(mu_j-mu_(j+1))sum_(i=1)^(j+1)(x_i-mu_j)
(20)
=(mu_j-mu_(j+1))[sum_(i=1)^(j)(x_i-mu_j)+(x_(j+1)-mu_j)]
(21)
=(mu_j-mu_(j+1))(x_(j+1)-mu_j-jmu_j+sum_(i=1)^(j)x_i).
(22)

但是

sum_(i=1)^jx_i=jmu_j,
(23)

所以

(mu_j-mu_(j+1))(x_(j+1)-mu_j)=(mu_j-mu_(j+1))(j+1)(mu_(j+1)-mu_j)
(24)
=-(j+1)(mu_j-mu_(j+1))^2.
(25)

最後，將 (◇)、(◇) 和 (◇) 代入 (◇)，

js_(j+1)^2=[(j-1)s_j^2+(j+1)^2(mu_(j+1)-mu_j)^2]+[(j+1)(mu_j-mu_(j+1))^2]+2[-(j+1)(mu_j-mu_(j+1))^2]
(26)
=(j-1)s_j^2+(j+1)^2(mu_(j+1)-mu_j)^2-(j+1)(mu_j-mu_(j+1))^2
(27)
=(j-1)s_j^2+(j+1)[(j+1)-1](mu_(j+1)-mu_j)^2
(28)
=(j-1)s_j^2+j(j+1)(mu_(j+1)-mu_j)^2,
(29)

給出了 s_(j+1) 關於 s_j、mu_(j+1) 和 mu_j 的所需表示式，

s_(j+1)^2=(1-1/j)s_j^2+(j+1)(mu_(j+1)-mu_j)^2.
(30)

另請參閱

樣本, 樣本方差, 樣本方差分佈, 方差

使用 探索

引用為

Weisstein, Eric W. “樣本方差計算。” 來自 Web 資源。 https://mathworld.tw/SampleVarianceComputation.html

主題分類