Ch6 추정

Kamator0·2026년 6월 15일

통계적 추정

  • 추정 (estimation)
    모집단으로부터 추출한 표본에서 얻은 통계량을 이용하여 모집단의 특성을 나타내는 값인 모수를 파악하는 것을 추정이라고 한다.

  • 점추정 (point estimation)
    표본이 달라지면 추정 결과도 달라진다. 따라서 추정에 사용하는 통계량은 하나의 확률변수로, 이 통계량을 추정량이라고 한다. 추정량은 가변적이므로 최적의 추적량을 설정하여 가장 보편타당한 추정값을 얻어야한다. 최적의 추정값을 구하는 과정

  • 구간추정 (interval estimation)
    그러나 점추정에 의한 모수의 추정은 표본이 어떻게 선정되는지에 따라서 잘못 추정할 수 있어 오류를 범할 수 있다. 어느 정도 확신을 가지고 모수의 참값이 포함되었을 구간을 추정하는 방법을 구간추정이라고 한다.

추정량과 그 성질

  • 미지의 모수 θ\theta로 추정하는 데 통계량 θ^\hat\theta 을 사용했다면 θ^\hat\thetaθ\theta추정량(estigmator) 라고 하고 θ\theta의 통곗값을 θ\theta추정값(estimate)이라 한다.

비편향추정량

  • 모수 θ\theta 의 점추량 θ^\hat\theta 에 대해여 E(θ^)=θE(\hat\theta) = \theta이면 추정량 θ^\hat\theta비편향성(unbiasedness)을 갖는다하고 θ^\hat\thetaθ\theta비편향추정량(unbiased estimator) 또는 불편추정량이라고 한다. 비편향 추정량이 아닌 추정량은 편향추정량(biased estimator) 이라고 한다.
    θ^\hat\thetaθ\theta의 편향추정량이면 b(θ^)=E(θ^)θb(\hat\theta) = E(\hat\theta) - \theta편향(bias) 라고 한다.

모평균 μ\mu 와 모분산 σ2\sigma^2의 비편향추정량

모평균 μ\mu 와 모분산 σ2\sigma^2의 점추정량인 표본평균과 포본분산은 각각 비편향추정량이다.

  • 표본평균 Xˉ\bar{X} 는 모평균 μ\mu의 비편향추정량
  • 표본분산 S2S^2은 모분산 σ2\sigma^2의 비편향추정량

(1)

모평균이 μ\mu인 모집단에서 취한 크기 nn인 확률표본을 X1,X2,,XnX_1, X_2, \cdots, X_n이라 하면 이 확률변수들은 서로 독립이며 동일한 분포를 이루므로

E(Xi)=μ(i=1,2,,n)E(X_i)=\mu \quad (i=1,2,\cdots,n)

이다. 한편, 표본평균은

X=1ni=1nXi\overline{X}=\frac{1}{n}\sum_{i=1}^{n}X_i

이므로 다음이 성립한다.

E(X)=E(1ni=1nXi) =1ni=1nE(Xi) =1ni=1nμ =1n(nμ) =μ\begin{aligned} E(\overline{X}) &=E\left(\frac{1}{n}\sum_{i=1}^{n}X_i\right) \ &=\frac{1}{n}\sum_{i=1}^{n}E(X_i) \ &=\frac{1}{n}\sum_{i=1}^{n}\mu \ &=\frac{1}{n}(n\mu) \ &=\mu \end{aligned}

따라서 표본평균 X\overline{X}는 모평균 μ\mu의 비편향추정량이다.

(2)

표본분산 S2S^2을 다음과 같이 변형할 수 있다.

E(S2)=1n1E{i=1n(Xiμ)2n(Xμ)2}=1n1[i=1nE{(Xiμ)2}nE{(Xμ)2}]=1n1{nσ2nVar(X)}=1n1(nσ2σ2)=σ2\begin{aligned} E(S^2) &=\frac{1}{n-1} E\left\{ \sum_{i=1}^{n}(X_i-\mu)^2 - n(\overline{X}-\mu)^2 \right\} \\ &=\frac{1}{n-1} \left[ \sum_{i=1}^{n}E\left\{(X_i-\mu)^2\right\} - nE\left\{(\overline{X}-\mu)^2\right\} \right] \\ &=\frac{1}{n-1} \left\{ n\sigma^2 - n\operatorname{Var}(\overline{X}) \right\} \\ &=\frac{1}{n-1} \left(n\sigma^2-\sigma^2\right) \\ &=\sigma^2 \end{aligned}

그러므로 표본분산 S2S^2은 모분산 σ2\sigma^2의 비편향추정량이다.

표본표준편차 SS는 모표준편차 σ\sigma의 비편향추정량이 아니다. 그러나 표본의 크기가 크면 (n>=30)(n>=30) 표본표준편차의 편향은 무시할 수 있을 정도로 작아지므로 보통 σ^=S\hat\sigma = S를 사용한다.

유효추정량

모수 θ\theta의 두 비편향추정량을 θ^1\hat{\theta}_1, θ^2\hat{\theta}_2이라 하고, 각각의 분산을 Var(θ^1)\text{Var}(\hat{\theta}_1), Var(θ^2)\text{Var}(\hat{\theta}_2)이라 할 때, Var(θ^1)<Var(θ^2)\text{Var}(\hat{\theta}_1) < \text{Var}(\hat{\theta}_2)이면 θ^1\hat{\theta}_1θ^2\hat{\theta}_2보다 더 유효한 추정량(more efficient estimator)이라 한다.

그리고 분산비 Var(θ^1)Var(θ^2)\dfrac{\text{Var}(\hat{\theta}_1)}{\text{Var}(\hat{\theta}_2)}θ^1\hat{\theta}_1에 대한 θ^2\hat{\theta}_2효율(efficiency)이라 한다. 모수 θ\theta의 가능한 모든 비편향추정량 중 분산이 가장 작은 추정량은 최소분산비편향추정량(minimum variance unbiased estimator) 또는 최량유효추정량(the most efficient estimator)이라고 한다.

Var(θ^1)<Var(θ^2)<Var(θ^3)\text{Var}(\hat{\theta}_1) < \text{Var}(\hat{\theta}_2) < \text{Var}(\hat{\theta}_3)인 모수 θ\theta의 비편향추정량 θ^1\hat{\theta}_1, θ^2\hat{\theta}_2, θ^3\hat{\theta}_3의 확률분포를 나타낸다. 이때 θ^1\hat{\theta}_1이 최소분산비편향추정량이다.

일치추정량(consistent estimator)

모수 θ\theta의 추정량을 θ^\hat{\theta}이라 할 때, 임의의 양수 ϵ\epsilon에 대하여 θ^\hat{\theta}θ\theta에 확률 수렴하면, 즉

limnP(θ^θ<ϵ)=1또는limnP(θ^θϵ)=0\lim_{n \to \infty} \text{P}(|\hat{\theta} - \theta| < \epsilon) = 1 \quad \text{또는} \quad \lim_{n \to \infty} \text{P}(|\hat{\theta} - \theta| \ge \epsilon) = 0

이 성립하면 θ^\hat{\theta}θ\theta일치추정량이라고 한다. 추정량의 일치성을 밝히는 데에는 체비셰프 부등식(Chebyshev inequality)이 유용하다. 체비셰프 부등식은 확률변수 XX의 평균과 분산이 각각 μ\muσ2\sigma^2일 때, 임의의 양수 ϵ\epsilon에 대하여 다음이 성립함을 의미한다.

P(Xμ<ϵ)1σ2ϵ2또는P(Xμϵ)σ2ϵ2\text{P}(|X - \mu| < \epsilon) \ge 1 - \frac{\sigma^2}{\epsilon^2} \quad \text{또는} \quad \text{P}(|X - \mu| \ge \epsilon) \le \frac{\sigma^2}{\epsilon^2}

모평균 μ\mu, 모분산 σ2\sigma^2, 모비율 pp의 일치추정량

모평균 μ\mu, 모분산 σ2\sigma^2, 모비율 pp의 점추정량인 표본평균, 표본분산, 표본비율은 각각 일치추정량이다.

