충분통계량

오류 제보
통계 수치해석 마지막 수정: 2026-08-11 04:28:41

○ 공개 검수 대기 표시 기준 보기

1. 개요

충분통계량
Sufficient Statistic
정의T(X) 가 주어졌을 때 X 의 조건부분포가 θ 에 의존하지 않음
판정피셔-네이만 분해정리
도입R. A. Fisher (1920, 1922)
핵심 정리라오-블랙웰 · 레만-셰페 · 바수 · 피트먼-쿠프먼-다르무아
대표 예정규: (Σx, Σx²) · 포아송: Σx · 균등 U(0,θ): max x
실무적 의미데이터를 고정 크기로 압축해도 추론이 손해 보지 않음

충분통계량(sufficient statistic)은 표본 X=(X1,…,Xn)X = (X_1,\dots,X_n) 을 요약한 함수 T(X)T(X) 중에서, TT 의 값을 알고 나면 원자료가 모수 θ\theta 에 대해 더 이상 아무 말도 해 주지 않는 것을 말한다. 형식적으로는 조건부분포

Pθ(X∈A∣T(X)=t)가 모든 A, t 에 대해 θ 에 의존하지 않는다P_\theta\big(X \in A \mid T(X) = t\big) \quad\text{가 모든 } A,\, t \text{ 에 대해 } \theta \text{ 에 의존하지 않는다}

는 조건으로 정의한다. 사람 말로 옮기면 이렇다. T=tT=t 만 받아 적고 데이터를 태워도, 그 tt 로부터 원자료와 통계적으로 구분 불가능한 가짜 데이터를 시뮬레이션해 낼 수 있다. 태운 것이 아깝지 않다는 게 “충분”의 의미다.

동전을 100번 던져 앞면이 37번 나왔다면, 앞면 확률 pp 를 추론하는 데 필요한 것은 37이라는 숫자뿐이고 “몇 번째에 앞이 나왔는가”라는 순서 정보는 통째로 버려도 된다. 실제로 ∑Xi=37\sum X_i = 37 을 조건으로 걸면 가능한 배열들은 pp 와 무관하게 전부 같은 확률을 갖는다. 이것이 충분성의 원형이며, 피셔가 1920년대에 이 개념을 만들면서 현대 추정론의 문법이 시작됐다.1

2. 피셔-네이만 분해정리

정의를 그대로 확인하려면 조건부분포를 계산해야 해서 성가시다. 실무에서 쓰는 것은 분해정리다. TT 가 충분할 필요충분조건은 우도가

f(x∣θ)=g(T(x), θ)  h(x)f(x \mid \theta) = g\big(T(x),\ \theta\big)\; h(x)

로 모수를 포함한 부분과 데이터만의 부분으로 쪼개지는 것이다. θ\theta 가 데이터를 오직 T(x)T(x) 를 통해서만 만나면 충분하다는, 거의 동어반복처럼 보이지만 판정 도구로는 압도적으로 편한 정리다.

이 정리를 지수족에 적용하면 왜 그 분포족이 특별한지가 한 줄로 나온다. 표준형 p(x∣η)=h(x)exp⁡(η⊤T(x)−A(η))p(x\mid\eta)=h(x)\exp(\eta^\top T(x) - A(\eta)) 를 nn 개 곱하면

p(x1:n∣η)=(∏ih(xi))exp⁡ ⁣(η⊤ ⁣∑iT(xi)−nA(η))p(x_{1:n}\mid\eta) = \Big(\textstyle\prod_i h(x_i)\Big)\exp\!\Big(\eta^{\top}\!\sum_i T(x_i) - nA(\eta)\Big)

이라 ∑iT(xi)\sum_i T(x_i) 가 곧바로 충분통계량이다. 표본이 몇 개든 차원이 늘지 않는다.

몇 가지 표준 예시.

  • 베르누이 Bern(p)\mathrm{Bern}(p): ∑Xi\sum X_i
  • 포아송 Pois(λ)\mathrm{Pois}(\lambda): ∑Xi\sum X_i
  • 정규 N(μ,σ2)N(\mu,\sigma^2), 둘 다 미지: (∑Xi, ∑Xi2)\big(\sum X_i,\ \sum X_i^2\big)
  • 균등 U(0,θ)U(0,\theta): max⁡iXi\max_i X_i
  • 균등 U(θ−12, θ+12)U(\theta-\tfrac12,\ \theta+\tfrac12): (min⁡iXi, max⁡iXi)\big(\min_i X_i,\ \max_i X_i\big)
  • 코시 위치모수: 순서통계량 전체. 압축 불가.

