📊 Machine Learning · Latent Variable Models · Unsupervised Learning
홈 aistudy 메인

기댓값 최대화 알고리즘 (Expectation-Maximization Algorithm : EM)

관측되지 않는 잠재 변수(Latent Variable) 나 결측치(Missing Data) 가 존재하는 확률 모델에서, E단계(기댓값 계산 )와 M단계(우도 극대화) 의 교대 반복을 통해 모수의 최대우도추정치(MLE) 또는 최대사후확률추정치(MAP)를 점근적으로 수렴해내는 통계적 기계학습의 핵심 메타 알고리즘 .... EM 알고리즘의 개념은 통계학에서 자주 쓰이는 "앞면이 나올 확률이 서로 다른 두 동전 A, B" 예시로 쉽게 이해할 수 있습니다.

EM의 접근법:
1. 초기화: 동전 A와 B의 앞면 확률을 대충 가정합니다 (예: A는 0.6, B는 0.4).
2. E 단계: 가정한 확률을 바탕으로, 관측된 앞/뒤 결과가 '동전 A에서 나왔을 확률'과 '동전 B에서 나왔을 확률'을 각각 계산(기댓값 할당)합니다.
3. M 단계: 가중치(확률)가 부여된 데이터를 바탕으로 동전 A와 B의 앞면 확률을 다시 더 정확하게 업데이트합니다.
이를 반복하면 점차 실제 두 동전의 숨겨진 확률로 수렴합니다.

📖

1. 원전 기록 및 고전적 정의 보존

aistudy 원전 기록 보존

EM 알고리즘은 불완전 정보 (missing data) 에 대한 Posterior 분포의 해석적 표현을 얻을 수 있으며, 완전 정보 (complete data) 의 최대우도 추정 (Maximum Likelihood Estimation) 의 해석적 표현이 가능할 때는 언제나 강력한 계산도구로 사용할 수 있으며, 추정치 (estimate) 가 안정적이다.

Wikipedia : Expectation-maximization algorithm : 통계 (Statistics) 에서 EM 알고리즘은 보이지 않는 잠재 변수 (latent variable) 에 의존하는 확률모델에서 모수 (parameter) 들의 최대우도 추정치 (maximum likelihood estimates of parameters) 를 찾고자하는 알고리즘이다.

즉 expectation (E) 단계에서는 잠재 변수의 기대치를 계산하고, maximization (M) 단계에서는 주어진 데이터와 기대치가 부여된 잠재 변수를 이용하여 모수들의 최대우도 추정치를 계산한다.

"Expectation-maximization" 은 특별한 알고리즘이 아니라 관련된 알고리즘의 종류에서 따온 이름이다 ; EM 은 특별한 알고리즘을 발명하기위해 사용되는 하나의 수단 (recipe) 또는 메타알고리즘 (meta-algorithm) 이다. Baum-Welch algorithm 은 hidden Markov models 에 응용된 EM 알고리즘의 한 예이다. 또다른 예는 mixture density model를 fitting 하기위한 EM 알고리즘이다.

🔍 잠재 변수 (Latent Variable)

실제 관측 데이터 이면에 존재하지만 직접 측정할 수 없는 숨은 변수(예: 클러스터 소속 라벨, HMM의 은닉 상태 등)를 통계적으로 다룹니다.

🔄 E-M 교대 최적화 (Alternating)

잠재 변수의 사후 기대치를 계산하는 E단계와, 이를 고정하고 모수를 극대화하는 M단계를 번갈아 수행하여 우도를 단조 증가시킵니다.

🛡️ 안정적 수렴성 (Convergence)

매 반복(Iteration)마다 관측 데이터의 로그 우도가 감소하지 않음을 수학적으로 보장하며 국소 최적해(Local Maximum)에 안정적으로 도달합니다.

📐

2. EM 알고리즘의 기본 개념과 수학적 작동 원리

(1) '닭과 달걀'의 딜레마 — 왜 직접 풀 수 없는가?

관측된 데이터 y 와 보이지 않는 잠재 변수 z 가 결합된 완전 정보(Complete Data)를 상정해 봅시다.

(2) E단계와 M단계의 수학적 정식화

EM 알고리즘은 임의의 초기 모수 θ_0 에서 출발하여, 모수 추정치를 θ_1, θ_2, ..., θ_N 으로 점진적으로 갱신합니다.

/* 1. 완전정보 로그우도의 기댓값 함수 Q(θ | θ_n) 정의 */
Q(θ | θ_n) = E_(z | y, θ_n) [ log p(y, z | θ) ] = ∑_z p(z | y, θ_n) × log p(y, z | θ)

/* 2. M단계 모수 갱신 (Q 함수를 극대화하는 모수 선택) */
θ_(n+1) = arg max_θ Q(θ | θ_n) = arg max_θ ∑_z p(z | y, θ_n) × log p(y, z | θ)