(1) 표본평균 X\overline{X}는 모평균 μ\mu의 일치추정량이다.

(2) 표본분산 S2S^2은 모분산 σ2\sigma^2의 일치추정량이다.

(3) 표본비율 p^\hat{p}은 모비율 pp의 일치추정량이다.

증명

(1) E(X)=μ\text{E}(\overline{X}) = \mu, Var(X)=σ2n\text{Var}(\overline{X}) = \dfrac{\sigma^2}{n}이므로 다음이 성립한다.

P(Xμϵ)Var(X)ϵ2=σ2nϵ2\text{P}(|\overline{X} - \mu| \ge \epsilon) \le \frac{\text{Var}(\overline{X})}{\epsilon^2} = \frac{\sigma^2}{n\epsilon^2}

limnP(Xμϵ)=limnσ2nϵ2=0\lim_{n \to \infty} \text{P}(|\overline{X} - \mu| \ge \epsilon) = \lim_{n \to \infty} \frac{\sigma^2}{n\epsilon^2} = 0

따라서 표본평균 X\overline{X}는 모평균 μ\mu의 일치추정량이다.

충분통계량(sufficient statistic)

이제 추정량의 또 다른 성질인 충분성을 알아보자. 충분성이란 추정량이 표본이 제공하는 모수의 모든 정보를 이용했다는 의미이다. 즉 해당 추정량이 표본의 모든 관측값을 이용할 뿐만 아니라, 관측값이 제공하는 정보도 이용했음을 의미한다.

모수 θ\theta를 추정하는 데 있어서 통계량 θ^\hat{\theta}θ\theta의 추정량으로 사용할 때 통계량 θ^\hat{\theta}이 표본에 포함된 θ\theta의 모든 정보를 갖고 있으면, 통계량 θ^\hat{\theta}θ\theta충분통계량이라고 한다.

모수 θ\theta를 갖는 분포로부터 확률표본 X1,X2,,XnX_1, X_2, \cdots, X_n을 추출했을 때, 그 추정량을 θ^\hat{\theta}이라 하자. 만일 θ^\hat{\theta}이 주어졌을 때, X1,X2,,XnX_1, X_2, \cdots, X_n의 조건부분포

f(x1,x2,,xnθ^)=f(x1,x2,,xn,θ^)g(θ^)f(x_1, x_2, \cdots, x_n \mid \hat{\theta}) = \frac{f(x_1, x_2, \cdots, x_n, \hat{\theta})}{g(\hat{\theta})}

θ\theta가 독립이면 θ^\hat{\theta}θ\theta의 충분통계량이라고 한다.

점추정

점추정량을 구하는 방법에는 적률법, 최대우도법, 최소제곱법 등이 있다. 여기서는 적률법과 최대우도법을 살펴보자.

적률법

모집단으로부터 추출한 확률표본의 표본적률을 모집단에 대응하는 적률과 같게 둠으로써 방정식을 얻을 수 있다. 이 방정식을 미지의 모수에 관해 풀어서 추정값을 구하는 방법을 적률법(method of moments)이라고 한다.

최대우도법

미지의 모수 θ\theta를 갖는 확률변수 XX의 모집단 확률분포를 f(x;θ)f(x;\theta)라 하자. 이 모집단으로부터 추출한 크기 nn인 확률표본이 X1,X2,,XnX_1, X_2, \cdots, X_n일 때, X1,X2,,XnX_1, X_2, \cdots, X_n의 결합확률밀도함수는 다음과 같다.

f(x1,x2,,xn;θ)=f(x1;θ)f(x2;θ)f(xn;θ)f(x_1, x_2, \cdots, x_n;\theta) = f(x_1;\theta)\,f(x_2;\theta)\cdots f(x_n;\theta)

결합확률밀도함수는 θ\theta의 함수로 볼 수 있으므로 다음과 같이 표기한다.

L(θ)=f(x1;θ)f(x2;θ)f(xn;θ)L(\theta) = f(x_1;\theta)\,f(x_2;\theta)\cdots f(x_n;\theta)

L(θ)L(\theta)를 확률표본의 우도함수(likelihood function) 또는 가능도함수라고 한다. 즉 우도함수와 결합확률밀도함수는 같으나, 함수를 모수 θ\theta의 함수로 볼 때 이를 우도함수라고 한다. 우도함수 L(θ)L(\theta)θ\theta 대신 통계량 θ^\hat{\theta}을 대입했을 때 L(θ)L(\theta)가 최대이면, θ^\hat{\theta}을 모수 θ\theta최대우도추정량(maximum likelihood estimator) 또는 최대가능도추정량이라고 한다. 최대우도 추정량 θ^\hat{\theta}은 다음과 같은 방정식으로 구할 수 있다.

L(θ)θ=0(또는lnL(θ)θ=0)\frac{\partial L(\theta)}{\partial \theta} = 0 \quad \left(\text{또는} \quad \frac{\partial \ln L(\theta)}{\partial \theta} = 0\right)

이 방정식을 최대우도방정식(maximum likelihood equation)이라고 한다. 그런데 L(θ)L(\theta)를 최대로 만드는 θ\theta를 구하기보다 lnL(θ)\ln L(\theta)를 최대로 만드는 θ\theta를 구하는 것이 쉬울 때가 많다. 이처럼 우도함수를 이용하여 주어진 표본의 값이 나올 확률을 최대로 만드는 추정량을 구하는 방법을 최대우도법(maximum likelihood method)이라고 한다.

일반적으로 모집단의 모수인 모평균, 모분산, 모표준편차, 모비율을 추정하는 추정량으로 보통 표본평균, 표본분산, 표본표준편차, 표본비율을 각각 사용하는데 이들은 앞에서 살펴본 추정량의 성질을 대부분 만족한다.

구간추정

구간추정은 미지의 모수 θ\theta의 참값이 속할 것이라 믿는 구간을 추정하는 것이다. 미지의 모수 θ\theta의 구간추정을 위해 표본에서 추정량 θ^1\hat{\theta}_1, θ^2\hat{\theta}_2을 정한다. 주어진 양수 α\alpha (0<α<10 < \alpha < 1)에 대해

P(θ^1<θ<θ^2)=1α\text{P}(\hat{\theta}_1 < \theta < \hat{\theta}_2) = 1 - \alpha

를 만족할 때, 구간 (θ^1,θ^2)(\hat{\theta}_1, \hat{\theta}_2)θ\theta100(1α)100(1-\alpha)% 신뢰구간(confidence interval) 또는 구간추정량(interval estimator)이라 한다. 이때 1α1 - \alpha신뢰계수(confidence coefficient) 또는 신뢰도(degree of confidence)라 하고, θ^1\hat{\theta}_1θ^2\hat{\theta}_2신뢰한계(confidence limit)라 하며, θ^1\hat{\theta}_1신뢰하한(confidence lower limit), θ^2\hat{\theta}_2신뢰상한(confidence upper limit)이라 한다.

이때 α\alpha는 오차확률이고, 표본 결과를 이용하여 모집단의 모수를 추측하기 때문에 필연적으로 오차가 발생한다. 오차를 무릅쓰고라도 의미 있는 신뢰구간을 추정해야 하므로 보통 α\alpha는 0.01, 0.05, 0.1을 많이 사용한다.

모평균의 추정

모평균의 점추정

모평균 μ\mu의 추정량으로 여러 가지 통계량을 사용할 수 있는데 가장 대표적인 것이 표본평균이다. 모집단으로부터 추출한 크기 nn인 확률표본을 X1,X2,,XnX_1, X_2, \cdots, X_n이라 할 때, 모평균의 추정량인 표본평균은 다음과 같이 정의한다.

X=1ni=1nXi\overline{X} = \frac{1}{n}\sum_{i=1}^{n} X_i

모평균이 μ\mu이고 모분산이 σ2\sigma^2인 무한 모집단으로부터 크기 nn인 표본을 복원추출할 때, 표본평균 X\overline{X}는 다음과 같은 성질을 갖는다.

E(X)=μ,Var(X)=σ2n\text{E}(\overline{X}) = \mu, \quad \text{Var}(\overline{X}) = \frac{\sigma^2}{n}