마지막 두 개가 중요하다. 균등분포는 지수족이 아닌데도 1차원(또는 2차원) 충분통계량을 갖는다 — 지지집합이 θ\theta 에 의존하는 덕에 정리의 가정을 빠져나가기 때문이다. 코시는 정반대로, 표본 1억 개를 요약할 방법이 원리적으로 없다.

3. 최소충분통계량

충분통계량은 유일하지 않다. 극단적으로 T(X)=XT(X) = X 자체가 항상 충분하고(아무것도 안 버렸으니 당연히), 충분통계량의 가역변환도 여전히 충분하다. 그러니 의미 있는 질문은 “가장 많이 압축한 것은 무엇인가”다.

최소충분통계량(minimal sufficient statistic)은 다른 모든 충분통계량의 함수로 표현되는 충분통계량이다. 레만과 셰페가 준 실용적 판정법은 우도비를 보는 것이다.

f(x∣θ)f(y∣θ) 가 θ 에 무관  ⟺  T(x)=T(y)\frac{f(x\mid\theta)}{f(y\mid\theta)} \text{ 가 } \theta \text{ 에 무관} \iff T(x) = T(y)

즉 “두 데이터셋이 모든 모수에서 같은 우도 곡선을 그린다”는 관계로 표본공간을 분할하면, 그 분할이 최소충분통계량이다. 이 관점은 우도원리와 곧장 연결된다. 추론이 우도 곡선의 모양에만 의존해야 한다면, 최소충분통계량 밖의 정보는 정의상 볼 필요가 없다.

정규분포에서 (∑Xi,∑Xi2)(\sum X_i, \sum X_i^2) 가 최소충분이고, 표본평균 하나만으로는 σ2\sigma^2 정보를 버리므로 충분하지 않다. 반대로 σ2\sigma^2 가 기지라면 ∑Xi\sum X_i 만으로 최소충분이다. 무엇이 미지인가에 따라 압축률이 달라진다는 점이 종종 헷갈리는 지점.

4. 라오-블랙웰 정리 — 조건부 기댓값은 손해를 안 본다

충분성이 단순한 정리 놀음이 아닌 이유가 이 정리다. δ(X)\delta(X) 가 임의의 추정량이고 TT 가 충분통계량이면

δ∗(T)=E[δ(X)∣T]\delta^{*}(T) = \mathbb{E}\big[\delta(X) \mid T\big]

은 여전히 추정량이며(충분성 덕에 θ\theta 가 들어가지 않는다 — 이게 결정적이다), 볼록 손실함수에 대해

R(θ,δ∗)≤R(θ,δ)∀ θR(\theta, \delta^{*}) \le R(\theta, \delta) \qquad \forall\,\theta

를 만족한다. 제곱오차 손실이면 젠센 부등식 한 방이고, 분산 분해 Var(δ)=Var(E[δ∣T])+E[Var(δ∣T)]\mathrm{Var}(\delta) = \mathrm{Var}(\mathbb{E}[\delta\mid T]) + \mathbb{E}[\mathrm{Var}(\delta\mid T)] 에서 두 번째 항이 통째로 사라지는 것이 개선분이다. 편향은 반복기댓값 법칙으로 그대로 보존된다.

실용적 번역: 엉성한 불편추정량 하나를 아무렇게나 만든 다음 충분통계량으로 조건부 기댓값을 취하면 무조건 더 좋아진다. 라오(1945)와 블랙웰(1947)이 독립적으로 얻은 결과이며, 오늘날 몬테카를로 분산감소 기법인 “라오-블랙웰화”의 어원이기도 하다. 입자 필터에서 선형-가우시안 부분을 칼만 필터로 해석적으로 적분하고 나머지만 입자로 다루는 Rao-Blackwellised particle filter가 정확히 이 정리를 몬테카를로 방법에 적용한 것이다.2

5. 완비성과 레만-셰페 — UMVU로 가는 길

