지수족은 다음과 같은 형태의 확률 분포를 말한다.
지수족의 원소로는 정규분포, 베르누이 분포, 지수분포, 베타분포, 감마분포, 푸아송분포, 카이제곱분포 등등 우리가 알고 있는 대부분의 분포 그러나 student-T분포는 지수족이 아니다.
또한, 시행횟수 N이 고정된 이항분포들의 집합은 지수족의 원소가 된다. n이 고정되어야 하고 n과 p가 모두 모수이면 exponential family가 안 된다.
https://youtu.be/3ZFNo43Bt-o?si=WZZCIO3DYua5vx6h
이항분포
이항분포가 지수족에 속하는지를 증명해보자.
\[f(x) = (nCx) p^{x}\times (1-p)^{n-x}\]
1)take the log of both sides
\[ \log{f(x)} = \log(nCx) + x\log{p} +(n-x)\log{(1-p)} \]
2) take the exponential of both sides
\[ f(x) = \exp{(x\log{p}-x\log{(1-p)}+n\log{(1-p)}+\log{(nCx)})}\]
\[= \exp{(x \log{(p/(1-p)}+n\log{(1-p)}+\log{(nCx)})} \]
Compare to -> \[f(y) = \exp((\theta(y) - b\theta)/a(\phi) +c(y;\phi)) \]
let \[ \theta = \log{p/(1-p)}\], then \[e^{\theta} = p/(1-p)\]
\[p = e^{\theta}/(1+e^{\theta}), 1-p = 1/(1+e^{\theta}) = (1+e^{\theta})^{-1}\]
\[n \log{(1-p)}= n\log{((1+e^{\theta})^{-1})} = -n\log{(1+e^{\theta})}\]
\[f(x) = \exp(x\theta- n\log{(1+e^{\theta})}+\log {(nCx)}) \]
\[b(\theta) = n\log{(1+e^{\theta})} , a(\phi) = 1, c(x;\phi) = \log{(nCx)} \]
다음으로 정규분포가 지수족에 속하는지를 증명해보자.
정규분포
\[f(x|\mu,\sigma^{2}) = \frac{1}{\sqrt{2\pi\sigma^{2}}}\exp(\frac{-(x-\mu)^{2}}{2\sigma^{2}})\]
1) 양변에 로그 취하기
\[\log{f(x|\mu,\sigma^{2})} = -\frac{1}{2}\log{2\pi\sigma^{2}}-(\frac{(x-\mu)^{2}}{2\sigma^{2}})\]
2) 양변에 지수 취하기
\[f(x|\mu,\sigma^{2}) = \exp(-\frac{1}{2}\log{2\pi\sigma^{2}}-\frac{(x-\mu)^{2}}{2\sigma^{2}})\]
\[ f(x|\mu, \sigma^{2}) = \exp{(\frac{-x^{2}+2x\mu-\mu^{2}}{2\sigma^{2}}}- \frac{1}{2}\log{(2\pi\sigma^{2}))} \]
\[ = \exp(\frac{2x\mu-\mu^{2}}{2\pi^{2}} - \frac{x^{2}}{2\sigma^{2}}-\frac{1}{2}\log{(2\pi\sigma^{2})}) \]
\[=\exp(\frac{x\mu-\frac{1}{2}\mu^{2}}{\sigma^{2}} - \frac{ x^{2}}{2\sigma^{2}} -\frac{1}{2}log{(2\pi\sigma^{2})}) \]
\[\theta = \mu, b(\theta) = \frac{1}{2}\theta^{2}, a(\phi)= \sigma^{2}\]
https://youtube.com/playlist?list=PLvCaViUWnpSn1YmTT6_3P_At1UapqaLGy&si=XVJ17mnAZoyd-M9m
'Stat' 카테고리의 다른 글
Factor Analysis 요인분석 R로 해보기 (0) | 2024.10.27 |
---|---|
충분통계량 (1) | 2023.12.10 |
Efficient Estimator 유효추정량 (0) | 2023.11.28 |
Unbiased Estimator 불편추정량 (0) | 2023.11.28 |