이때 추정량의 표준편차 σn\dfrac{\sigma}{\sqrt{n}}는 어떤 불편추정량이 얼마나 좋은지를 나타내는 추정량의 정도를 의미하는 측도로서 표준오차(standard error)라 하고, 기호 SE(X)\text{SE}(\overline{X})로 표기한다. 그러나 일반적으로 σ\sigma는 미지의 수이므로 σ\sigma 대신 표본표준편차 SS를 대입하여 Sn\dfrac{S}{\sqrt{n}}를 사용하는데, 이를 표준오차의 추정값이라 한다.

한편, 정규 모집단으로부터 추출한 크기가 nn인 확률표본의 표본평균 X\overline{X}는 정규분포 N ⁣(μ,σ2n)\text{N}\!\left(\mu, \dfrac{\sigma^2}{n}\right)을 따르므로 이를 표준화한 다음 ZZ-통계량은 표준정규분포 N(0,1)\text{N}(0, 1)을 따른다.

Z=Xμσ/nZ = \frac{\overline{X} - \mu}{\sigma / \sqrt{n}}

따라서 표준정규분포의 성질로부터 임의의 양수 α\alpha (0<α<10 < \alpha < 1)에 대해 다음이 성립한다.

P(zα/2<Z<zα/2)=1α\text{P}(-z_{\alpha/2} < Z < z_{\alpha/2}) = 1 - \alpha

이 결과를 그래프로 나타내면 그림과 같고, 구간추정에서 zα/2z_{\alpha/2} 값으로 z0.05=1.645z_{0.05} = 1.645, z0.025=1.96z_{0.025} = 1.96, z0.005=2.58z_{0.005} = 2.58을 주로 사용한다.

이때 절댓값을 이용하면 P(Z<zα/2)=1α\text{P}(|Z| < z_{\alpha/2}) = 1 - \alpha이고, ZZ-통계량을 대입하면 다음과 같이 정리할 수 있다.

P ⁣(Xμσ/n<zα/2)=1α또는P ⁣(Xμ<zα/2σn)=1α\text{P}\!\left(\left|\frac{\overline{X} - \mu}{\sigma / \sqrt{n}}\right| < z_{\alpha/2}\right) = 1 - \alpha \quad \text{또는} \quad \text{P}\!\left(|\overline{X} - \mu| < z_{\alpha/2}\frac{\sigma}{\sqrt{n}}\right) = 1 - \alpha

여기서 Xμ|\overline{X} - \mu|는 |(추정값) - (참값)|으로 오차를 나타낸다. 이때 zα/2σnz_{\alpha/2}\dfrac{\sigma}{\sqrt{n}}오차한계(limit of error)이므로, 모평균 μ\mu에 대한 추정값의 100(1α)100(1-\alpha)% 오차한계는 zα/2σnz_{\alpha/2}\dfrac{\sigma}{\sqrt{n}}이다.

한편, 오차한계의 크기는 표본의 크기에 따라 달라지므로 오차한계는 조사 결과의 신뢰성과 정확성을 판단하는 근거가 된다.

모평균의 점추정 정리

모평균의 점추정은 다음과 같이 정리할 수 있다.

(1) 모평균 μ\mu의 점추정량은 표본평균 X\overline{X}이다.

(2) 추정량의 표준오차는 추정량 X\overline{X}의 표준편차 σn\dfrac{\sigma}{\sqrt{n}}이다. 즉 다음이 성립한다.

SE(X)=σn\text{SE}(\overline{X}) = \frac{\sigma}{\sqrt{n}}

(3) 모평균 μ\mu에 대한 추정값의 100(1α)100(1-\alpha)% 오차한계는 zα/2σnz_{\alpha/2}\dfrac{\sigma}{\sqrt{n}}이다.

모평균의 구간추정

모분산 σ2\sigma^2을 아는 경우와 모르는 경우로 나누어 모평균의 구간추정을 살펴보자.

모분산 σ2\sigma^2을 아는 경우

정규 모집단 N(μ,σ2)\text{N}(\mu, \sigma^2)에서 모분산 σ2\sigma^2을 알 때, 미지의 모평균 μ\mu의 신뢰구간을 구해보자.

정규 모집단으로부터 추출한 크기 nn인 확률표본의 표본평균 X\overline{X}는 정규분포 N ⁣(μ,σ2n)\text{N}\!\left(\mu, \dfrac{\sigma^2}{n}\right)을 따른다. 따라서 이를 표준화한 다음 ZZ-통계량은 표준정규분포 N(0,1)\text{N}(0, 1)을 따른다.

Z=Xμσ/nZ = \frac{\overline{X} - \mu}{\sigma / \sqrt{n}}

그리고 앞에서 살펴보았듯이 표준정규분포의 성질로부터 임의의 양수 α\alpha (0<α<10 < \alpha < 1)에 대하여 다음이 성립한다.

P ⁣(Xμ<zα/2σn)=1α\text{P}\!\left(|\overline{X} - \mu| < z_{\alpha/2}\frac{\sigma}{\sqrt{n}}\right) = 1 - \alpha

그러므로 모평균 μ\mu에 대한 추정값의 100(1α)100(1-\alpha)% 오차한계는 zα/2σnz_{\alpha/2}\dfrac{\sigma}{\sqrt{n}}이고, 위 식의 좌변 괄호 속 부등식을 μ\mu에 대해 풀면 다음과 같다.

P ⁣(Xzα/2σn<μ<X+zα/2σn)=1α\text{P}\!\left(\overline{X} - z_{\alpha/2}\frac{\sigma}{\sqrt{n}} < \mu < \overline{X} + z_{\alpha/2}\frac{\sigma}{\sqrt{n}}\right) = 1 - \alpha

모분산을 아는 경우 모평균의 구간추정 정리

모분산 σ2\sigma^2을 알 때, 모평균 μ\mu100(1α)100(1-\alpha)% 신뢰구간은 다음과 같다.

(Xzα/2σn,X+zα/2σn)\left(\overline{X} - z_{\alpha/2}\frac{\sigma}{\sqrt{n}}, \quad \overline{X} + z_{\alpha/2}\frac{\sigma}{\sqrt{n}}\right)

즉 모평균 μ\mu100(1α)100(1-\alpha)% 신뢰구간은 다음과 같다.

추정값±100(1α)% 오차한계\text{추정값} \pm 100(1-\alpha)\% \ \text{오차한계}

이는 모평균 μ\mu가 위 구간에 포함됨을 100(1α)100(1-\alpha)% 신뢰할 수 있다는 의미이다. 또한 정규분포를 따르지 않는 모집단에서도 모분산 σ2\sigma^2을 알고 있다면, nn이 대표본(n30n \ge 30)인 경우 중심극한정리에 따라 위 식을 미지의 모평균 μ\mu의 신뢰구간으로 그대로 사용한다.


모분산 σ2\sigma^2을 모르는 경우

정규 모집단 N(μ,σ2)\text{N}(\mu, \sigma^2)에서 모분산 σ2\sigma^2을 모를 때, 미지의 모평균 μ\mu의 신뢰구간을 구해보자.

우선 정규 모집단으로부터 추출한 표본의 크기 nn이 대표본(n30n \ge 30)일 경우 모평균 μ\mu100(1α)100(1-\alpha)% 신뢰구간은 모표준편차 σ\sigma 대신 표본표준편차 SS를 대입하여 다음과 같이 구할 수 있다.

모분산을 모르는 경우 모평균의 구간추정 정리

모분산 σ2\sigma^2을 모를 때, nn이 대표본(n30n \ge 30)일 경우 모평균 μ\mu100(1α)100(1-\alpha)% 신뢰구간은 다음과 같다. 여기서 SS는 표본표준편차이다.

(Xzα/2Sn,X+zα/2Sn)\left(\overline{X} - z_{\alpha/2}\frac{S}{\sqrt{n}}, \quad \overline{X} + z_{\alpha/2}\frac{S}{\sqrt{n}}\right)


정규분포를 따르지 않는 모집단에서 nn이 대표본일 경우, 중심극한정리에 의해 X\overline{X}는 근사적으로 정규분포를 따르므로 모표준편차 σ\sigma는 표본표준편차 SS로 추정할 수 있다. 따라서 모평균 μ\mu의 구간추정은 식을 그대로 사용할 수 있다.

