4.1 적률추정법
일반적으로 모적률 μr′은 모수벡터 θ=(θ1,…,θk)의 함수이므로 k개의 연립방정식
mj′=μj′(θ1,…,θk),j=1,…,k
를 풀면 (θ1^,…,θk^)를 구할 수 있을 것이며, 이를 (θ1,…,θk)에 대한 적률추정량이라고 한다.
이는 대수 법칙에 의하여 r차 모적률이 존재할 때 r차 표본적률이 r차 모적률로 표본의 크기가 커짐에 따라 확률적으로 수렴한다는 사실에 기초한 방법으로서, n이 크다면 표본적률이 확률적으로 모적률에 매우 가깝게 되므로 표본적률을 이용하여 모적률을 추정하는 것은 합리적인 방법이라고 할 수 있다.
참고
모수벡터 θ=(θ1,…,θn)를 갖는 확률분포에서 X=(X1,…,Xn)의 표본을 뽑는다고 하자.
모분포의 r차 적률:
μr′=E(Xr)
r차 표본적률:
mr′=i=1∑n(Xr)/n
예 4.2
X1,…,Xn을 N(μ,σ2)으로부터 얻은 랜덤표본이라고 할 때, μ와 σ2을 적률을 사용하는 방법으로 추정.
정규분포의 모적률 값
- 각 적률의 값은 적률생성함수로부터 구함.
μ1′=E(X1)=μμ2′=E(X2)=μ2+σ2
표본적률
m1′=i=1∑nXi1/nm2′=i=1∑nXi2/n
적률 추정
1) μ^추정
m1′→i=1∑nXi1/n∴μ^=μ1′(μ,σ2)=μ=μ^=Xˉn
2) σ^2추정
m2′→i=1∑nXi2/n∴σ^2=μ2′(μ,σ2)=μ2+σ2=μ^2+σ^2=i=1∑nXi2/n−μ^2=i=1∑nXi2/n−(Xˉn)2(∵ μ^=Xˉn)=n∑i=1nXi2−n(Xˉn)2=n∑i=1nXi2−∑i=1n(Xˉn)2=n∑i=1n(Xi−Xˉn)2(∵ E[X2]−[E[X]]2=E[X−E(X)]2)
예 4.3
X1,…,Xn을 GAM(k,θ)으로부터 얻은 랜덤표본이라고 할 때, μ와 σ2을 적률을 사용하는 방법으로 추정.
정규분포의 모적률 값
- 각 적률의 값은 적률생성함수로부터 구함.
μ1′=E(X1)=kθμ2′=E(X2)=kθ2+k2θ2
표본적률
m1′=i=1∑nXi1/nm2′=i=1∑nXi2/n
적률 추정
1) k^추정
m1′→i=1∑nXi1/n∴k^=μ1′(k,θ)=kθ=k^θ^=θ^Xˉn
아래 추정된 θ^를 이용하여 k^ 추정
k^=Xnˉ/nXnˉ∑i=1n(Xi−Xnˉ)2=∑i=1n(Xi−Xnˉ)2(nXnˉ)2
2) θ^추정
m2′→i=1∑nXi2/n∴θ^=μ2′(k,θ)=kθ2+k2θ2=kθ(θ+kθ)=k^θ^(θ^+k^θ^)=m1′(θ^+m1′)=m1′m2′−(m1′)2=Xˉn∑i=1nXi2/n−(Xˉn)2=nXˉn∑i=1nXi2−n(Xˉn)2=nXˉn∑i=1nXi2−∑i=1n(Xˉn)2=nXˉn∑i=1n(Xi−Xˉn)2
[참고문헌]