라오-블랙웰은 “더 좋아진다”만 말하지 “이게 최선이다”를 말하지 않는다. 시작 추정량이 다르면 다른 곳에 도착할 수도 있다. 이 마지막 구멍을 막는 것이 완비성(completeness)이다.

Eθ[g(T)]=0    ∀θ⟹g(T)=0  (거의 확실히)\mathbb{E}_\theta\big[g(T)\big] = 0 \;\;\forall\theta \quad\Longrightarrow\quad g(T) = 0 \ \ \text{(거의 확실히)}

를 만족하면 TT 를 완비라 한다. 풀어 쓰면 ”TT 의 함수 중에 항등적으로 0을 불편추정하는 것은 0뿐”이라는 뜻이고, 결국 TT 의 함수로서 어떤 양의 불편추정량은 많아야 하나라는 유일성 진술이다.

여기서 레만-셰페 정리가 나온다. TT 가 완비충분통계량이고 δ∗(T)\delta^*(T) 가 τ(θ)\tau(\theta) 의 불편추정량이면, δ∗\delta^* 는 유일한 UMVU 추정량(uniformly minimum variance unbiased)이다. 실무 레시피는 두 줄로 끝난다.

  1. 완비충분통계량 TT 를 찾는다(대개 지수족이면 자동).
  2. TT 의 함수 중 불편인 것을 하나 찾는다. 그게 답이다.

정규분포의 Xˉ\bar X 와 S2S^2(분모 n−1n-1)가 각각 μ\mu, σ2\sigma^2 의 UMVU인 것이 이 레시피의 결과다. U(0,θ)U(0,\theta) 에서는 max⁡Xi\max X_i 가 완비충분이고 E[max⁡Xi]=nn+1θ\mathbb{E}[\max X_i] = \frac{n}{n+1}\theta 이므로 n+1nmax⁡Xi\frac{n+1}{n}\max X_i 가 UMVU다 — 표본평균의 2배(=2Xˉ=2\bar X, 역시 불편)보다 분산이 압도적으로 작고, 심지어 nn 이 커지면 오차가 1/n1/n 로 줄어 크라메르-라오 하한의 1/n1/\sqrt n 을 “뚫는” 것처럼 보인다. 물론 뚫은 게 아니라 정칙조건 밖에서 논 것.3

바수 정리도 완비성의 부산물이다. 완비충분통계량은 임의의 부수통계량(ancillary statistic, 분포가 θ\theta 에 무관한 통계량)과 독립이다. 정규표본에서 Xˉ\bar X 와 S2S^2 의 독립성을 계산 없이 세 줄로 증명할 수 있고, 이 독립성이 t-분포 유도의 출발점이다.

완비성이 깨지는 사례도 알아 둘 가치가 있다. 굽은 지수족 N(θ,θ2)N(\theta,\theta^2) 에서는 (∑Xi,∑Xi2)(\sum X_i, \sum X_i^2) 가 최소충분이지만 완비가 아니다 — 두 통계량이 같은 θ\theta 를 서로 다른 방식으로 추정하므로, 그 차이가 0을 불편추정하는 비자명한 함수가 된다. 완비가 아니면 UMVU 유일성 논증이 통째로 무너진다. 최소충분 ≠ 완비임을 기억해 두면 시험에서 한 문제 건진다.

6. 피트먼-쿠프먼-다르무아 정리

충분통계량 이야기의 정점은 그 역이다.

지지집합이 모수에 의존하지 않는 정칙 분포족 중에서, 표본 크기 nn 이 커져도 차원이 고정된 충분통계량을 갖는 것은 지수족뿐이다.

다르무아(1935), 쿠프먼(1936), 피트먼(1936)이 거의 동시에 독립적으로 얻었다. 함의가 상당히 잔인하다.

  • 데이터를 유한 개 숫자로 요약해 두고 원자료를 버리는 관행 — 온라인 추정, 스트리밍 집계, 분산 학습에서 각 노드가 부분합만 보내는 방식 — 은 지수족 모형에서만 정확히 정당화된다. 그 밖에서는 근사다.4
  • 코시나 스튜던트 t 같은 두꺼운 꼬리 모형에서는 요약이 원리적으로 불가능하다. 로버스트 통계가 유독 계산량이 많은 이유의 절반은 여기에 있다.
  • 정리의 가정(지지집합 조건)을 빼먹으면 U(0,θ)U(0,\theta) 라는 반례가 즉시 등장한다. 세미나에서 이 정리를 인용할 때 가정을 안 붙이면 반드시 손이 올라온다.

