적률추정법

choyunjeong·2024년 12월 12일

4.1 적률추정법

일반적으로 모적률 μr\mu_r'은 모수벡터 θ=(θ1,,θk)\theta=(\theta_1,\ldots,\theta_k)의 함수이므로 kk개의 연립방정식

mj=μj(θ1,,θk),j=1,,km_j'=\mu_j'(\theta_1,\ldots,\theta_k),\quad j=1,\ldots,k

를 풀면 (θ1^,,θk^)(\hat{\theta_1},\ldots,\hat{\theta_k})를 구할 수 있을 것이며, 이를 (θ1,,θk)(\theta_1,\ldots,\theta_k)에 대한 적률추정량이라고 한다.

이는 대수 법칙에 의하여 rr차 모적률이 존재할 때 rr차 표본적률이 rr차 모적률로 표본의 크기가 커짐에 따라 확률적으로 수렴한다는 사실에 기초한 방법으로서, nn이 크다면 표본적률이 확률적으로 모적률에 매우 가깝게 되므로 표본적률을 이용하여 모적률을 추정하는 것은 합리적인 방법이라고 할 수 있다.

참고
모수벡터 θ=(θ1,,θn)\theta=(\theta_1,\ldots,\theta_n)를 갖는 확률분포에서 X=(X1,,Xn)X=(X_1,\ldots,X_n)의 표본을 뽑는다고 하자.

모분포의 rr차 적률:

μr=E(Xr)\mu_r'=E(X^r)

rr차 표본적률:

mr=i=1n(Xr)/nm_r'=\sum_{i=1}^{n}(X^r)/n

\\[20pt]

예 4.2
X1,,XnX_1,\ldots,X_nN(μ,σ2)N(\mu,\sigma^2)으로부터 얻은 랜덤표본이라고 할 때, μ\muσ2\sigma^2을 적률을 사용하는 방법으로 추정.

정규분포의 모적률 값

  • 각 적률의 값은 적률생성함수로부터 구함.
    μ1=E(X1)=μμ2=E(X2)=μ2+σ2\mu_1'=E(X^1)=\mu \\[10pt] \mu_2'=E(X^2)=\mu^2+\sigma^2

표본적률

m1=i=1nXi1/nm2=i=1nXi2/nm_1'=\sum_{i=1}^{n}X_i^1/n \\[10pt] m_2'=\sum_{i=1}^{n}X_i^2/n

적률 추정

1) μ^\hat{\mu}추정

m1=μ1(μ,σ2)=μi=1nXi1/n=μ^μ^=Xˉn\begin{aligned} m_1'&=\mu_1'(\mu,\sigma^2)=\mu \\[10pt] \rightarrow \sum_{i=1}^{n}X_i^1/n &= \hat{\mu} \\[15pt] \therefore \hat{\mu} &= \bar{X}_n \end{aligned}

2) σ^2\hat{\sigma}^2추정

m2=μ2(μ,σ2)=μ2+σ2i=1nXi2/n=μ^2+σ^2σ^2=i=1nXi2/nμ^2=i=1nXi2/n(Xˉn)2( μ^=Xˉn)=i=1nXi2n(Xˉn)2n=i=1nXi2i=1n(Xˉn)2n=i=1n(XiXˉn)2n( E[X2][E[X]]2=E[XE(X)]2)\begin{aligned} m_2'&=\mu_2'(\mu,\sigma^2)=\mu^2+\sigma^2 \\[10pt] \rightarrow \sum_{i=1}^{n}X_i^2/n &= \hat{\mu}^2 + \hat{\sigma}^2\\[10pt] \therefore \hat{\sigma}^2 &= \sum_{i=1}^{n}X_i^2/n - \hat{\mu}^2 \\[15pt] &=\sum_{i=1}^{n}X_i^2/n-\left(\bar{X}_n\right)^2 \quad \left(\because\ \hat{\mu}=\bar{X}_n\right) \\[15pt] &=\dfrac{\sum_{i=1}^{n}X_i^2-n\left(\bar{X}_n\right)^2}{n} \\[10pt] &=\dfrac{\sum_{i=1}^{n}X_i^2-\sum_{i=1}^{n}\left(\bar{X}_n\right)^2}{n} \\[10pt] &=\dfrac{\sum_{i=1}^{n}\left(X_i-\bar{X}_n\right)^2}{n} \quad \left(\because\ E[X^2]-[E[X]]^2=E[X-E(X)]^2\right) \end{aligned}