식을 이용하여 신뢰구간을 구하는 경우, 표본에 따라 X\overline{X}SS가 변하므로 신뢰구간도 변한다. 따라서 신뢰구간을 반복하여 구할 때 모평균 μ\mu를 포함하지 않는 신뢰구간도 있지만 전체의 약 100(1α)100(1-\alpha)% 정도는 모평균 μ\mu를 포함할 것으로 기대할 수 있다.

그러나 표본의 크기 nn이 소표본(n<30n < 30)일 경우에는 모집단이 정규분포를 따르더라도 σ\sigma 대신 SS를 대입한 결과가 정규분포와는 크게 다를 수 있다. 이 경우에는 다른 분포가 필요하다. 따라서 ZZ-통계량에서 σ\sigma 대신 SS를 대입한 다음과 같은 TT-통계량을 이용한다.

T=XμS/nT = \frac{\overline{X} - \mu}{S / \sqrt{n}}

TT-통계량은 자유도 (n1)(n-1)tt-분포를 따른다. 그러므로 임의의 양수 α\alpha (0<α<10 < \alpha < 1)에 대하여 다음이 성립한다.

P(tα/2(n1)<T<tα/2(n1))=1α\text{P}(-t_{\alpha/2}(n-1) < T < t_{\alpha/2}(n-1)) = 1 - \alpha

TT-통계량을 대입하면 다음과 같다.

P ⁣(tα/2(n1)<XμS/n<tα/2(n1))=1α\text{P}\!\left(-t_{\alpha/2}(n-1) < \frac{\overline{X} - \mu}{S / \sqrt{n}} < t_{\alpha/2}(n-1)\right) = 1 - \alpha

여기서 tα/2(n1)t_{\alpha/2}(n-1)은 자유도 (n1)(n-1)tt-분포의 오른쪽 꼬리 넓이가 α2\dfrac{\alpha}{2}인 점을 나타낸다. 위 식의 좌변 괄호 속 부등식을 μ\mu에 대해 풀면 다음과 같이 쓸 수 있다.

P ⁣(Xtα/2(n1)Sn<μ<X+tα/2(n1)Sn)=1α\text{P}\!\left(\overline{X} - t_{\alpha/2}(n-1)\frac{S}{\sqrt{n}} < \mu < \overline{X} + t_{\alpha/2}(n-1)\frac{S}{\sqrt{n}}\right) = 1 - \alpha

따라서 다음과 같이 모평균의 구간추정을 구할 수 있다.

모분산을 모르는 경우 모평균의 구간추정 (소표본)**

모분산 σ2\sigma^2을 모를 때, nn이 소표본(n<30n < 30)일 경우 모평균 μ\mu100(1α)100(1-\alpha)% 신뢰구간은 다음과 같다. 여기서 SS는 표본표준편차이다.

(Xtα/2(n1)Sn,X+tα/2(n1)Sn)\left(\overline{X} - t_{\alpha/2}(n-1)\frac{S}{\sqrt{n}}, \quad \overline{X} + t_{\alpha/2}(n-1)\frac{S}{\sqrt{n}}\right)

모평균 차의 추정

각각 정규분포 N(μ1,σ12)N(\mu_1, \sigma_1^2), N(μ2,σ22)N(\mu_2, \sigma_2^2)을 따르면서 서로 독립인 두 모집단에서 두 모평균의 차인 μ1μ2\mu_1-\mu_2의 신뢰구간을 두 모분산 σ12\sigma_1^2σ22\sigma_2^2을 아는 경우와 모르는 경우로 나누어서 살펴보자. 또한 두 모평균을 비교할 때 두 모집단이 서로 독립이 아닌 경우에 사용하는 짝비교(paired comparison) 방법을 이용하여 모평균 차의 신뢰구간을 구해보자.

두 모분산 σ12\sigma_1^2σ22\sigma_2^2을 아는 경우

각각 정규분포 N(μ1,σ12)N(\mu_1, \sigma_1^2), N(μ2,σ22)N(\mu_2, \sigma_2^2)을 따르는 두 모집단에 대하여 크기가 각각 n1n_1, n2n_2이고 서로 독립인 두 확률표본의 표본평균을 각각 X1\overline{X}_1, X2\overline{X}_2라 하면, 통계량 X1X2\overline{X}_1-\overline{X}_2의 기대값과 분산은 각각 다음과 같다.

E(X1X2)=μ1μ2E(\overline{X}_1-\overline{X}_2)=\mu_1-\mu_2
Var(X1X2)=σ12n1+σ22n2\operatorname{Var}(\overline{X}_1-\overline{X}_2)= \frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}

X1X2\overline{X}_1-\overline{X}_2를 표준화한 ZZ-통계량은 다음과 같이 표준정규분포 N(0,1)N(0,1)을 따른다.

Z=(X1X2)(μ1μ2)σ12n1+σ22n2N(0,1)Z= \frac{(\overline{X}_1-\overline{X}_2)-(\mu_1-\mu_2)} {\sqrt{\frac{\sigma_1^2}{n_1}+\frac{\sigma_2^2}{n_2}}} \sim N(0,1)

따라서 표준정규분포의 성질로부터 임의의 양수 α(0<α<1)\alpha(0<\alpha<1)에 대하여 다음이 성립한다.

P(zα/2<Z<zα/2)=1αP(-z_{\alpha/2}<Z<z_{\alpha/2})=1-\alpha

위 식에 ZZ-통계량을 대입하여 좌변 괄호 속 부등식을 μ1μ2\mu_1-\mu_2에 대해 풀면 다음과 같이 쓸 수 있다.

P((X1X2)zα/2σ12n1+σ22n2<μ1μ2<(X1X2)+zα/2σ12n1+σ22n2)=1αP\left( (\overline{X}_1-\overline{X}_2)- z_{\alpha/2} \sqrt{ \frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2} } < \mu_1-\mu_2 < (\overline{X}_1-\overline{X}_2) + z_{\alpha/2} \sqrt{ \frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2} } \right) =1-\alpha

따라서 모평균 차의 구간추정은 다음과 같이 구할 수 있다.

두 모분산을 아는 경우 모평균 차의 구간추정

모분산 σ12\sigma_1^2σ22\sigma_2^2을 알 때, 모평균 차 μ1μ2\mu_1-\mu_2100(1α)%100(1-\alpha)\% 신뢰구간은 다음과 같다.

((X1X2)zα/2σ12n1+σ22n2,(X1X2)+zα/2σ12n1+σ22n2)\left((\overline{X}_1-\overline{X}_2)- z_{\alpha/2} \sqrt{ \frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2} }, \quad (\overline{X}_1-\overline{X}_2) + z_{\alpha/2} \sqrt{ \frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2} } \right)

또한 정규분포를 따르지 않는 두 모집단에서 모분산 σ12\sigma_1^2σ22\sigma_2^2을 알 때, 모평균 차 μ1μ2\mu_1-\mu_2의 신뢰구간은 n1n_1, n2n_2가 대표본(n130,n230)(n_1 \ge 30, n_2 \ge 30)일 경우 중심극한정리에 의해 [정리 6-7]의 식을 그대로 사용한다.


두 모분산 σ12\sigma_1^2σ22\sigma_2^2을 모르는 경우

각각 정규분포 N(μ1,σ12)N(\mu_1, \sigma_1^2), N(μ2,σ22)N(\mu_2, \sigma_2^2)을 따르는 두 모집단에 대하여 두 모분산 σ12\sigma_1^2σ22\sigma_2^2을 모를 때, 모평균 차 μ1μ2\mu_1-\mu_2의 신뢰구간을 구해보자. 우선 두 정규 모집단으로부터 추출한 표본의 크기 n1n_1, n2n_2가 대표본(n130,n230)(n_1 \ge 30, n_2 \ge 30)일 경우 μ1μ2\mu_1-\mu_2100(1α)%100(1-\alpha)\% 신뢰구간은 모표준편차 σ1\sigma_1σ2\sigma_2 대신 표본표준편차 S1S_1S2S_2를 대입하여 다음과 같이 구할 수 있다.

