점추정량의 성질과 추정방법

yst3147·2022년 8월 13일

수리통계학

목록 보기
9/10

Wackerly et al., Mathematical Statistics with Applications을 읽고
책 내용을 정리하였습니다.

9.1 서론

9.2 상대 효율

상대적으로 효율적

θ^1,θ^2\hat{\theta}_1, \hat{\theta}_2이 같은 모수 θ\theta에 대한 두 개의 비편향 추정량이라면
V(θ^2)>V(θ^1)V(\hat{\theta}_2) > V(\hat{\theta}_1)이라면 θ^1\hat{\theta}_1θ^2\hat{\theta}_2보다 상대적으로 더 효율적(relatively more efficient)이다.

효율

분산이 각각 V(θ^1)V(\hat{\theta}_1), V(θ^2)V(\hat{\theta}_2)인 모수 θ\theta의 두 비편향 추정량 θ^1\hat{\theta}_1θ^2\hat{\theta}_2이 주어졌다면 θ^2\hat{\theta}_2에 대한 θ^1\hat{\theta}_1효율(efficiency), eff(θ^1,θ^2\hat{\theta}_1, \hat{\theta}_2)은
eff(θ^1,θ^2\hat{\theta}_1, \hat{\theta}_2) = V(θ^2)V(θ^1)\frac{V(\hat{\theta}_2)}{V(\hat{\theta}_1)}

9.3 일치성

수렴(convergence)

Y/nY/n은 확률변수이므로 확률용어로 pp접근(closeness)라 표현

확률수렴

n이 크다면 아래 식의 확률은 1에 가까워야만 한다
P(Ynpϵ)P(|\frac{Y}{n} - p| \leq \epsilon)
nn \to \infin일 때 이 확률이 1로 가면 Y/nY/npp일치 추정량(consistent estimator) 혹은 Y/nY/npp확률수렴(convergence in probability to p)한다고 한다.

일치추정량

만약 임의의 양수 ϵ\epsilon에 대하여
limnP(θn^θϵ)=1lim_{n \to \infin}P(|\hat{\theta_n} - \theta| \leq \epsilon) = 1
혹은
limnP(θn^θ>ϵ)=0lim_{n \to \infin}P(|\hat{\theta_n} - \theta| > \epsilon) = 0
이라면
추정량 θn^\hat{\theta_n}일치 추정량(consistent estimator)이라 한다.

비편향 추정량 -> 일치 추정량

θ\theta 에 대한 비편향 추정량 θn^\hat{\theta_n}은 만약
limnV(θn^)=0lim_{n \to \infin}V(\hat{\theta_n}) =0
라면 일치 추정량이다.

추정량의 일치성 증명을 위한 확률수렴 관련 정리

θn^\hat{\theta_n}θ\theta로 확률수렴하고 θn^\hat{\theta_n}'θ\theta'로 확률수렴한다고 가정
a. θn^+θn^\hat{\theta_n} + \hat{\theta_n}'θ+θ\theta + \theta'로 확률수렴
b. θn^×θn^\hat{\theta_n} \times \hat{\theta_n}'θ×θ\theta \times \theta'로 확률수렴
c. θ0\theta' \neq 0이라면 θn^/θn^\hat{\theta_n} / \hat{\theta_n}'θ/θ\theta / \theta'로 확률수렴
d. g()g(\cdot)θ\theta에서 연속인 실수 값 함수라면 g(θn^)g(\hat{\theta_n})g(θ)g(\theta)로 확률수렴

슬러츠키의 정리

UnU_nnn \to \infin에 따라 표준정규분포로 수렴하는 분포함수를 갖는다고 하자
만약 WnW_n이 1로 확률수렴한다면 Un/WnU_n/W_n의 분포는 표준정규분포함수로 수렴

9.4 충분성

충분 통계량

충분성을 갖는 통계량을 충분 통계량(sufficient statistics)이라고 한다.

충분

YY가 주어진 경우에 X1,X2,...XnX_1, X_2, ... X_n의 조건부 분포가 pp에 의존하지 않고
YYpp에 관한 모든 정보를 포함하면 YYpp에 대하여 충분(sufficient)라고 한다.

Y1,Y2,...YnY_1, Y_2, ... Y_n을 모수 θ\theta를 갖는 확률분포로부터의 확률표본이라 하면
통계량 U=g(Y1,Y2,...,Yn)U = g(Y_1, Y_2, ..., Y_n)UU가 주어진 경우에 Y1,Y2,...YnY_1, Y_2, ... Y_n의 조건부 분포가 θ\theta에 의존하지 않는다면 θ\theta에 대하여 충분하다고 한다.