\\[20pt]

예 4.3
X1,,XnX_1,\ldots,X_nGAM(k,θ)\text{GAM}(k,\theta)으로부터 얻은 랜덤표본이라고 할 때, μ\muσ2\sigma^2을 적률을 사용하는 방법으로 추정.

정규분포의 모적률 값

  • 각 적률의 값은 적률생성함수로부터 구함.
    μ1=E(X1)=kθμ2=E(X2)=kθ2+k2θ2\mu_1'=E(X^1)=k\theta\\[10pt] \mu_2'=E(X^2)=k\theta^2+k^2\theta^2

표본적률

m1=i=1nXi1/nm2=i=1nXi2/nm_1'=\sum_{i=1}^{n}X_i^1/n \\[10pt] m_2'=\sum_{i=1}^{n}X_i^2/n

적률 추정

1) k^\hat{k}추정

m1=μ1(k,θ)=kθi=1nXi1/n=k^θ^k^=Xˉnθ^\begin{aligned} m_1'&=\mu_1'(k,\theta)=k\theta \\[10pt] \rightarrow \sum_{i=1}^{n}X_i^1/n &= \hat{k}\hat{\theta} \\[15pt] \therefore \hat{k} &= \dfrac{\bar{X}_n}{\hat{\theta}} \end{aligned}

아래 추정된 θ^\hat{\theta}를 이용하여 k^\hat{k} 추정

k^=Xnˉ/i=1n(XiXnˉ)2nXnˉ=(nXnˉ)2i=1n(XiXnˉ)2\begin{aligned} \hat{k} &=\bar{X_n}/\dfrac{\sum_{i=1}^{n}(X_i-\bar{X_n})^2}{n\bar{X_n}} \\[15pt] &=\dfrac{(n\bar{X_n})^2}{\sum_{i=1}^{n}(X_i-\bar{X_n})^2} \end{aligned}

2) θ^\hat{\theta}추정

m2=μ2(k,θ)=kθ2+k2θ2=kθ(θ+kθ)i=1nXi2/n=k^θ^(θ^+k^θ^)=m1(θ^+m1)θ^=m2(m1)2m1=i=1nXi2/n(Xˉn)2Xˉn=i=1nXi2n(Xˉn)2nXˉn=i=1nXi2i=1n(Xˉn)2nXˉn=i=1n(XiXˉn)2nXˉn\begin{aligned} m_2'&=\mu_2'(k,\theta)= k\theta^2+k^2\theta^2 = k\theta(\theta+k\theta) \\[10pt] \rightarrow \sum_{i=1}^{n}X_i^2/n &= \hat{k}\hat{\theta}(\hat{\theta}+\hat{k}\hat{\theta}) = m_1'(\hat{\theta}+m_1') \\[10pt] \therefore \hat{\theta} &= \dfrac{m_2'-(m_1')^2}{m_1'} \\[10pt] &=\dfrac{\sum_{i=1}^{n}X_i^2/n-\left(\bar{X}_n\right)^2}{\bar{X}_n} \\[10pt] &=\dfrac{\sum_{i=1}^{n}X_i^2-n\left(\bar{X}_n\right)^2}{n\bar{X}_n} \\[10pt] &=\dfrac{\sum_{i=1}^{n}X_i^2-\sum_{i=1}^{n}\left(\bar{X}_n\right)^2}{n\bar{X}_n} \\[10pt] &=\dfrac{\sum_{i=1}^{n}(X_i-\bar{X}_n)^2}{n\bar{X}_n} \end{aligned}

[참고문헌]

  • 수리통계학 제 5판 - 송성주, 전명식

0개의 댓글