7. 계산과 시뮬레이션에서의 충분성

교과서 밖에서 이 개념이 실제로 일하는 곳들.

  • 기댓값 최대화 알고리즘: 완전자료 로그우도가 충분통계량에 선형이면 E-단계는 “충분통계량의 조건부 기댓값 몇 개”로 붕괴하고, M-단계는 완전자료 MLE에 그 기댓값을 대입하는 것으로 끝난다. 가우시안 혼합 모형이나 은닉 마르코프 모형의 갱신식이 깔끔한 이유가 이것이다.
  • 온라인·스트리밍 추정: 충분통계량 벡터 하나만 들고 다니면 되므로 메모리가 nn 과 무관하다. 미니배치마다 누적합을 갱신하는 확률적 EM, 분산 환경에서 부분합을 합치는 집계 방식이 전부 여기에 기댄다.
  • 필터링의 믿음 상태: 과거 이력 전체를 요약하는 충분통계량이 존재하면 상태공간이 유한 차원으로 닫힌다. 칼만 필터에서 그것이 (평균, 공분산) 두 개이고, 부분관측 마르코프 결정 과정에서는 믿음 상태 전체다.
  • 근사 베이즈 계산(ABC): 우도를 못 쓰는 시뮬레이터 기반 추론에서는 관측과 시뮬레이션 결과를 요약통계량으로 비교한다. 그 요약통계량이 충분하지 않으면 정보 손실이 그대로 사후분포 왜곡으로 나타난다 — ABC가 “근사”인 첫 번째 이유가 요약통계량의 불충분성이다. 지수족이 아닌 복잡한 시뮬레이터에서 충분통계량은 존재하지 않으므로, 이 손실은 원리적으로 피할 수 없고 잘 고른 요약통계량으로 줄이는 수밖에 없다.
  • 신경망의 병목층: 표현학습에서 “정보를 잃지 않는 압축”을 말할 때 참조하는 이상적 기준이 충분성이다. 물론 실제 신경망 표현이 충분하다는 보장은 어디에도 없다.

8. 관련 문서

9. Footnotes

  1. 피셔는 “sufficient”라는 단어를 1922년 논문에서 확정했는데, 같은 논문에서 “efficiency”, “likelihood”, “estimation”까지 한꺼번에 정의했다. 한 편으로 학문 용어집을 새로 쓴 셈. 덕분에 후대 통계학도는 그 논문 하나를 이해하려다 학기를 보낸다. ↩

  2. 이름이 재미있는 것은 라오-블랙웰이 원래 “추정량을 개선하는 정리”인데 몬테카를로 동네에서는 “적분할 수 있는 건 적분해라”라는 실천 지침으로 쓰인다는 점이다. 샘플링으로 때울 수 있다고 다 때우면 분산이 커진다 — 손으로 풀 수 있는 부분은 손으로 푸는 게 이득이라는, 수치해석 전반에 통하는 국룰의 통계학 버전. ↩

  3. 지지집합이 모수에 의존하는 문제를 통칭 비정칙(non-regular) 문제라 하고, 수렴률이 n−1n^{-1} 이나 그보다 빠른 것들이 흔히 튀어나온다. 처음 보면 “공짜로 정밀도가 좋아졌다”는 착각을 하는데, 대신 극한분포가 정규가 아니라서 표준오차·신뢰구간 공식이 전부 못 쓰게 된다. 세상에 공짜는 없다. ↩

  4. 그래서 “우리 파이프라인은 평균과 분산만 저장합니다”라는 문장은 은근히 강한 모형 가정을 깔고 있다. 두꺼운 꼬리가 의심되는 로그 데이터에 이 관행을 적용해 놓고 나중에 분포를 다시 보고 싶어지면, 이미 태워 버린 원자료를 되살릴 방법이 없다. ↩

이 문서가 문제를 푸는 데 도움이 됐나요? 한 번의 클릭이 다음 검수 순서를 정합니다.

의견을 반영할게요. 감사합니다.

다음으로 읽기 최대우도추정 약 5분 더 탐색 →