5️⃣

5장 Consistency and Limiting Distributions

왜 저희는 5장을 배워야 하나요?

모집단과 표본이 있다면 저희는 표본을 이용해서 모집단을 추정하는 것이 목적입니다.

정리하자면 모수 θ\theta 는 모집단 내에서 저희가 궁금한 정보이며 (예: 모평균, 모분산 등등)

이 값을 구하기 위해서 모집단 내에 있는 모든 값들을 일일이 비교하는 것은 불가능하기 때문에

표본 (Sample) 에서 뽑아서 표본 X1,⋯ ,XnX_1, \cdots, X_n 을 이용해서 모수 θ\theta 를 추정 및 검정을 하는 겁니다.

그렇다면 이 표본을 통해서 구체적으로 어떻게 하면 모수를 추정할 수 있을까요? 바로 통계량 (Statistic) 을

만드는 것입니다! 그럼 통계량은 무엇일까요? 모수를 나타내기 위해서 표본들을 이용한 함수라고 생각하시면 됩니다.

Tn=u(X1,X2,⋯ ,Xn)T_n = u(X_1, X_2, \cdots, X_n)

근데 이 통계량 자체가 어떠한 분포를 따라가는지, 즉 어떤 극한 분포 (Limiting Distribution) 을 가지는지를 알아야만 모수에 대한 추정 + 검정 이 가능합니다!

Convergence in Probability
Definition 정의
{Xn}: Sequence of Random VariablesX: Random Variable defined on a sample space We say that Xnconverges in Probability to X,Xn→PX, if, for all ϵ>0,lim⁡n→∞P[∣Xn−X∣≥ϵ]=0orlim⁡n→∞P[∣Xn−X∣<ϵ]=1\{X_n\} : \text{ Sequence of Random Variables}\\ X : \text{ Random Variable defined on a sample space}\\~\\ \text{We say that } X_n \text{converges in Probability to } X, X_n \xrightarrow{P} X, \text{ if, for all }\epsilon >0,\\ \lim _{n \rightarrow \infty} P\left[\left|X_n-X\right| \geq \epsilon\right]=0 \\ \text{or} \\ \lim _{n \rightarrow \infty} P\left[\left|X_n-X\right|<\epsilon\right]=1

Example of Convergence in Probability →\rightarrow Weak Law of Large Numbers (WLLN)

Let {Xn}\{X_n\} be a sequence of i.i.d random variables, with common mean μ\mu and variance σ2<∞\sigma^2 < \infty

Then the following is True

Xˉn=1n∑i=1nXi Xˉn→Pμ\bar{X}_n = \frac{1}{n} \sum_{i=1}^{n} X_i \\~\\ \bar{X}_n \xrightarrow{P} \mu

증명 :

It is given already that E⁡[Xˉn]=μ,Var⁡[Xˉn]=σ2n\operatorname{E}\left[\bar{X}_n \right] = \mu, \quad \operatorname{Var}\left[\bar{X}_n \right] = \frac{\sigma^2}{n}, then by using Chebyshev’s Inequality,

Chebyshev’s Inequality :

∀a>0,P⁡[∣X−μ∣≥a⋅σ]≤1a2\forall a >0, \quad\operatorname{P}\left[|X - \mu| \ge a \cdot \sigma \right] \le \frac{1}{a^2}

we can say the following is True

∀ϵ>0,P⁡[∣Xˉn−μ∣≥ϵ]=P⁡[∣Xˉn−μ∣≥(ϵn/σ)⋅(σ/n)]=P⁡[∣Xˉn−μ∣≥a⋅(σ/n)]≤1a2=1(ϵnσ)2=σ2n⋅ϵ2→n→∞0 Hence, P⁡[∣Xˉn−μ∣≥ϵ]=0,Xn Converges in Probability to μ!\forall \epsilon > 0, \quad \operatorname{P}\left[|\bar{X}_n - \mu| \ge \epsilon \right] = \operatorname{P}\left[|\bar{X}_n - \mu| \ge \left(\epsilon \sqrt{n} / \sigma \right) \cdot \left(\sigma / \sqrt{n} \right) \right] \\ = \operatorname{P}\left[|\bar{X}_n - \mu| \ge a \cdot \left(\sigma / \sqrt{n}\right) \right] \le \frac{1}{a^2} = \frac{1}{\left(\frac{\epsilon \sqrt{n}}{\sigma} \right)^2} = \frac{\sigma^2}{n \cdot \epsilon^2} \xrightarrow{n \rightarrow \infty} 0 \\~\\ \text{Hence, } \operatorname{P}\left[|\bar{X}_n - \mu| \ge \epsilon \right] = 0, \quad X_n \text{ Converges in Probability to } \mu !
Xn+Yn→PX+YX_n + Y_n \xrightarrow{P} X + Y

Suppose that Xn→PXX_n \xrightarrow{P} X, and Yn→PYY_n \xrightarrow{P} Y, then the following is True

Xn+Yn→PX+YX_n + Y_n \xrightarrow{P} X + Y

증명 :

Let ∀ϵ>0 be True, Triangle Inequality is as follows ∣a+b∣≤∣a∣+∣b∣\text{Let } \forall \epsilon >0 \text{ be True, Triangle Inequality is as follows } \\ |a + b| \le |a| + |b|

그럼 다음과 같이 대입을 시켜서 정리가 됩니다.