두 모분산을 모르는 경우 모평균 차의 구간추정

모분산 σ12\sigma_1^2σ22\sigma_2^2을 모를 때, n1n_1, n2n_2가 대표본(n130,n230)(n_1 \ge 30, n_2 \ge 30)일 경우 모평균 차 μ1μ2\mu_1-\mu_2100(1α)%100(1-\alpha)\% 신뢰구간은 다음과 같다.

((X1X2)zα/2S12n1+S22n2,(X1X2)+zα/2S12n1+S22n2)\left( (\overline{X}_1-\overline{X}_2)- z_{\alpha/2} \sqrt{ \frac{S_1^2}{n_1} + \frac{S_2^2}{n_2} }, \quad (\overline{X}_1-\overline{X}_2) + z_{\alpha/2} \sqrt{ \frac{S_1^2}{n_1} + \frac{S_2^2}{n_2} } \right)

또한 정규분포를 따르지 않는 두 모집단에서 n1n_1, n2n_2가 대표본의 경우, 중심극한정리에 의해 X1X2\overline{X}_1-\overline{X}_2는 근사적으로 정규분포를 따르므로 모표준편차 σ1\sigma_1σ2\sigma_2는 표본표준편차 S1S_1S2S_2로 추정할 수 있다. 따라서 μ1μ2\mu_1-\mu_2의 구간추정은 [정리 6-8]의 식을 그대로 사용할 수 있다.

한편 표본의 크기 n1n_1, n2n_2가 소표본(n1<30,n2<30)(n_1<30, n_2<30)일 경우, 두 모평균 차의 추정에서는 X1X2\overline{X}_1-\overline{X}_2를 표준화한 ZZ-통계량을 사용할 수 없고 tt-분포를 이용해야 한다. 따라서 σ12=σ22\sigma_1^2=\sigma_2^2인 경우와 σ12σ22\sigma_1^2 \ne \sigma_2^2인 경우로 나누어서 신뢰구간을 구해보자.


두 모분산 σ12\sigma_1^2σ22\sigma_2^2을 모르지만 σ12=σ22\sigma_1^2=\sigma_2^2인 경우

두 모집단의 모분산을 모르지만 같다고(σ12=σ22)(\sigma_1^2=\sigma_2^2) 가정할 수 있을 때, 공통분산의 추정은 다음과 같은 합동표본분산 Sp2S_p^2을 이용한다.

Sp2=(n11)S12+(n21)S22n1+n22S_p^2= \frac{ (n_1-1)S_1^2+(n_2-1)S_2^2 } {n_1+n_2-2}

합동표본분산 Sp2S_p^2을 이용하여 X1X2\overline{X}_1-\overline{X}_2를 표준화한 TT-통계량은 다음과 같다.

T=(X1X2)(μ1μ2)Sp1n1+1n2T= \frac{ (\overline{X}_1-\overline{X}_2)-(\mu_1-\mu_2) } { S_p \sqrt{ \frac{1}{n_1} + \frac{1}{n_2} } }

통계량 TT는 자유도 (n1+n22)(n_1+n_2-2)tt-분포를 따른다. 따라서 임의의 양수 α(0<α<1)\alpha(0<\alpha<1)에 대하여 다음이 성립한다.

P(tα/2(n1+n22)<(X1X2)(μ1μ2)Sp1n1+1n2<tα/2(n1+n22))=1αP\left( -t_{\alpha/2}(n_1+n_2-2) < \frac{ (\overline{X}_1-\overline{X}_2)-(\mu_1-\mu_2) } { S_p \sqrt{ \frac{1}{n_1} + \frac{1}{n_2} } } < t_{\alpha/2}(n_1+n_2-2)\right)= 1-\alpha

위 식의 좌변 괄호 속 부등식을 μ1μ2\mu_1-\mu_2에 대해 풀면 다음과 같이 쓸 수 있다.

P((X1X2)tα/2(n1+n22)Sp1n1+1n2<μ1μ2<(X1X2)+tα/2(n1+n22)Sp1n1+1n2)=1αP\left((\overline{X}_1-\overline{X}_2)- t_{\alpha/2}(n_1+n_2-2) S_p \sqrt{ \frac{1}{n_1} + \frac{1}{n_2} } < \mu_1-\mu_2 < (\overline{X}_1-\overline{X}_2) + t_{\alpha/2}(n_1+n_2-2) S_p \sqrt{ \frac{1}{n_1} + \frac{1}{n_2}}\right)= 1-\alpha

두 모분산을 모르지만 σ12=σ22\sigma_1^2=\sigma_2^2인 경우 모평균 차의 구간추정

두 모분산 σ12\sigma_1^2, σ22\sigma_2^2을 모르지만 σ12=σ22\sigma_1^2=\sigma_2^2일 때, n1n_1, n2n_2가 소표본(n1<30,n2<30)(n_1<30, n_2<30)이면 두 모평균 차 μ1μ2\mu_1-\mu_2(1α)100%(1-\alpha)100\% 신뢰구간은 다음과 같다.

((X1X2)tα/2(n1+n22)Sp1n1+1n2,(X1X2)+tα/2(n1+n22)Sp1n1+1n2)\left( (\overline{X}_1-\overline{X}_2)- t_{\alpha/2}(n_1+n_2-2) S_p \sqrt{ \frac{1}{n_1} + \frac{1}{n_2} }, \quad (\overline{X}_1-\overline{X}_2) + t_{\alpha/2}(n_1+n_2-2) S_p \sqrt{ \frac{1}{n_1} + \frac{1}{n_2} } \right)

두 모분산 σ12\sigma_1^2σ22\sigma_2^2을 모르지만 σ12σ22\sigma_1^2 \ne \sigma_2^2인 경우

σ12σ22\sigma_1^2 \ne \sigma_2^2일 때 다음 TT-통계량은 근사적으로 tt-분포를 따른다는 성질을 이용한다.

T=(X1X2)(μ1μ2)S12n1+S22n2T= \frac{ (\overline{X}_1-\overline{X}_2)-(\mu_1-\mu_2) } { \sqrt{ \frac{S_1^2}{n_1} + \frac{S_2^2}{n_2} } }

두 모분산을 모르지만 σ12σ22\sigma_1^2 \ne \sigma_2^2인 경우 모평균 차의 구간추정

두 모분산 σ12\sigma_1^2, σ22\sigma_2^2을 모르지만 σ12σ22\sigma_1^2 \ne \sigma_2^2일 때, n1n_1, n2n_2가 소표본(n1<30,n2<30)(n_1<30, n_2<30)이면 두 모평균 차 μ1μ2\mu_1-\mu_2100(1α)%100(1-\alpha)\% 신뢰구간은 다음과 같다.

((X1X2)tα/2(n)S12n1+S22n2,(X1X2)+tα/2(n)S12n1+S22n2)\left( (\overline{X}_1-\overline{X}_2)- t_{\alpha/2}(n^*) \sqrt{ \frac{S_1^2}{n_1} + \frac{S_2^2}{n_2} }, \quad (\overline{X}_1-\overline{X}_2) + t_{\alpha/2}(n^*) \sqrt{ \frac{S_1^2}{n_1} + \frac{S_2^2}{n_2} } \right)

이때 자유도 nn^*은 근사적으로 다음과 같다. 만약 nn^*가 정수가 아니라면 nn^*를 넘지 않는 최대정수를 자유도로 선택한다.

n=(S12n1+S22n2)2(S12n1)2n11+(S22n2)2n21n^*= \frac{ \left( \frac{S_1^2}{n_1} + \frac{S_2^2}{n_2} \right)^2 } { \frac{ \left( \frac{S_1^2}{n_1} \right)^2 } {n_1-1} + \frac{ \left( \frac{S_2^2}{n_2} \right)^2 } {n_2-1} }

짝비교에 의한 모평균 차의 구간추정

