WinBUGS 베이지안 통계 프로그램

2017. 7. 10. 23:02

WinBUGS 베이지안 통계 프로그램

베이지안을 바탕으로 Markove chain Monte Carlo (MCMC)방법을 이용한 통계 프로그램입니다. 석사 페이퍼쓰면서 이용하는 중인데 자료 찾은걸 바탕으로 정리하고 싶어서 포스팅합니다. exe파일이라서 윈도우에서 설치할 수 있어요.

WinBUGS 설치방법

사이트 : https://www.mrc-bsu.cam.ac.uk/software/bugs/the-bugs-project-winbugs/

설치방법은 동영상보고 참고했습니다.

WinBUGS에서 중요한 점

가끔 오류나고 그러는데 (에러메시지는 제가 정리를 안해둬서요) 관리자모드로 실행해야합니다. R2WinBUGS패키지 이용해서 WinBugs프로그램 불러오는 경우 역시 R을 실행할때도 관리자모드로 실행해야 오류가 안뜹니다.

전 HIV 데이타를 가지고 WinBugs을 이용해 parameters을 구하는 중인데요.

모델을 적고, 이니셜벨류를 지정하고 (혹은 알아서 generate할수도 있습니다) 데이터를 입력하게되면 stat을 이용해서 parameters값을 보거나 trace기능을 통해 계산되는 과정을 볼 수 있습니다.

제가 참고한 사이트 pdf버전 올려드리니 참고하면 될 듯 합니다.

http://homepage.stat.uiowa.edu/~gwoodwor/BBIText/AppendixBWinbugs.pdf

http://www.stats.ox.ac.uk/~cholmes/Courses/BDA/Winbugs/winbugs-help.pdf

R2WinBUGS 패키지

R을 이용해서 WinBUGS를 사용할 수 있고요. 예를들어, 시뮬레이션을 통해 데이터를 생성한 후, parameters의 평균값을 구할때 보통 R을 이용해 코드를 짜고 거기에 WinBUGS를 불러와서 각 데이터의 parameters값을 구할 수 있는데요. 아래 이미지 보면 제 컴퓨터 캡쳐화면이고요. Console보면 7, 8, 9 등... 7번째 8번째 9번째 데이터를 만들어서 거기에 따른 값을 WinBUGS을 통해 구해진걸 보여지는 상황입니다. 200번 구해서 평균값을 구할 예정이고요.

제가 참고한 자료 올립니다.

http://www.dme.ufrj.br/mcmc/R2WinBUGS.pdf

http://plantecology.syr.edu/fridley/bio793/bayes.html

https://www.unc.edu/courses/2010fall/ecol/563/001/docs/lectures/lecture13.htm