ϵ≤∣(Xn+Yn)−(X+Y)∣≤∣Xn−X∣+∣Yn−Y∣\epsilon \le |(X_n + Y_n) - (X + Y)| \le |X_n - X| + |Y_n - Y|

Convergence in Probability 정의를 가져오면 다음과 같습니다.

P⁡(∣Xn+Yn−X−Y∣≥ϵ)≤P⁡(∣Xn−X∣+∣Yn−Y∣≥ϵ)≤P⁡(∣Xn−X∣≥ϵ/2)+P⁡(∣Yn−Y∣≥ϵ/2)\operatorname{P}\left(|X_n + Y_n - X - Y| \ge \epsilon \right) \le \operatorname{P}\left(|X_n - X| + |Y_n - Y| \ge \epsilon \right) \le \\ \operatorname{P}\left(|X_n - X| \ge \epsilon/2\right) + \operatorname{P}\left(|Y_n - Y| \ge \epsilon/2 \right)

하지만 저희는 이미 Xn→PXX_n \xrightarrow{P} X 와 Yn→PYY_n \xrightarrow{P} Y 이 사실이라는 것을 알고 있습니다. 즉

∀ϵ>0,P⁡(∣Xn−X∣≥ϵ/2)=0 and P⁡(∣Yn−Y∣≥ϵ/2)=0\forall \epsilon >0, \quad \operatorname{P}\left(|X_n - X| \ge \epsilon/2\right) = 0 \text{ and } \operatorname{P}\left(|Y_n - Y| \ge \epsilon/2\right) = 0

그렇기 때문에

P⁡(∣Xn−X∣≥ϵ/2)+P⁡(∣Yn−Y∣≥ϵ/2)=0≥P⁡(∣Xn−X∣+∣Yn−Y∣≥ϵ)≥P⁡(∣Xn+Yn−X−Y∣≥ϵ) ∴Xn+Yn→PX+Y\operatorname{P}\left(|X_n - X| \ge \epsilon/2\right) + \operatorname{P}\left(|Y_n - Y| \ge \epsilon/2 \right) = 0 \\ \ge \operatorname{P}\left(|X_n - X| + |Y_n - Y| \ge \epsilon \right) \\ \ge \operatorname{P}\left(|X_n + Y_n - X - Y| \ge \epsilon \right) \\~\\ \therefore X_n + Y_n \xrightarrow{P} X + Y
a⋅Xn→Pa⋅Xa \cdot X_n \xrightarrow{P} a \cdot X

Suppose that aa is a constant, and that Xn→PXX_n \xrightarrow{P} X, then the following is True

a⋅Xn→Pa⋅Xa \cdot X_n \xrightarrow{P} a \cdot X

증명:

a=0a = 0 인 경우와 a≠0a \ne 0 인 경우들을 살펴야 합니다.

전자의 경우 식에 대입만 시키면 됩니다.

후자의 경우 식에 대입시키면 다음과 같이 Convergence in Probability 의 정의를 활용해서 다음과 같이 식으로 정리가 가능합니다.

P[∣aXn−aX∣≥ϵ]=P[∣a∣∣Xn−X∣≥ϵ]P\left[\left|a X_n-a X\right| \geq \epsilon\right] = P\left[|a|\left|X_n-X\right| \geq \epsilon\right]

aa 를 앞으로 갖고오는 것은 당연히 가능합니다! 양쪽에다가 나눠주기만 한다면 다음과 같이 나오게 됩니다.

P[∣a∣∣Xn−X∣≥ϵ]=P[∣Xn−X∣≥ϵ/∣a∣]P\left[|a|\left|X_n-X\right| \geq \epsilon\right] = P\left[\left|X_n-X\right| \geq \epsilon /|a|\right]

우측 식은 사실상 Xn→PXX_n \xrightarrow{P} X 를 의미하고 있기 때문에 해당 식은 True 이며 증명을 마치겠습니다.

Continuous Mapping Theorem

Suppose that Xn→PaX_n \xrightarrow{P} a and gg is a continuous function at aa, then the following is True

g(Xn)→Pg(a)g(X_n) \xrightarrow{P} g(a)

증명 :

Let ∀ϵ>0\forall \epsilon >0, we know that gg is continuous at aa, 그리고 이것이 의미하는 것은 바로 다른 ∀δ>0\forall \delta >0 과 활용해서 표현하면 다음과 같이 표현이 가능합니다.

if ∣x−a∣<δ,then ∣g(x)−g(a)∣<ϵ if ∣g(x)−g(a)∣≥ϵ,then ∣x−a∣≥δ\text{if } |x-a| < \delta, \quad \text{then } |g(x) - g(a)| < \epsilon \\~\\ \text{if } |g(x) - g(a)| \ge \epsilon, \quad \text{then } |x-a| \ge \delta

이것이 의미하는 것은 다음과 같습니다.

이것을 저희의 식에 적용시키게 된다면 다음과 같이 나오게 됩니다.

P⁡[∣g(Xn)−g(a)∣≥ϵ]≤P⁡[∣Xn−a∣≥δ]=0\operatorname{P}\left[|g(X_n) - g(a)| \ge \epsilon \right] \le \operatorname{P}\left[|X_n - a| \ge \delta \right] = 0

