标准方差是方差无偏估计的证明——编辑版

1、问题描述:假设有一批独立同分布的样本{xi,i=1,2,3...,n−1,n{x_i,i=1,2,3...,n-1,n}xi,i=1,2,3...,n1,n}。
标准方差的公式为:
S2=1n−1∑i=1n(xi−x‾)2 S^{2} = \frac{1}{n-1}\sum_{i=1}^{n}(x_i - \overline{x})^{2} S2=n11i=1n(xix)2
,其中x‾=1n∑i=1nxi\overline{x}=\frac{1}{n}\sum_{i=1}^{n}x_ix=n1i=1nxi


另外假设均值E(xi)=μE(x_i) = \muE(xi)=μ,方差D(xi)=σ2D(x_i) = \sigma^{2}D(xi)=σ2


证明:标准方差是方差的无偏估计。

【分析】:要想证明标准方差是方差的无偏估计,只需证明E(S2)=σ2=D(xi)E(S^2) = \sigma^{2} = D(x_i)E(S2)=σ2=D(xi),其中比较麻烦的地方是xi∗x‾x_i * \overline{x}xix,所以要把x‾\overline{x}x中的xix_ixi分离处理。

2、证明开始:
E(S2)=E(1n−1∑i=1n(xi−x‾)2)=1n−1E(∑i=1n(xi−1nxi−1n∑j≠i,j=1nxj)2)=1n−1E(∑i=1n((n−1)2n2xi2−2(n−1)n2xi∗∑j≠i,j=1nxj+1n2(∑j≠i,j=1nxj)2)) \begin{aligned} E(S^{2}) &= E(\frac{1}{n-1}\sum_{i=1}^{n}(x_i - \overline{x})^{2}) \\ &= \frac{1}{n-1}E(\sum_{i=1}^{n}(x_i - \frac{1}{n}x_i - \frac{1}{n}\sum_{j\neq{i},j=1}^{n}x_j)^{2})\\ &= \frac{1}{n-1}E(\sum_{i=1}^{n}(\frac{(n-1)^2}{n^2}x_i^{2} - \frac{2(n-1)}{n^{2}}x_i *\sum_{j\neq{i},j=1}^{n}x_j + \frac{1}{n^2}(\sum_{j\neq{i},j=1}^{n}x_j)^{2})) \\ \end{aligned} E(S2)=E(n11i=1n(xix)2)=n11E(i=1n(xin1xin1j=i,j=1nxj)2)=n11E(i=1n(n2(n1)2xi2n22(n1)xij=i,j=1nxj+n21(j=i,j=1nxj)2))
下面将E(⋅)E(\centerdot)E()中的三部分分别计算,按照从左往右的顺序,系数1n−1\frac{1}{n-1}n11暂时不带,最后算完结果一起带入。
第一部分:
E(∑i=1n((n−1)2n2)xi2)=(n−1)2n2∑i=1nE(xi2)=(n−1)2n2∑i=1n(E2(xi)+D(xi))=(n−1)2n2∑i=1n(μ2+σ2)=(n−1)2n(μ2+σ2) \begin{aligned} E(\sum_{i=1}^{n}(\frac{(n-1)^2}{n^2})x_i^{2}) &= \frac{(n-1)^2}{n^2}\sum_{i=1}^{n}E(x_i^{2}) \\ &= \frac{(n-1)^2}{n^2}\sum_{i=1}^{n}(E^2(x_i) + D(x_i)) \\ &= \frac{(n-1)^2}{n^2}\sum_{i=1}^{n}(\mu^{2} + \sigma^{2}) \\ &= \frac{(n-1)^2}{n}(\mu^{2} + \sigma^{2}) \end{aligned} E(i=1n(n2(n1)2)xi2)=n2(n1)2i=1nE(xi2)=n2(n1)2i=1n(E2(xi)+D(xi))=n2(n1)2i=1n(μ2+σ2)=n(n1)2(μ2+σ2)
第二部分:
E(∑i=1n(2(n−1)n2xi∗∑j≠i,j=1nxj))=2(n−1)n2∑i=1nE(xi∗∑j≠i,j=1nxj))=2(n−1)n2∑i=1n(E(xi)∗∑j≠i,j=1nE(xj))=2(n−1)n2∑i=1n(μ∗(n−1)μ)=2(n−1)n2∗n∗(n−1)∗μ2=2(n−1)2∗μ2n \begin{aligned} E(\sum_{i=1}^{n}( \frac{2(n-1)}{n^{2}}x_i *\sum_{j\neq{i},j=1}^{n}x_j)) &= \frac{2(n-1)}{n^{2}}\sum_{i=1}^{n}E(x_i *\sum_{j\neq{i},j=1}^{n}x_j)) \\ &= \frac{2(n-1)}{n^{2}}\sum_{i=1}^{n}(E(x_i) * \sum_{j\neq{i},j=1}^{n}E(x_j)) \\ &= \frac{2(n-1)}{n^{2}}\sum_{i=1}^{n}(\mu * (n-1)\mu) \\ &= \frac{2(n-1)}{n^{2}} * n*(n-1)*\mu^{2} \\ & = \frac{2(n-1)^{2}*\mu^{2}}{n} \end{aligned} E(i=1n(n22(n1)xij=i,j=1nxj))=n22(n1)i=1nE(xij=i,j=1nxj))=n22(n1)i=1n(E(xi)j=i,j=1nE(xj))=n22(n1)i=1n(μ(n1)μ)=n22(n1)n(n1)μ2=n2(n1)2μ2
第三部分:
E(∑i=1n(1n2(∑j≠i,j=1nxj)2))=1n2∑i=1nE(((∑j≠i,j=1nxj)2))=1n2∑i=1n(E2(∑j≠i,j=1nxj)+D(∑j≠i,j=1nxj))=1n2∑i=1n((∑j≠i,j=1nE2(xj))+∑j≠i,j=1nD(xj))=1n2∑i=1n(((∑j≠i,j=1nμ)2)+∑j≠i,j=1nσ2)=1n2∑i=1n((n−1)2∗μ2+(n−1)∗σ2)=(n−1)2n2∗n∗μ2+n−1n2∗nσ2=(n−1)2n∗μ2+n−1nσ2 \begin{aligned} E(\sum_{i=1}^{n}(\frac{1}{n^2}(\sum_{j\neq{i},j=1}^{n}x_j)^{2})) &= \frac{1}{n^2}\sum_{i=1}^{n}E(((\sum_{j\neq{i},j=1}^{n}x_j)^{2})) \\ &= \frac{1}{n^2}\sum_{i=1}^{n}(E^{2}(\sum_{j\neq{i},j=1}^{n}x_j) +D(\sum_{j\neq{i},j=1}^{n}x_j)) \\ &= \frac{1}{n^2}\sum_{i=1}^{n}((\sum_{j\neq{i},j=1}^{n}E^{2}(x_j))+\sum_{j\neq{i},j=1}^{n}D(x_j)) \\ &= \frac{1}{n^2}\sum_{i=1}^{n}(((\sum_{j\neq{i},j=1}^{n}\mu)^{2})+\sum_{j\neq{i},j=1}^{n}\sigma^{2}) \\ &= \frac{1}{n^2}\sum_{i=1}^{n}((n-1)^{2}*\mu^2+(n-1)*\sigma^2) \\ &= \frac{(n-1)^2}{n^2} *n *\mu^2+\frac{n-1}{n^2} *n\sigma^2\\ &= \frac{(n-1)^2}{n} *\mu^2+\frac{n-1}{n} \sigma^2 \end{aligned} E(i=1n(n21(j=i,j=1nxj)2))=n21i=1nE(((j=i,j=1nxj)2))=n21i=1n(E2(j=i,j=1nxj)+D(j=i,j=1nxj))=n21i=1n((j=i,j=1nE2(xj))+j=i,j=1nD(xj))=n21i=1n(((j=i,j=1nμ)2)+j=i,j=1nσ2)=n21i=1n((n1)2μ2+(n1)σ2)=n2(n1)2nμ2+n2n1nσ2=n(n1)2μ2+nn1σ2
终于到了相加的时候了,将1-3部分的内容带上系数相加如下:
(n−1)E(⋅)=(n−1)2n(μ2+σ2) −2(n−1)2∗μ2n+(n−1)2n∗μ2+n−1nσ2=(n−1)2nσ2+n−1nσ2=n2−nnσ2=(n−1)σ2 \begin{aligned} (n-1)E(\centerdot) & =\frac{(n-1)^2}{n}(\mu^{2} + \sigma^{2}) \ - \frac{2(n-1)^{2}*\mu^{2}}{n} + \frac{(n-1)^2}{n} *\mu^2+\frac{n-1}{n} \sigma^2 \\ &= \frac{(n-1)^2}{n}\sigma^{2} + \frac{n-1}{n} \sigma^2 \\ &= \frac{n^2-n}{n}\sigma^{2} \\ &= (n-1)\sigma^{2} \end{aligned} (n1)E()=n(n1)2(μ2+σ2) n2(n1)2μ2+n(n1)2μ2+nn1σ2=n(n1)2σ2+nn1σ2=nn2nσ2=(n1)σ2
所以
E(⋅)=σ2 E(\centerdot) = \sigma^2 E()=σ2
得证。
注:标准方差前面的系数1n−1\frac{1}{n-1}n11是为了让其是方差的无偏估计而加上的。

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值