두 모평균의 비교에 있어서 두 모집단이 서로 독립이 아닌 경우에 짝비교 방법을 이용한다. 예를 들어 실험의 효과를 알아보기 위한 실험 전과 후의 비교, 왼손과 오른손의 악력 비교 등 한 개체에서 그 처리 효과를 비교할 때 짝비교 방법을 사용한다. 짝비교의 경우 하나의 짝의 실험 단위에는 처리 1(X1)1(X_1)을, 다른 실험 단위에는 처리 2(X2)2(X_2)를 적용하고, 각 짝에서 얻은 관측값의 차(D=X1X2)(D=X_1-X_2)를 이용하여 두 처리를 비교한다. 짝비교에서의 자료 구조는 다음과 같다.

두 확률변수 X1X_1X2X_2가 정규분포를 따른다고 가정할 때, Di=X1iX2iD_i=X_{1i}-X_{2i}는 정규분포 N(μ1μ2,σD2)N(\mu_1-\mu_2, \sigma_D^2)을 따른다. 이때 σD2\sigma_D^2은 알고 있다고 가정한다. 또한 짝 (X1i,X2i)(X_{1i}, X_{2i})는 서로 독립적으로 추출하며 X1iX_{1i}X2iX_{2i}는 일반적으로 강한 양의 상관관계를 갖는다.

D1,D2,,DnD_1, D_2, \cdots, D_n의 표본평균

D=1ni=1nDi\overline{D}= \frac{1}{n} \sum_{i=1}^{n}D_i

는 정규분포

N(μ1μ2,σD2n)N\left(\mu_1-\mu_2, \frac{\sigma_D^2}{n}\right)

을 따르고, 이를 표준화한 다음 ZZ-통계량은 표준정규분포 N(0,1)N(0,1)을 따른다.

Z=D(μ1μ2)σD/nZ= \frac{ \overline{D}-(\mu_1-\mu_2) } { \sigma_D/\sqrt{n} }

따라서 표준정규분포의 성질로부터 임의의 양수 α(0<α<1)\alpha(0<\alpha<1)에 대하여 다음이 성립한다.

P(zα/2<Z<zα/2)=1αP(-z_{\alpha/2}<Z<z_{\alpha/2})=1-\alpha

위 식에 ZZ-통계량을 대입하여 좌변 괄호 속 부등식을 μ1μ2\mu_1-\mu_2에 대해 풀면 다음과 같이 쓸 수 있다.

P(Dzα/2σDn<μ1μ2<D+zα/2σDn)=1αP\left( \overline{D}- z_{\alpha/2} \frac{\sigma_D}{\sqrt{n}} < \mu_1-\mu_2 < \overline{D} + z_{\alpha/2} \frac{\sigma_D}{\sqrt{n}} \right)= 1-\alpha

짝비교에 의한 모평균 차의 구간추정

두 모평균 차 μ1μ2\mu_1-\mu_2100(1α)%100(1-\alpha)\% 신뢰구간은 다음과 같다.

(Dzα/2σDn,D+zα/2σDn)\left( \overline{D}- z_{\alpha/2} \frac{\sigma_D}{\sqrt{n}}, \quad \overline{D} + z_{\alpha/2} \frac{\sigma_D}{\sqrt{n}} \right)

실제로 σD2\sigma_D^2을 모르는 경우에는 다음 TT-통계량을 이용한다.

T=D(μ1μ2)SD/nT= \frac{ \overline{D}-(\mu_1-\mu_2) } { S_D/\sqrt{n} }

통계량 TT는 자유도 (n1)(n-1)tt-분포를 따른다. 따라서 같은 방법으로 두 모평균 차 μ1μ2\mu_1-\mu_2100(1α)%100(1-\alpha)\% 신뢰구간을 구하면 다음과 같다.

(Dtα/2(n1)SDn,D+tα/2(n1)SDn)\left( \overline{D}- t_{\alpha/2}(n-1) \frac{S_D}{\sqrt{n}}, \quad \overline{D} + t_{\alpha/2}(n-1) \frac{S_D}{\sqrt{n}} \right)

여기서

SD2=1n1i=1n(DiD)2S_D^2= \frac{1}{n-1} \sum_{i=1}^{n} (D_i-\overline{D})^2

이다.

모비율의 추정

어떤 특정한 속성의 모비율이 pp 인 무한 모집단으로부터 n개의 확률표본을 추출했을 때, 특정한 속성을 갖는 개체의 총수를 XX라 하면 확률변수 XX는 이항분포 B(n,p)B(n,p)를 따른다. 이때 모비율의 추정량인 표본비율은 다음과 같이 정의했다.

p^=Xn\hat p = \frac{X}{n}

확률변수 XX가 이항분포 B(n,p)\text{B}(n, p)를 따를 때, 다음이 성립한다.

E(X)=np\text{E}(X) = np

Var(X)=np(1p)\text{Var}(X) = np(1-p)

따라서 표본비율 p^\hat{p}의 평균과 분산은 각각 다음과 같다.

E(p^)=p\text{E}(\hat{p}) = p

Var(p^)=p(1p)n\text{Var}(\hat{p}) = \frac{p(1-p)}{n}

그러므로 p^\hat{p}은 모비율 pp의 불편추정량이고, p^\hat{p}의 표준오차는 p(1p)n\sqrt{\dfrac{p(1-p)}{n}}이며 표준오차의 추정값은 p^(1p^)n\sqrt{\dfrac{\hat{p}(1-\hat{p})}{n}}임을 알 수 있다.

이항 모집단의 모비율 pp의 점추정량이 p^=Xn\hat{p} = \dfrac{X}{n}이고 표본의 크기 nn이 대표본일 때, 확률변수 XX의 분포는 근사적으로 정규분포 N(np,np(1p))\text{N}(np, \, np(1-p))에 가까워진다. 따라서 표본비율 p^\hat{p}의 분포는 근사적으로 정규분포 N ⁣(p,p(1p)n)\text{N}\!\left(p, \, \dfrac{p(1-p)}{n}\right)를 따른다는 사실이 알려져 있다.

이때 대표본의 기준은 흔히 np^>5n\hat{p} > 5, n(1p^)>5n(1-\hat{p}) > 5nn이 충분히 큰 경우를 의미하므로

Z=Xnpnp^(1p^)=p^pp^(1p^)nZ = \frac{X - np}{\sqrt{n\hat{p}(1-\hat{p})}} = \frac{\hat{p} - p}{\sqrt{\dfrac{\hat{p}(1-\hat{p})}{n}}}

는 근사적으로 표준정규분포를 따른다. 따라서 임의의 양수 α\alpha (0<α<10 < \alpha < 1)에 대하여 다음이 성립한다.

P ⁣(zα/2<p^pp^(1p^)n<zα/2)=1α\text{P}\!\left(-z_{\alpha/2} < \frac{\hat{p} - p}{\sqrt{\dfrac{\hat{p}(1-\hat{p})}{n}}} < z_{\alpha/2}\right) = 1 - \alpha

좌변의 괄호 속 부등식을 pp에 대해 풀면 다음과 같다.

P ⁣(p^zα/2p^(1p^)n<p<p^+zα/2p^(1p^)n)=1α\text{P}\!\left(\hat{p} - z_{\alpha/2}\sqrt{\frac{\hat{p}(1-\hat{p})}{n}} < p < \hat{p} + z_{\alpha/2}\sqrt{\frac{\hat{p}(1-\hat{p})}{n}}\right) = 1 - \alpha

따라서 모비율의 구간추정은 다음과 같이 구할 수 있다.

모비율의 구간추정

nn이 대표본(n30n \ge 30)일 경우 모비율 pp100(1α)100(1-\alpha)% 신뢰구간은 다음과 같다.

(p^zα/2p^(1p^)n,p^+zα/2p^(1p^)n)\left(\hat{p} - z_{\alpha/2}\sqrt{\frac{\hat{p}(1-\hat{p})}{n}}, \quad \hat{p} + z_{\alpha/2}\sqrt{\frac{\hat{p}(1-\hat{p})}{n}}\right)


모비율 차의 추정

두 이항 모집단으로부터 추출한 크기가 n1n_1, n2n_2이고 서로 독립인 확률표본의 표본비율을 각각 p^1\hat{p}_1, p^2\hat{p}_2이라 하면 통계량 p^1p^2\hat{p}_1 - \hat{p}_2의 기댓값과 분산은 다음과 같다.