우측 식의 경우 Xn→PaX_n \xrightarrow{P} a 에 의해서 0 의 값을 가지기 때문에 좌측 식도 0으로 되면서 Convergence in Probability 가 성립이 됩니다.

여기에서 그럼 Inequality 가 왜 해당이 될까요?

이것이 해당 되는 이유는 어차피 δ\delta 랑 ϵ\epsilon 은 Arbitrary 한 값으로 저것을 만족하는 값인 Arbitrary 값을 가져오면 됩니다. Convergence in Probability 를 증명하기 위해서 저렇게 사용됐다~ 라고 생각하시면 됩니다!

그럼 이것이 사용될 만한 예시들을 보여드리겠습니다.

추가적으로 이것보다 더 어려운 강의들에서는 다음이 만족 된다는 것을 증명하며 해당 강의에서는 그저 이것이 해당이 되구나~ 로만 이해해도 무방합니다. (위와는 다르게 이번에는 확률변수입니다)

If Xn→PX, and g is a continuous function, then the following is Trueg(Xn)→Pg(X)\text{If } X_n \xrightarrow{P} X, \\\text{ and } g \text{ is a continuous function, then the following is True} \\ g(X_n) \xrightarrow{P} g(X)
Xn⋅Yn→PXYX_n \cdot Y_n \xrightarrow{P} XY

이전의 증명들을 이용하면 다음과 같이 식을 사용할 수 있습니다.

Xn⋅Yn=12Xn2+12Yn2−12(Xn−Yn)2X_n \cdot Y_n = \frac{1}{2} X_n^2 + \frac{1}{2}Y_n^2 - \frac{1}{2} \left(X_n - Y_n \right)^2

이거는 그저 xy 를 (x-y)^2 와 연관 지어서 푼 것이며, 이것이 가능한 이유는 다음과 같습니다.

xy=12x2+12y2−12(x−y)2=12x2+12y2−12(x2−2xy+y2)=12x2+12y2−12x2+xy−12y2=12x2−12x2+12y2−12y2+xy=xyxy = \frac{1}{2}x^2 + \frac{1}{2}y^2 - \frac{1}{2}(x - y)^2 = \frac{1}{2}x^2 + \frac{1}{2}y^2 - \frac{1}{2}\left(x^2 - 2xy + y^2 \right) \\ = \frac{1}{2}x^2 + \frac{1}{2}y^2 - \frac{1}{2}x^2 + xy - \frac{1}{2} y^2 = \frac{1}{2}x^2 -\frac{1}{2}x^2 +\frac{1}{2}y^2 -\frac{1}{2}y^2 + xy = xy

그럼 기존 식을 다시 살피겠습니다. 해당 식에서 상수를 빼서 유심히 살피면 Xn2,Yn2,(Xn−Yn)2X_n^2, Y_n^2, (X_n - Y_n)^2 요소들이 있습니다.

Xn2X_n^2 와 Yn2Y_n^2 의 경우 :

gg 를 g(x)=x2g(x) = x^2 으로 잡는다면 Continuous Mapping Theorem 에 의해서 Convergence in Probability 가 해당이 된다는 것을 알 수가 있습니다.

(Xn−Yn)2\left(X_n - Y_n \right)^2 의 경우 :

우선은 Xn−Yn→PX−YX_n - Y_n \xrightarrow{P} X - Y 이라는 것을 알 수가 있습니다. 이것 자체를 표현을 조금 바꾸면 다음과 같습니다 : Zn→PZZ_n \xrightarrow{P} Z

gg 를 g(x)=x2g(x) = x^2 으로 잡는다면 Continuous Mapping Theorem 에 의해서 다음이 해당이 됩니다.

g(Zn)→Pg(Z),↔(Zn)2→P(Z)2,↔(Xn−Yn)2→P(X−Y)2g(Z_n) \xrightarrow{P} g(Z), \leftrightarrow \left(Z_n\right)^2 \xrightarrow{P} \left(Z \right)^2, \leftrightarrow \left(X_n - Y_n \right)^2 \xrightarrow{P} \left(X - Y \right)^2

정리를 하면 다음이 해당이 됩니다.

Xn⋅Yn=12Xn2+12Yn2−12(Xn−Yn)2→P12X2+12Y2−12(X−Y)2=X⋅YX_n \cdot Y_n = \frac{1}{2} X_n^2 + \frac{1}{2}Y_n^2 - \frac{1}{2} \left(X_n - Y_n \right)^2 \xrightarrow{P} \frac{1}{2}X^2 + \frac{1}{2}Y^2 - \frac{1}{2}\left(X - Y\right)^2 = X \cdot Y

중간의 Convergence in Probability 가 성립이 되기에 증명을 마치겠습니다.

Consistency

우선은 Unbiasedness 에 대해서 복습을 하겠습니다.

Unbiasedness

Statistic, 통계량, TT 가 있을 때 만약에 E⁡[T]=θ\operatorname{E}[T] = \theta 가 만족이 되면 Unbiased Estimator 라고 합니다.

Consistency

Statistic, 통계량, TT 가 있을 때 만약에 Tn→PθT_n \xrightarrow{P} \theta 가 만족이 되면 Consistent Estimator 라고 합니다.

이것이 의미하는 것은 표본 개수가 오르면 오를수록 모수랑 같아진다는 의미를 가집니다.

Sample Variance

