베이즈정리와 정규분포의 곱Bayesian theorem and normal distribution multiplication
베이지안 선형회귀bayesian linear regression를 공부하다 보면 계수벡터 \(w\)에 대한 사후 확률분포posterior를 구하게 된다. 이때 데이터의 분포를 정규분포로 가정한 가능도와 \(w\)에 대한 공액 사전분포conjugate prior의 곱을 완전제곱completing the square 형태로 바꾸는 테크닉을 사용하여 사후 확률분포의 평균과 분산을 간단히 구하는 과정을 거치게 된다.
정규분포를 계수벡터의 사전분포로 선택하면 그것의 사후분포도 정규분포가 되는 성질을 이용한 것인데 이 글에서는 베이즈룰을 완전히 손으로 계산하여 일변수 정규분포에 대해서 \(\frac{p(D|\mu,\sigma^2)p(\mu)}{p(D)}\)가 정규분포가 되는지 확인해본다.
이를 위해 정규분포의 곱을 정리하여 간단한 프로그램으로 곱의 의미를 확인 해보고 그 결과를 이용하여 가능도와 공액 사전분포를 곱하여 정리하였다.
그후 \(p(D,\mu)\)를 \(\mu\)에 대해 주변화시키는 과정인 적분을 직접하여 베이즈정리에 의해 구해지는 사후 확률분포가 정확히 정규분포가 되는지 확인하였다. 매우 긴 계산과정으로 인해 많은 자료에서 생략된 부분이라 구체적인 계산이 궁금한 분들에게 도움이 되리라 생각한다.
전체 글은 jupyter notebook으로 작성되어 있어서 아래 링크를 통해 nbviewer로 공유한다.