Monte Carlo Integration
- 기대값을 샘플의 평균으로 근사
Importance Sampling
- 효율적인 접근 방법
- sample의 w(s)로 한다.
- 평균은 weighted sample의 평균이다.
- 관심있는 영역이나, 중요한 영역에서 샘플하길 원한다.
- 중요한 영역에 over-weight을 주는 분산으로 부터 샘플 추출
- 왜 이 것을 쓰는가?
- A영역 밖에 랜덤 변수들은 거의 0에 가깝다면, 샘플링 실패 확률이 높아진다.
Gibbs sampling
- 확률적 알고리즘: 결합 확률 분포 => 일련의 표본 생성
- 2개 이상 확률 변수
- MCMC의 대중적인 sampler
- 다 변수
- Bayesian 방법으로 인기가 있다.
- 자주 conditional expression 사용되기 때문
- 느리다.
- 높은 확률을 가지는 쪽에 빠질 수 있다.
- heavy tails나 multiple mode일 경우에 약하다.
Random Generation Number
- randn()
- 표준정규분포(standard normal distribution)
- 평균이 0
- 표준 편차 1
- normrnd(mu,sigma)=> r = randn() * sigma + mu
- http://theclevermachine.wordpress.com/tag/gibbs-sampler/ 참조