Sn2=1n−1∑i=1n(Xi−Xˉn)2=nn−1(∑i=1nXi2−Xˉn2)→P1⋅[E⁡(X12)−μ2]=Var⁡(X1)=σ2S^2_n = \frac{1}{n-1} \sum_{i=1}^{n} \left(X_i - \bar{X}_n \right)^2 = \frac{n}{n-1} \left(\sum_{i=1}^{n}X_i^2 - \bar{X}_n^2 \right) \\ \xrightarrow{P} 1 \cdot \left[\operatorname{E}\left(X_1^2 \right) - \mu^2 \right] = \operatorname{Var}(X_1) = \sigma^2

Hence Sample Variance is a Consistent Estimator

Convergence in Distribution
Definition 정의
{Xn}: Sequence of Random VariablesX: Random Variable defined on a sample spaceFXn: cdf of XnFX: cdf of XC(FX): The set of all points where FX is continuous We say that Xnconverges in Distribution to X,Xn→DX, if, for all x∈C(FX), lim⁡n→∞FXn(x)=FX(x)\{X_n\} : \text{ Sequence of Random Variables}\\ X : \text{ Random Variable defined on a sample space}\\ F_{X_n} : \text{ cdf of } X_n\\ F_{X} : \text{ cdf of } X\\ C(F_X) : \text{ The set of all points where } F_X \text{ is continuous} \\~\\ \text{We say that } X_n \text{converges in Distribution to } X, X_n \xrightarrow{D} X, \\ \text{ if, for all } x \in C(F_X),\\~\\ \lim_{n \rightarrow \infty} F_{X_n}(x) = F_X(x)

사실상 다른게 아니라 XnX_n 과 XX 의 CDF 가 같으면 그냥 그것을 Convergence in Distribution 이라고 보는 것입니다.

그리고 이게 의미하는게 무엇인지 헷갈릴 수도 있는데, Convergence in Distribution 이

Xn→DXX_n \xrightarrow{D} X 과 같이 있다면 해당 Sequence 는 XX 의 분포를 따라간다는 뜻이며, 이것의 예시를 갖고오면 이해가 편할 것으로 생각되어 다음과 같이 쓰겠습니다.

예:

Xn→DN⁡(0,1)X_n \xrightarrow{D} \operatorname{N}(0, 1)

여기에서 저희는 표현을 할 때에 XX 는 Sequence {Xn}\{X_n\} 의 Limiting Distribution / Asymptotic Distribution 이라고 부르기도 합니다.

Convergence in Probability and Distribution

If Xn→PXX_n \xrightarrow{P} X, then Xn→DXX_n \xrightarrow{D} X is True.

즉 Convergence in Probability 하다면 Convergence in Distribution 도 해당이 된다는 뜻입니다.

이것은 if 관계를 가지며, 반대는 항상 성립이 되지는 않기 때문에 if and only if 관계가 성립이 되지는 않습니다.

증명 :

Upper Bound 과 Lower Bound 를 다음과 같이 구할 수 있습니다.

∀ϵ>0,FXn(x)=P⁡[Xn≤x]=P⁡[{Xn≤x}∩{∣Xn−X∣<ϵ}]+P⁡[{Xn≤x}∩{∣Xn−X∣≥ϵ}]≤P⁡[X≤x+ϵ]+P⁡[∣Xn−X∣≥ϵ]\forall \epsilon >0, \quad F_{X_n}(x) = \operatorname{P}\left[X_n \le x \right] \\ = \operatorname{P}\left[\{X_n \le x\} \cap \{|X_n - X| < \epsilon\}\right]+ \operatorname{P}\left[\{X_n \le x\} \cap \{|X_n - X| \ge \epsilon\} \right] \\ \le \operatorname{P}\left[X \le x + \epsilon \right] + \operatorname{P}\left[|X_n - X| \ge \epsilon \right]

Convergence in Probability 가 해당이 되기에 마지막 식의 오른쪽 부분은 0 으로 수렴하며 좌측 부분만 이용해서 Upper Bound 을 다음과 같이 정리할 수 있습니다.

lim‾⁡n→∞FXn(x)=P⁡[Xn≤x]≤P⁡[X≤x+ϵ]=FX(x+ϵ)lim‾⁡n→∞FXn(x)≤FX(x+ϵ)\varlimsup_{n \rightarrow \infty} F_{X_n}(x) = \operatorname{P}\left[X_n \le x \right] \le \operatorname{P}\left[X \le x + \epsilon \right] = F_X(x + \epsilon) \\ \varlimsup_{n \rightarrow \infty}F_{X_n}(x) \le F_X(x+\epsilon)

Lower Bound 는 똑같이 위에서 했던거랑 하면 되며, 이번에는 반대로 다음과 같이 구하겠습니다.

1−FXn(x)=P⁡[Xn>x]=P⁡[{Xn>x}∩{∣Xn−X∣<ϵ}]+P⁡[{Xn>x}∩{∣Xn−X∣≥ϵ}]≥P⁡[X≥x−ϵ]+P⁡[∣Xn−X∣≥ϵ]=P⁡[X≥x−ϵ]1- F_{X_n}(x) = \operatorname{P}\left[X_n > x\right] \\ = \operatorname{P}\left[\{X_n > x\} \cap \{|X_n - X| < \epsilon\}\right]+ \operatorname{P}\left[\{X_n > x\} \cap \{|X_n - X| \ge \epsilon\} \right] \\ \ge \operatorname{P}\left[X \ge x - \epsilon \right] + \operatorname{P}\left[|X_n - X| \ge \epsilon \right] = \operatorname{P}\left[X \ge x-\epsilon \right]