가능도

y1,y2,...yny_1, y_2, ... y_n을 모수 θ\theta에 의존하는 대응하는 확률변수 Y1,Y2,...YnY_1, Y_2, ... Y_n의 표본관측값이라 하자
그러면 만약 Y1,Y2,...YnY_1, Y_2, ... Y_n가 이산확률변수라면 표본 가능도 L(y1,y2,...ynθ)L(y_1, y_2, ... y_n|\theta)y1,y2,...yny_1, y_2, ... y_n의 결합확률로 정의된다
Y1,Y2,...YnY_1, Y_2, ... Y_n가 연속확률변수라면 표본 가능도 L(y1,y2,...ynθ)L(y_1, y_2, ... y_n|\theta)y1,y2,...yny_1, y_2, ... y_n에서 평가된 결합밀도로 정의된다

분해정리

UU를 확률표본 Y1,Y2,...YnY_1, Y_2, ... Y_n에 기초한 통계량이라 한다. 그러면 UU가 모수 θ\theta의 추정에 대하여 충분 통계량일 필요충분조건은
가능도 L(θ)=L(y1,y2,...,ynθ)L(\theta) = L(y_1, y_2, ..., y_n | \theta)가 두 개의 음이 아닌 함수로 분해, 즉 g(u,θ)g(u, \theta)가 오직 uuθ\theta의 함수이고
h(y1,y2,..,yn)h(y_1, y_2, .., y_n)θ\theta의 함수가 아닐 때,
L(y1,y2,...,ynθ)=g(u,θ)×h(y1,y2,...,yn)L(y_1, y_2, ..., y_n|\theta) = g(u, \theta) \times h(y_1, y_2, ..., y_n)
이 된다.

9.5 라오-블랙웰 정리와 최소-분산비편향추정

라오-블랙웰 정리

θ^\hat{\theta}V(θ^)<V(\hat{\theta}) < \infinθ\theta에 대한 비편향 추정량이라 한다.
UUθ\theta에 대한 충분 통계량이라면
θ^=E(θ^U)\hat{\theta}^* = E(\hat{\theta} | U)라 정의한다.
그러면 모든 θ\theta에 대하여
E(θ^)=θ,V(θ^)V(θ^)E(\hat{\theta}^*) = \theta, V(\hat{\theta}^*) \leq V(\hat{\theta}) 이다.

최소충분통계량

모수 θ\theta에 대하여 자료에 있는 정보를 가장 잘 요약하는 통계량 UU최소충분통계량(minimal sufficient statistics)이라 한다.

최소-분산 비편향 추정량

UU를 사용하면 작은 분산을 갖는 비편향 추정량 뿐만 아니라 실제로 최소분산을 갖는 θ\theta에 대한 비편향 추정량을 얻을 수 있다는 것을 보장하는 다른 성질(완비성)을 가진 추정량

9.6 적률법

적률법

tt를 추정해야 하는 모수의 수라 할 때, 식 μk=mk,k=1,2,...t\mu_k' = m_k', k=1,2,... t의 해가 되는 모수의 값을 추정값으로 선택한다.

장점

  • 적용하기 쉽고 일치성을 갖는다
  • 사용하기 쉽고 때로는 합리적인 성질을 갖는 추정량을 가져온다

단점

  • 이 방법으로 유도된 추정량은 충분 통계량의 함수가 아닐 수 있다.
    -> 많은 경우에 적률법 추정량은 편향되어 있다.

9.7 최대가능도법

최대가능도법

가능도 함수가 kk개의 모수 θ1,θ2,...θk\theta_1, \theta_2, ... \theta_k에 의존한다고 가정
가능도 L(y1,y2,...,ynθ1,θ2,...,θk)L(y_1, y_2, ..., y_n| \theta_1, \theta_2, ..., \theta_k)
최대로 하는 모수의 값을 추정값으로 선택

장점

  • 일치성을 갖는다
  • 비편향이 되도록 조정된다면 최소분산 비편향추정량을 유도함

MLE의 불변성

t(θ)t(\theta)θ\theta의 일대일 함수이고 θ^\hat{\theta}θ\theta의 MLE라면 t(θ)t(\theta)의 MLE는
t(θ^)=t(θ^)t(\widehat{\theta}) = t(\hat{\theta})

9.8 최대가능도추정량의 대표본 성질(선택)

9.9 요약

0개의 댓글