[원전 수식 표기 보존]

θ_{n+1} = \arg\max_{\theta} \sum_z p(z \,|\, y, \theta_n) \times \log p(y, z \,|\, \theta)
Q(\theta) = \sum_z p(z \,|\, y, \theta_n) \times \log p(y, z \,|\, \theta)
💡 expectation (E) 단계에 대한 원전의 정밀한 통찰 :
"expectation (E) 단계에 대해 말하자면 조금 잘못된 명칭 (misnomer) 이다. E단계에서 실제로 계산되는 것은 함수 Q 의 고정된 (fixed), 데이터 의존 (data-dependent) 모수들(사후확률)이다. 일단 Q 의 모수들이 알려지면, Q 함수는 완전히 결정되고 maximization (M) 단계에서 최대화된다."

(3) 젠센 부등식(Jensen's Inequality)과 하한선(ELBO)의 기하학적 의미

EM 알고리즘의 본질은 로그 우도 함수 log p(y | θ) 의 '증거 하한(Evidence Lower Bound : ELBO)' 곡면을 계속해서 만들어 들어 올리는 과정입니다.

📈 E단계의 기하학적 의미

현재의 모수 추정치 θ_n 에서 실제 로그 우도 곡선에 완벽히 접하는 보조 함수(Lower Bound)를 구합니다. 이 점에서 실제 우도와 하한 함수의 차이(KL Divergence)는 정확히 0이 됩니다.

🚀 M단계의 기하학적 의미

E단계에서 만든 보조 하한 곡선의 정점(Maximum)으로 모수를 이동시켜 새로운 θ_(n+1) 을 찾습니다. 하한이 높아졌으므로 실제 우도 또한 반드시 같거나 높아집니다.

📌 단조 수렴성 증명 (Dempster, Laird, Rubin, 1977) :
EM 알고리즘의 매 반복마다 관측 데이터의 로그 우도는 감소하지 않으며, 단조 증가(Monotonic Increase)합니다 : log p(y | θ_(n+1)) ≥ log p(y | θ_n). 따라서 유계된 우도 함수 하에서 EM 알고리즘은 언제나 국소 최적해(Local Maximum)에 안정적으로 수렴합니다.
🧩

3. 대표적 응용 : 가우시안 혼합 모델(GMM)과 K-평균 군집화

EM 알고리즘이 가장 널리 쓰이는 전형적인 사례는 비지도 학습의 대표 모델인 가우시안 혼합 모델(Gaussian Mixture Model : GMM)입니다.

비교 항목 K-평균 군집화 (K-Means Clustering) GMM + EM 알고리즘 (Soft Clustering)
군집 할당 방식 하드 군집화 (Hard Clustering)
데이터가 특정 클러스터에 100% 속함 (0 또는 1)
소프트 군집화 (Soft Clustering)
각 클러스터에 속할 사후 확률(소속도, Responsibility) 부여
클러스터 형태 원형 / 구형 (Spherical)
모든 클러스터의 분산과 방향이 동일하다고 가정
타원형 (Elliptical)
각 클러스터마다 고유한 평균벡터와 공분산 행렬을 가짐
알고리즘 연관성 EM 알고리즘의 극단적인 특수 형태
(분산이 0으로 수렴할 때의 하드 EM)
일반화된 확률적 군집화 프레임워크
E단계 동작 가장 가까운 군집 중심점에 점들을 배정 현재 모수로 각 점이 각 가우시안에 속할 사후 확률 계산
M단계 동작 배정된 점들의 산술 평균으로 중심점 재계산 소속도 가중치를 반영하여 평균, 공분산, 혼합비율 갱신
🚀

4. EM 알고리즘이 AI · 기계학습 · 데이터마이닝 발전에 기여한 점

1977년 뎀스터(Dempster), 레어드(Laird), 루빈(Rubin)에 의해 체계화된 EM 알고리즘은 통계학을 넘어 인공지능, 기계학습, 데이터마이닝의 패러다임을 바꾼 혁명적 도구였습니다.

🎬

추천 강의 : The Expectation-Maximization (EM) Algorithm 원리 해설

잠재 변수 모델과 가우시안 혼합 모델(GMM)을 통해 EM 알고리즘의 E단계와 M단계 수렴 원리를 직관적이고 명쾌하게 설명하는 강의

📚

5. 원전 참고 문헌 및 국내외 학술 논문 보존

📖 해외 고전 학술 문헌 (Classic Papers) :

📄 국내 연구 및 학술 아카이브 논문 (Korean Papers) :

🌐 외부 권위 사전 및 최적화 기법 링크 :

🌐

6. aistudy 수학 · 기계학습 지식 네트워크

EM 알고리즘과 유기적으로 연결되는 aistudy 핵심 지식 문서들을 탐색할 수 있습니다.