위의 식을 정리하면 Limit 에 대한 Lower Bound 도 구할 수가 있습니다.

lim‾⁡n→∞FXn(x)≥FX(x−ϵ)\varliminf_{n \rightarrow \infty} F_{X_n}(x) \geq F_X(x-\epsilon)

두 식을 합치면 다음과 같이 Sandwich Theorem 과 ϵ\epsilon 값을 아주 작게 구한다면 다음과 같은 결론을 지을 수가 있습니다.

FX(x−ϵ)≤lim‾⁡n→∞FXn(x)≤lim‾⁡n→∞FXn(x)≤FX(x+ϵ) Xn→DXF_X(x-\epsilon) \leq \varliminf_{n \rightarrow \infty} F_{X_n}(x) \leq \varlimsup_{n \rightarrow \infty} F_{X_n}(x) \leq F_X(x+\epsilon) \\~\\ X_n \xrightarrow{D} X
Xn→DbX_n \xrightarrow{D} b

If Xn→DbX_n \xrightarrow{D} b, where bb is a constant, then Xn→PbX_n \xrightarrow{P} b

Let ∀ϵ>0\forall \epsilon > 0, then

lim⁡n→∞P⁡[∣Xn−b∣≤ϵ]=lim⁡n→∞FXn(b+ϵ)−lim⁡n→∞FXn(b−ϵ)=1−0=1\lim_{n \rightarrow \infty}\operatorname{P}\left[|X_n - b| \le \epsilon \right] = \lim_{n\rightarrow \infty} F_{X_n}(b+\epsilon) - \lim_{n \rightarrow \infty} F_{X_n}(b-\epsilon) = 1 - 0 = 1

설명을 해보겠습니다. Convergence in Distribution 이기 때문에 다음과 같은 식이 성립이 됩니다.

lim⁡n→∞FXn(t)=Fb(t)\lim _{n \rightarrow \infty} F_{X_n}(t)=F_b(t)

만약에 t<bt < b 이면 Fb(t)=0F_b(t) = 0, 아니면 Fb(t)=1F_b(t) = 1 입니다. 이 이유는 CDF 의 성질을 생각해보면 됩니다. Constant 이기 때문에 bb 시점을 지나면 1, 아니면 0 인 CDF 를 가지겠죠.

그러면 t<bt <b 와 t≥bt\ge b 으로 나눌 수가 있으며, bb 시점에 대해서는 우측 Limit 와 좌측 Limit 를 구할 수가 있게 됩니다. 이것은 각각 맨 위의 식의 우측과 좌측입니다. 각각 0 과 1 값을 가지기에 맨 위의 식이 해당이 되는 것입니다.

Xn→DX,  Yn→P0,  Xn+Yn→DXX_n \xrightarrow{D} X,\; Y_n \xrightarrow{P} 0, \;X_n + Y_n \xrightarrow{D} X

위의 것이랑 사실상 똑같으며, 이것은 Slutsky’s Theorem 의 기초가 된다.

Continuous Mapping Theorem

Xn→DXX_n \xrightarrow{D} X, 그리고 gg 는 XX 의 Support 내에서 Continuous 하겠다고 두겠습니다. 그러면 다음과 같이 Continuous Mapping Theorem 을 적용시킬 수 있습니다.

g(Xn)→Dg(X)g(X_n) \xrightarrow{D} g(X)

이것을 활용시켜서 표준정규분포의 제곱이 카이제곱분포를 따라간다는 것을 구체적으로 알 수가 있습니다.

Slutsky’s Theorem

아마 이 절에서 가장 활용하기 어려운 정리인 것 같습니다.

If Xn→DX,An→Pa, and Bn→Pb, then  An+BnXn→Da+bX\text{If } X_n \xrightarrow{D} X, A_n \xrightarrow{P} a \text {, and } B_n \xrightarrow{P} b \text {, then } \\~\\ A_n+B_n X_n \xrightarrow{D} a+b X

증명은 따로 알 필요 없으며, 이것을 통해서 Convergence in Distribution 을 증명할 수 있다는 것만 알면 됩니다. 이것을 수리통계학 1 에서도 잠시 나오지만 CLT 와 활용해서 극한분포를 알아내는 것이 가능합니다.

Bounded in Probability

Sequence {Xn}\{X_n\} is Bounded in Probability, if ∀ϵ>0\forall \epsilon >0, there exists a constant Bϵ>0B_{\epsilon} >0 and an integer NϵN_{\epsilon} such that

n≥Nϵ→P⁡[∣Xn∣≤Bϵ]≥1−ϵn \ge N_\epsilon \rightarrow \operatorname{P}\left[|X_n| \le B_\epsilon \right] \ge 1 - \epsilon

이것은 Convergence 랑은 비슷한 개념이지만 엄밀히 따지자면 다릅니다. 이것의 의미는 nn 값이 크면 클 수록 수열이 발산하지 않는다는 뜻입니다.

Convergence 의 경우 1개의 값으로 수렴한다는 뜻을 가지지만, Bounded in Probability 는 1개의 값으로 수렴한다는 의미를 가지지는 않고 있습니다.