E(p^1p^2)=p1p2\text{E}(\hat{p}_1 - \hat{p}_2) = p_1 - p_2

Var(p^1p^2)=p1(1p1)n1+p2(1p2)n2\text{Var}(\hat{p}_1 - \hat{p}_2) = \frac{p_1(1-p_1)}{n_1} + \frac{p_2(1-p_2)}{n_2}

표본의 크기 n1n_1, n2n_2가 충분히 클 때, p^1p^2\hat{p}_1 - \hat{p}_2을 표준화한 다음 ZZ-통계량은 근사적으로 표준정규분포를 따른다.

Z=(p^1p^2)(p1p2)p^1(1p^1)n1+p^2(1p^2)n2Z = \frac{(\hat{p}_1 - \hat{p}_2) - (p_1 - p_2)}{\sqrt{\dfrac{\hat{p}_1(1-\hat{p}_1)}{n_1} + \dfrac{\hat{p}_2(1-\hat{p}_2)}{n_2}}}

따라서 임의의 양수 α\alpha (0<α<10 < \alpha < 1)에 대하여 다음이 성립한다.

P ⁣(zα/2<(p^1p^2)(p1p2)p^1(1p^1)n1+p^2(1p^2)n2<zα/2)=1α\text{P}\!\left(-z_{\alpha/2} < \frac{(\hat{p}_1 - \hat{p}_2) - (p_1 - p_2)}{\sqrt{\dfrac{\hat{p}_1(1-\hat{p}_1)}{n_1} + \dfrac{\hat{p}_2(1-\hat{p}_2)}{n_2}}} < z_{\alpha/2}\right) = 1 - \alpha

그러므로 좌변의 괄호 속 부등식을 p1p2p_1 - p_2에 대해 풀면 다음과 같다.

P ⁣((p^1p^2)zα/2p^1(1p^1)n1+p^2(1p^2)n2<p1p2<(p^1p^2)+zα/2p^1(1p^1)n1+p^2(1p^2)n2)=1α\text{P}\!\left((\hat{p}_1 - \hat{p}_2) - z_{\alpha/2}\sqrt{\frac{\hat{p}_1(1-\hat{p}_1)}{n_1} + \frac{\hat{p}_2(1-\hat{p}_2)}{n_2}} < p_1 - p_2 < (\hat{p}_1 - \hat{p}_2) + z_{\alpha/2}\sqrt{\frac{\hat{p}_1(1-\hat{p}_1)}{n_1} + \frac{\hat{p}_2(1-\hat{p}_2)}{n_2}}\right) = 1 - \alpha

모비율 차의 구간추정 정리

n1n_1, n2n_2가 대표본(n130n_1 \ge 30, n230n_2 \ge 30)일 경우 두 모비율 차 p1p2p_1 - p_2100(1α)100(1-\alpha)% 신뢰구간은 다음과 같다.

((p^1p^2)zα/2p^1(1p^1)n1+p^2(1p^2)n2,(p^1p^2)+zα/2p^1(1p^1)n1+p^2(1p^2)n2)\left((\hat{p}_1 - \hat{p}_2) - z_{\alpha/2}\sqrt{\frac{\hat{p}_1(1-\hat{p}_1)}{n_1} + \frac{\hat{p}_2(1-\hat{p}_2)}{n_2}}, \quad (\hat{p}_1 - \hat{p}_2) + z_{\alpha/2}\sqrt{\frac{\hat{p}_1(1-\hat{p}_1)}{n_1} + \frac{\hat{p}_2(1-\hat{p}_2)}{n_2}}\right)


모분산의 추정

모평균이 μ\mu이고 모분산이 σ2\sigma^2인 모집단으로부터 추출한 크기 nn인 확률표본을 X1,X2,,XnX_1, X_2, \cdots, X_n이라 할 때, 모분산의 추정량인 표본분산은 다음과 같이 정의한다.

S2=1n1i=1n(XiX)2S^2 = \frac{1}{n-1}\sum_{i=1}^{n}(X_i - \overline{X})^2

무한 모집단으로부터 확률표본을 추출하여 얻은 표본분산 S2S^2은 다음과 같은 성질을 갖는다.

E(S2)=σ2\text{E}(S^2) = \sigma^2

그러므로 S2S^2은 모분산 σ2\sigma^2의 불편추정량임을 알 수 있다. 모표준편차 σ\sigma의 추정량인 표본표준편차(sample standard deviation)는 다음과 같이 정의한다.

S=1n1i=1n(XiX)2S = \sqrt{\frac{1}{n-1}\sum_{i=1}^{n}(X_i - \overline{X})^2}

이때 E(S)<σ\text{E}(S) < \sigma이므로 표본표준편차 SS는 모표준편차 σ\sigma의 불편추정량이 아니다. 그러나 표본의 크기 nn이 크면(n>30n > 30) 편향이 작아서 SSσ\sigma의 추정량으로 사용할 수 있다.

모평균 μ\mu를 아는 경우와 모르는 경우로 나누어 모분산의 구간추정을 살펴보자.


모평균 μ\mu를 아는 경우 모분산 σ2\sigma^2의 구간추정

정규분포 N(μ,σ2)\text{N}(\mu, \sigma^2)을 따르는 모집단에서 모평균 μ\mu를 알 때, 미지의 모분산 σ2\sigma^2의 신뢰구간을 구해보자. 다음과 같이 정의한 χ2\chi^2-통계량은 자유도가 nnχ2\chi^2-분포를 따른다.

χ2=i=1n(Xiμ)2σ2\chi^2 = \frac{\sum_{i=1}^{n}(X_i - \mu)^2}{\sigma^2}

따라서 임의의 양수 α\alpha (0<α<10 < \alpha < 1)에 대하여 다음이 성립한다.

P ⁣(χ1α/22(n)<χ2<χα/22(n))=1α\text{P}\!\left(\chi^2_{1-\alpha/2}(n) < \chi^2 < \chi^2_{\alpha/2}(n)\right) = 1 - \alpha

여기서 χ1α/22(n)\chi^2_{1-\alpha/2}(n)χα/22(n)\chi^2_{\alpha/2}(n)은 각각 100(α/2)100(\alpha/2)% 점과 100(1α/2)100(1-\alpha/2)% 점이다. 앞의 식에 χ2\chi^2-통계량을 대입하면 다음과 같다.

P ⁣(χ1α/22(n)<i=1n(Xiμ)2σ2<χα/22(n))=1α\text{P}\!\left(\chi^2_{1-\alpha/2}(n) < \frac{\sum_{i=1}^{n}(X_i - \mu)^2}{\sigma^2} < \chi^2_{\alpha/2}(n)\right) = 1 - \alpha

좌변의 괄호 속 부등식을 σ2\sigma^2에 대해 풀면 다음과 같이 쓸 수 있다.

P ⁣(i=1n(Xiμ)2χα/22(n)<σ2<i=1n(Xiμ)2χ1α/22(n))=1α\text{P}\!\left(\frac{\sum_{i=1}^{n}(X_i - \mu)^2}{\chi^2_{\alpha/2}(n)} < \sigma^2 < \frac{\sum_{i=1}^{n}(X_i - \mu)^2}{\chi^2_{1-\alpha/2}(n)}\right) = 1 - \alpha

따라서 모평균을 아는 경우 모분산의 신뢰구간은 다음과 같이 구할 수 있다.

모평균을 아는 경우 모분산의 구간추정 정리

모평균 μ\mu를 알 때, 모분산 σ2\sigma^2100(1α)100(1-\alpha)% 신뢰구간은 다음과 같다.

(i=1n(Xiμ)2χα/22(n),i=1n(Xiμ)2χ1α/22(n))\left(\frac{\sum_{i=1}^{n}(X_i - \mu)^2}{\chi^2_{\alpha/2}(n)}, \quad \frac{\sum_{i=1}^{n}(X_i - \mu)^2}{\chi^2_{1-\alpha/2}(n)}\right)


모평균 μ\mu를 모르는 경우 모분산 σ2\sigma^2의 구간추정