추가적으로 만약에 해당 수열이 Random Variables, 확률변수들이면 다음이 성립이 됩니다.

If Xn→DX, then {Xn} is bounded in probability\text{If } X_n \xrightarrow{D} X, \text{ then } \{X_n\} \text{ is bounded in probability}
XnYn→P0X_n Y_n \xrightarrow{P} 0

{Xn}\{X_n\} 이 Bounded in Probability 하며

{Yn}\{Y_n\} 이 Convergent to 0 이면

XnYn→P0X_n Y_n \xrightarrow{P} 0 가 성립이 됩니다.

증명 :

저희가 증명할 것은 다음과 같습니다.

lim⁡n→∞P⁡(∣Xn⋅Yn−0∣≥ϵ)\lim_{n \rightarrow \infty}\operatorname{P}(|X_n \cdot Y_n - 0| \ge \epsilon)

Bounded in Probability 하다는 것은 다음을 만족하는 Bϵ>0B_\epsilon >0 과 integer NϵN_\epsilon 이 존재하다는 것입니다.

n≥Nϵ→P⁡[∣Xn∣≤Bϵ]≥1−ϵ ↔P⁡[∣Xn∣>Bϵ]≤ϵn \ge N_\epsilon \rightarrow \operatorname{P}\left[|X_n| \le B_\epsilon \right] \ge 1 - \epsilon \\~\\ \leftrightarrow \operatorname{P}\left[|X_n| > B_\epsilon \right] \le \epsilon

Convergence in Probability to 0 하다는 것은 다음이 해당이 됩니다.

lim⁡n→∞P⁡(∣Yn−0∣≥ϵ)=0\lim_{n \rightarrow \infty}\operatorname{P}(|Y_n - 0| \ge \epsilon) = 0

위의 2가지 식을 응용하면 다음과 같은 결론을 지을 수가 있습니다.

lim‾⁡n→∞P⁡(∣Xn⋅Yn−0∣≥ϵ)=lim‾⁡n→∞P⁡(∣Xn⋅Yn∣≥ϵ)=lim‾⁡n→∞P⁡({∣Xn⋅Yn∣≥ϵ}∩{∣Xn∣≤Bϵ})+lim‾⁡n→∞P⁡({∣Xn⋅Yn∣≥ϵ}∩{∣Xn∣>Bϵ})≤lim‾⁡n→∞P⁡(∣Yn∣≥ϵ/Bϵ)+ϵ=ϵ\varlimsup_{n \rightarrow \infty}\operatorname{P}(|X_n \cdot Y_n - 0| \ge \epsilon) = \varlimsup_{n \rightarrow \infty}\operatorname{P}(|X_n \cdot Y_n| \ge \epsilon) \\ = \varlimsup_{n \rightarrow \infty}\operatorname{P}(\{|X_n \cdot Y_n| \ge \epsilon\} \cap \{|X_n| \le B_\epsilon\}) + \\ \varlimsup_{n \rightarrow \infty}\operatorname{P}(\{|X_n \cdot Y_n| \ge \epsilon\} \cap \{|X_n| > B_\epsilon\}) \\ \le \varlimsup_{n \rightarrow \infty}\operatorname{P}(|Y_n| \ge \epsilon / B_\epsilon) + \epsilon = \epsilon

좌측 Limit :

{∣Xn⋅Yn∣≥ϵ}∩{∣Xn∣≤Bϵ}\{|X_n \cdot Y_n| \ge \epsilon\} \cap \{|X_n| \le B_\epsilon\} : 우측이 만족되기 위해서는 좌측의 식에 YnY_n 과 ϵ\epsilon 을 추가해서 생각해야 됩니다. 즉 ∣Yn∣≥ϵ/∣XN∣≥ϵ/Bϵ|Y_n| \ge \epsilon / |X_N| \ge \epsilon / B_\epsilon 의 조건이 만족됩니다. 그러면 다음이 만족이 됩니다.

lim‾⁡n→∞P⁡({∣Xn⋅Yn∣≥ϵ}∩{∣Xn∣≤Bϵ})≤lim‾⁡n→∞P⁡(∣Yn∣≥ϵ/Bϵ)\varlimsup_{n \rightarrow \infty}\operatorname{P}(\{|X_n \cdot Y_n| \ge \epsilon\} \cap \{|X_n| \le B_\epsilon\}) \le \varlimsup_{n \rightarrow \infty}\operatorname{P}(|Y_n| \ge \epsilon / B_\epsilon)

우측 Limit :

{∣Xn⋅Yn∣≥ϵ}∩{∣Xn∣>Bϵ}\{|X_n \cdot Y_n| \ge \epsilon\} \cap \{|X_n| > B_\epsilon\} : 좌측보다 우측이 더 작기 때문에 우측의 확률을 구하면 되며, 이것은 Bounded in Probability 정의를 갖고오면 다음과 같이 정리가 가능합니다.

lim‾⁡n→∞P⁡({∣Xn⋅Yn∣≥ϵ}∩{∣Xn∣>Bϵ})=lim‾⁡n→∞P⁡(∣Xn∣>Bϵ)=ϵ\varlimsup_{n \rightarrow \infty}\operatorname{P}(\{|X_n \cdot Y_n| \ge \epsilon\} \cap \{|X_n| > B_\epsilon\}) = \varlimsup_{n \rightarrow \infty}\operatorname{P}(|X_n| > B_\epsilon) = \epsilon
Delta Method
g(y)=g(x)+g′(x)(y−x)+o(∣y−x∣)a=o(b), iff ab→0, as b→0Yn=Op(Xn) iff YnXn is bounded in probability as n→∞g(y) = g(x) + g^{\prime}(x) (y-x) + o(|y-x|) \\ a = o(b), \text{ iff } \frac{a}{b} \rightarrow 0, \text{ as } b \rightarrow 0 \\ Y_n = O_p(X_n) \text{ iff } \frac{Y_n}{X_n} \text{ is bounded in probability as } n \rightarrow \infty

첫번째 식이 의미하는 것은 다른게 아닌 Taylor-Series Approximation 입니다.

두번째랑 세번째 식은 o-notation 들이 뭔지에 대한 설명이며, 간단히 말하자면 0 으로 수렴하는 에러 값들이라고 생각해도 됩니다.

이것을 이용해서 다음과 같은 정리를 사용할 수 있습니다.

{Yn}:Bounded in ProbabilityXn=op(Yn)Then, Xn→P0, as n→∞\{Y_n\} : \text{Bounded in Probability} \\ X_n = o_p(Y_n) \\ \text{Then, } X_n \xrightarrow{P} 0, \text{ as } n \rightarrow \infty

증명 :

P⁡[∣Xn∣≥ϵ]=P⁡[{∣Xn∣≥ϵ}∩{∣Yn∣≤Bϵ′}]+P⁡[{∣Xn∣≥ϵ}∩{∣Yn∣>Bϵ′}]≤P⁡[Xn∣Yn∣≥ϵBϵ′]+P⁡[∣Yn∣>Bϵ′]<ϵ′2+[1−P⁡(∣Yn∣≤Bϵ′)]<ϵ′2+[1−(1−ϵ′2)]=ϵ′2+ϵ′2=ϵ′\operatorname{P}[|X_n| \ge \epsilon] = \operatorname{P}[\{|X_n| \ge \epsilon\} \cap \{|Y_n| \le B_{\epsilon^\prime}\}] \\ + \operatorname{P}[\{|X_n| \ge \epsilon\} \cap \{|Y_n| > B_{\epsilon{^\prime}}\}]\\ \le \operatorname{P}\left[\frac{X_n}{|Y_n|} \ge \frac{\epsilon}{B_{\epsilon^{\prime}}} \right] + \operatorname{P}\left[|Y_n| > B_{\epsilon^{\prime}} \right] \\ < \frac{\epsilon^\prime}{2} + \left[1 - \operatorname{P}(|Y_n| \le B_{\epsilon^\prime}) \right] \\ < \frac{\epsilon^\prime}{2} + \left[1 - (1 - \frac{\epsilon^\prime}{2}) \right] = \frac{\epsilon^\prime}{2} + \frac{\epsilon^\prime}{2} = \epsilon^\prime
Delta Method for Normal Distributions

Let {Xn}\{X_n\} be a sequence of Random Variables such that

n(Xn−θ)→DN(0,σ2)\sqrt{n}\left(X_n-\theta\right) \xrightarrow{D} N\left(0, \sigma^2\right)

If the function gg is differentiable at θ\theta, and g′(θ)≠0g^\prime(\theta) \ne 0, then the following is True

n(g(Xn)−g(θ))→DN(0,σ2(g′(θ))2)\sqrt{n}\left(g\left(X_n\right)-g(\theta)\right) \xrightarrow{D} N\left(0, \sigma^2\left(g^{\prime}(\theta)\right)^2\right)

이게 저희가 이 장에서 Delta Method 를 배우는 가장 큰 이유라고 저는 개인적으로 생각하며, 이것을 활용, 및 증명하는 것들 역시 매우 중요하다고 생각합니다.

증명 :

g(Xn)=g(θ)+g′(θ)(Xn−θ)+op(∣Xn−θ∣)n(g(Xn)−g(θ))=g′(θ)n⋅(Xn−θ)+op(n⋅∣Xn−θ∣)g(X_n) = g(\theta) + g^\prime(\theta)(X_n - \theta) + o_p(|X_n - \theta|) \\ \sqrt{n}(g(X_n) - g(\theta)) = g^\prime(\theta) \sqrt{n} \cdot (X_n - \theta) + o_p(\sqrt{n} \cdot |X_n - \theta|)

식을 정리하자면 첫번째 식은 다름이 아니라 Delta Method 를 적용시키면 첫번째 식과 같이 표현을 할 수 있다는 것입니다.

g(θ)g(\theta) 를 좌측으로, 그리고 양쪽에 n\sqrt{n} 을 곱해주면 2번째 식과 같이 나오게 됩니다.

저희는 이미 n(Xn−θ)→DN(0,σ2)\sqrt{n}\left(X_n-\theta\right) \xrightarrow{D} N\left(0, \sigma^2\right) 이라는 것을 알고 있으며, 이것은 다른 의미로

n(Xn−θ)\sqrt{n}\left(X_n-\theta\right) 가 Bounded In Probability 이라는 것을 알고 있습니다. 그렇기 때문에

op(n⋅∣Xn−θ∣)→P0o_p(\sqrt{n} \cdot |X_n - \theta|) \xrightarrow{P} 0 이라는 것과