정규분포 N(μ,σ2)\text{N}(\mu, \sigma^2)을 따르는 모집단에서 모평균 μ\mu를 모를 때, 미지의 모분산 σ2\sigma^2의 신뢰구간을 구해보자. 다음과 같이 정의한 χ2\chi^2-통계량은 자유도가 (n1)(n-1)χ2\chi^2-분포를 따른다.

χ2=(n1)S2σ2\chi^2 = \frac{(n-1)S^2}{\sigma^2}

따라서 임의의 양수 α\alpha (0<α<10 < \alpha < 1)에 대하여 다음이 성립한다.

P ⁣(χ1α/22(n1)<(n1)S2σ2<χα/22(n1))=1α\text{P}\!\left(\chi^2_{1-\alpha/2}(n-1) < \frac{(n-1)S^2}{\sigma^2} < \chi^2_{\alpha/2}(n-1)\right) = 1 - \alpha

위 식에서 좌변의 괄호 속 부등식을 σ2\sigma^2에 대해 풀면 다음과 같이 쓸 수 있다.

P ⁣((n1)S2χα/22(n1)<σ2<(n1)S2χ1α/22(n1))=1α\text{P}\!\left(\frac{(n-1)S^2}{\chi^2_{\alpha/2}(n-1)} < \sigma^2 < \frac{(n-1)S^2}{\chi^2_{1-\alpha/2}(n-1)}\right) = 1 - \alpha

따라서 모평균을 모르는 경우 모분산의 신뢰구간은 다음과 같이 구할 수 있다.

모평균을 모르는 경우 모분산의 구간추정

모평균 μ\mu를 모를 때, 모분산 σ2\sigma^2100(1α)100(1-\alpha)% 신뢰구간은 다음과 같다.

((n1)S2χα/22(n1),(n1)S2χ1α/22(n1))\left(\frac{(n-1)S^2}{\chi^2_{\alpha/2}(n-1)}, \quad \frac{(n-1)S^2}{\chi^2_{1-\alpha/2}(n-1)}\right)


모분산 비의 구간추정

서로 독립인 두 정규 모집단 N(μ1,σ12)\text{N}(\mu_1, \sigma_1^2), N(μ2,σ22)\text{N}(\mu_2, \sigma_2^2)으로부터 크기 n1n_1, n2n_2인 확률표본을 추출할 때, 표본분산을 S12S_1^2, S22S_2^2이라고 하면 다음과 같은 통계량 FF는 자유도 (n11,n21)(n_1 - 1, \, n_2 - 1)FF-분포를 따른다.

F=S12/σ12S22/σ22F = \frac{S_1^2 / \sigma_1^2}{S_2^2 / \sigma_2^2}

이로부터 다음과 같은 신뢰구간을 얻을 수 있다. 임의의 양수 α\alpha (0<α<10 < \alpha < 1)에 대하여

P ⁣(F1α/2(n11,n21)<S12/σ12S22/σ22<Fα/2(n11,n21))=1α\text{P}\!\left(F_{1-\alpha/2}(n_1-1, \, n_2-1) < \frac{S_1^2/\sigma_1^2}{S_2^2/\sigma_2^2} < F_{\alpha/2}(n_1-1, \, n_2-1)\right) = 1 - \alpha

가 성립하므로 좌변의 괄호 속 부등식을 σ12σ22\dfrac{\sigma_1^2}{\sigma_2^2}에 대해 풀면 다음과 같이 쓸 수 있다.

P ⁣(S12S221Fα/2(n11,n21)<σ12σ22<S12S221F1α/2(n11,n21))=1α\text{P}\!\left(\frac{S_1^2}{S_2^2} \cdot \frac{1}{F_{\alpha/2}(n_1-1, \, n_2-1)} < \frac{\sigma_1^2}{\sigma_2^2} < \frac{S_1^2}{S_2^2} \cdot \frac{1}{F_{1-\alpha/2}(n_1-1, \, n_2-1)}\right) = 1 - \alpha

따라서 모분산 비의 신뢰구간은 다음과 같이 구할 수 있다.

(S12/S22Fα/2(n11,n21),S12/S22F1α/2(n11,n21))\left(\frac{S_1^2 / S_2^2}{F_{\alpha/2}(n_1-1, \, n_2-1)}, \quad \frac{S_1^2 / S_2^2}{F_{1-\alpha/2}(n_1-1, \, n_2-1)}\right)

이 결과는 다음과 같이 정리할 수 있다.

모분산 비의 구간추정

모분산 비 σ12σ22\dfrac{\sigma_1^2}{\sigma_2^2}100(1α)100(1-\alpha)% 신뢰구간은 다음과 같다.

(S12S22F1α/2(n21,n11),S12S22Fα/2(n21,n11))\left(\frac{S_1^2}{S_2^2}\,F_{1-\alpha/2}(n_2-1, \, n_1-1), \quad \frac{S_1^2}{S_2^2}\,F_{\alpha/2}(n_2-1, \, n_1-1)\right)


표본의 크기 결정

모수를 추정하기 위해 표본을 추출할 때 표본의 크기를 얼마로 해야 하는지를 결정해야 한다. 일반적으로 표본의 크기가 너무 크면 시간과 노력이 많이 소비되는 반면, 표본의 크기가 너무 작으면 오차가 큰 정보를 얻을 수 있다. 따라서 정해진 오차한계를 만족하는 표본의 크기 결정이 필요하다.

모평균 추정에서 표본의 크기 결정

모평균 μ\mu를 추정하는 데 오차한계를 정해진 값 dd보다 작게 하고자 할 때 표본의 크기를 결정해보자. 모분산 σ2\sigma^2을 알 때, 모평균 μ\mu100(1α)%100(1-\alpha)\% 신뢰구간

X±zα/2σn\overline{X} \pm z_{\alpha/2}\frac{\sigma}{\sqrt{n}}

에서

zα/2σnz_{\alpha/2}\frac{\sigma}{\sqrt{n}}

는 모수와 추정량 차이의 100(1α)%100(1-\alpha)\% 오차한계를 의미한다. 따라서 오차한계

zα/2σnz_{\alpha/2}\frac{\sigma}{\sqrt{n}}

dd보다 작다고 하고 nn에 대해 풀면 표본의 크기를 결정할 수 있다.

모평균 추정에서 표본의 크기 결정

모분산 σ2\sigma^2을 아는 경우, 모평균 μ\mu에 대한 표본의 크기 nn의 범위는 다음과 같다.

n>(zα/2σd)2n> \left( z_{\alpha/2}\frac{\sigma}{d} \right)^2

모표준편차 σ\sigma를 알고 있는 경우, 위 식을 이용하면 되지만 실제로는 σ\sigma를 모르는 경우가 대부분이다. 이럴 때는 예비표본을 추출한 후 σ\sigma를 추정하여 표본의 크기 nn을 구할 수 있다.


모비율 추정에서 표본의 크기 결정

모비율 pp100(1α)%100(1-\alpha)\% 신뢰구간

p^±zα/2p^(1p^)n\hat{p} \pm z_{\alpha/2} \sqrt{ \frac{\hat{p}(1-\hat{p})}{n} }

에서 오차한계

zα/2p^(1p^)nz_{\alpha/2} \sqrt{ \frac{\hat{p}(1-\hat{p})}{n} }

을 정해진 값 dd보다 작다고 두고, nn에 대해 풀면 다음과 같다.

n>p^(1p^)(zα/2d)2n> \hat{p}(1-\hat{p}) \left( \frac{z_{\alpha/2}}{d} \right)^2

위 식은 p^\hat{p}이 사전조사나 과거 경험으로 추정한 pp값과 일치하리라 기대할 수 있을 경우에만 사용할 수 있다. 그러나 pp값에 대한 사전정보가 전혀 없을 때는 p^(1p^)\hat{p}(1-\hat{p})의 최대변동인

p^=12\hat{p}=\frac{1}{2}

일 때 주로 위 식을 사용한다.

모비율 추정에서 표본의 크기 결정

모비율 pp에 대한 표본의 크기 nn의 범위는 다음과 같다.

n>14(zα/2d)2n> \frac{1}{4} \left( \frac{z_{\alpha/2}}{d} \right)^2

0개의 댓글