g′(θ)n⋅(Xn−θ)→Dg′(θ)⋅N⁡(0,σ2)↔N⁡(0,σ2⋅(g′(θ))2)g^\prime(\theta) \sqrt{n} \cdot (X_n - \theta) \xrightarrow{D} g^\prime(\theta) \cdot \operatorname{N}(0, \sigma^2) \leftrightarrow \operatorname{N}\left(0, \sigma^2 \cdot (g^\prime(\theta))^2\right)

이 사실이라고 알 수가 있습니다.

그렇기 때문에 위에서 2번째 식의 증명이 완료되었습니다.

MGF-Technique for Convergence in Distribution

Convergence in Distribution 을 증명하기 위해서는 MGF 들을 비교해서 증명하는 방법도 존재합니다. 우선은 MGF 가 뭔지 먼저 복습하자면, Moment-Generating Function 이며 다음과 같은 식을 가집니다.

MX(t)=E⁡[etX]M_X(t) = \operatorname{E}\left[e^{tX} \right]

{Xn}\{X_n\} : A sequence of R.V.’s with mgf MXn(t)M_{X_n}(t), that exists for −h<t<h-h < t < h

XX : A random variable with mgf MX(t)M_X(t), that exists for ∣t∣≤h1≤h|t| \le h_1 \le h

If lim⁡n→∞MXn(t)=M(t)\lim_{n \rightarrow \infty} M_{X_n}(t) = M(t) exists for ∣t∣≤h1|t| \le h_1,

Then Xn→DXX_n \xrightarrow{D} X

다시 말해 MGF 가 똑같다면 Convergence in Distribution 하다는 것과 같다는 의미입니다.

Central-Limit Theorem (CLT)
What is CLT?

CLT 의 의미는 표본 평균의 Limiting Distribution 이 정규분포를 따른다는 것입니다.

이것은 모든 통계 교수님들이 동의하실 통계학에서 가장 중요한 정리이기 때문에 이 정리가 무슨 의미를 가졌는지를 완전히 이해하시는 것도 향후 통계 과목/강의 들을 들으실 때에 큰 도움이 될 거라고 생각합니다.

이것이 중요한 이유는 표본이 어떠한 분포를 따르던간에 표본 평균이 정규분포를 따르기 때문에 그토록 더 중요한 것입니다.

그럼 구체적으로 식으로 살피겠습니다.

X1,⋯ ,XnX_1, \cdots, X_n, where it follows any distribution with mean μ\mu and variance σ2<∞\sigma^2 < \infty

Then, Yn=n⋅(Xˉn−μ)σY_n = \frac{\sqrt{n} \cdot (\bar{X}_n - \mu)}{\sigma} has the limiting distribution N⁡(0,1)\operatorname{N}(0, 1), or in other words,

Xˉn−E⁡(Xˉn)Var⁡(Xˉn)→DN⁡(0,1)\frac{\bar{X}_n - \operatorname{E}(\bar{X}_n)}{\sqrt{\operatorname{Var}(\bar{X}_n)}} \xrightarrow{D} \operatorname{N}(0,1)

좌측을 다른 식으로 표현하면 다음과 같습니다.

Xˉn−E⁡(Xˉn)Var⁡(Xˉn)=Xˉn−μσ2/n=n(Xˉn−μ)σ=n⋅(Xˉn−μ)n⋅σ=∑i=1nXi−n⋅μn⋅σ→DN⁡(0,1)\frac{\bar{X}_n - \operatorname{E}(\bar{X}_n)}{\sqrt{\operatorname{Var}(\bar{X}_n)}} = \frac{\bar{X}_n - \mu}{\sqrt{\sigma^2 / n}} = \frac{\sqrt{n} (\bar{X}_n - \mu)}{\sigma} = \frac{n \cdot (\bar{X}_n - \mu)}{\sqrt{n} \cdot\sigma} \\ = \frac{\sum_{i=1}^{n} X_i - n \cdot \mu}{\sqrt{n} \cdot \sigma} \xrightarrow{D} \operatorname{N}(0,1)

이것을 또 다르게 표현하면 다음과 같습니다.

Xˉn−E⁡(Xˉn)Var⁡(Xˉn)→DN⁡(0,1) n⋅(Xˉ−μ)σ→DN⁡(0,1) n(Xˉ−μ)→DN⁡(0,σ2) ∑i=1nXi−n⋅μσ/n→DN⁡(0,1)\frac{\bar{X}_n - \operatorname{E}(\bar{X}_n)}{\sqrt{\operatorname{Var}(\bar{X}_n)}} \xrightarrow{D} \operatorname{N}(0,1) \\~\\ \frac{\sqrt{n}\cdot (\bar{X} - \mu)}{\sigma} \xrightarrow{D} \operatorname{N}(0,1)\\~\\ \sqrt{n}(\bar{X} - \mu) \xrightarrow{D} \operatorname{N}(0, \sigma^2) \\~\\ \frac{\sum_{i=1}^{n} X_i - n\cdot \mu}{\sigma / \sqrt{n}} \xrightarrow{D} \operatorname{N}(0,1)

위의 식들은 사실상 다 똑같은 의미를 가지고 있습니다. 그리고 이것은 모두 다 Central-Limit Theorem 을 의미합니다.