Microsoft MVP성태의 닷넷 이야기
글쓴 사람
정성태 (techsharer at outlook.com)
홈페이지
첨부 파일
 

"Probabilistic Programming and Bayesian Methods for Hackers" 예제 코드 실행 방법

얼마 전 트위터에서 본 내용입니다.

확률이론과 베이즈추론법 프로그래밍에 관한 전자책 Bayesian Methods for Hackers. Python기반.
; https://twitter.com/sjoonk/status/344658745562914816

PDF로 다운로드 받아서 보면 다음과 같이 책에 직접 실행해 볼 수 있는 파이썬 코드가 들어 있습니다.

bayesian_python_code_0.png

윈도우에서 쉽게 이를 테스트 해보려면 "EPD(Enthought Python Distribution)"를 다운로드하시면 됩니다.

Enthought Python Distribution Free
; https://www.enthought.com/products/epd/free/

Download Canopy 1.0, 32-bit for Windows
; https://www.enthought.com/downloads/

위의 무료 버전을 설치하고 실행하면 다음과 같은 화면이 나옵니다.

bayesian_python_code_1.png

Canopy 도구에서는 다행히 "Probabilistic Programming and Bayesian Methods for Hackers" 책에서 요구하는 구성 요소(enstaller, ipython)를 미리 내장하고 있기 때문에 더 이상 별도로 다운로드는 하지 않아도 됩니다. 단지, "Package Manager"를 실행해서 "Updates" 항목에 새로 업데이트 받을 것이 있다면 기분상 해주시면 됩니다.

자... 이제 Editor 버튼을 누르고, 우측의 "In" 명령 프롬프트에서 차례로 본문의 코드를 입력해 주면,

%pylab inline

figsize( 11, 9)

import scipy.stats as stats
dist = stats.beta
n_trials = [0,1,2,3,4,5,8,15, 50, 500]
data = stats.bernoulli.rvs(0.5, size = n_trials[-1] )
x = np.linspace(0,1,100)

for k, N in enumerate(n_trials):
    sx = subplot( len(n_trials)/2, 2, k+1)
    plt.xlabel("$p$, probability of heads") if k in [0,len(n_trials)-1] else None
    plt.setp(sx.get_yticklabels(), visible=False)
    heads = data[:N].sum()
    y = dist.pdf(x, 1 + heads, 1 + N - heads )
    plt.plot( x, y, label= "observe %d tosses,\n %d heads"%(N,heads) )
    plt.fill_between( x, 0, y, color="#348ABD", alpha = 0.4 )
    plt.vlines( 0.5, 0, 4, color = "k", linestyles = "--", lw=1 )
    leg = plt.legend()
    leg.get_frame().set_alpha(0.4)
    plt.autoscale(tight = True)

plt.suptitle( "Bayesian updating of posterior probabilities",
y = 1.02,
fontsize = 14);
plt.tight_layout()

다음과 같이 실행이 되는 것을 확인할 수 있습니다.

bayesian_python_code_2.png

그런데, 이걸 코드 파일로 해서 실행해 볼 수는 없을까요? 내용을 다음과 같이 다소 바꿔주면 됩니다. ^^ (휴~~~ 해당하는 모듈 찾느라 고생했네요. ^^)

import scipy.stats as stats
import matplotlib.pyplot as plt
import numpy as np

plt.figsize(11, 9)

dist = stats.beta
n_trials = [0,1,2,3,4,5,8,15, 50, 500]
data = stats.bernoulli.rvs(0.5, size = n_trials[-1] )
x = np.linspace(0,1,100)

for k, N in enumerate(n_trials):
    sx = plt.subplot( len(n_trials)/2, 2, k+1)
    plt.xlabel("$p$, probability of heads") if k in [0,len(n_trials)-1] else None
    plt.setp(sx.get_yticklabels(), visible=False)
    heads = data[:N].sum()
    y = dist.pdf(x, 1 + heads, 1 + N - heads )
    plt.plot( x, y, label= "observe %d tosses,\n %d heads"%(N,heads) )
    plt.fill_between( x, 0, y, color="#348ABD", alpha = 0.4 )
    plt.vlines( 0.5, 0, 4, color = "k", linestyles = "--", lw=1 )
    leg = plt.legend()
    leg.get_frame().set_alpha(0.4)
    plt.autoscale(tight = True)
    
plt.suptitle( "Bayesian updating of posterior probabilities",y = 1.02,fontsize = 14);
plt.tight_layout()

아마 책의 나머지 예제도 위와 같은 규칙을 적용하면 코드 파일로 저장해서 실행할 수 있을 것입니다. 다음은 실제로 실행된 화면입니다. ^^

bayesian_python_code_3.png

그나저나... 정작 책은 언제 다 읽어볼런지...? ^^

참고로, 닷넷의 경우 Bayesian 추론 관련해서 Infer.NET이라는 라이브러리가 있습니다.

Infer.NET
; http://research.microsoft.com/en-us/um/cambridge/projects/infernet/





[이 글에 대해서 여러분들과 의견을 공유하고 싶습니다. 틀리거나 미흡한 부분 또는 의문 사항이 있으시면 언제든 댓글 남겨주십시오.]







[최초 등록일: ]
[최종 수정일: 1/28/2022]

Creative Commons License
이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.
by SeongTae Jeong, mailto:techsharer at outlook.com

비밀번호

댓글 작성자
 




... 181  182  183  184  185  186  187  [188]  189  190  191  192  193  194  195  ...
NoWriterDateCnt.TitleFile(s)
303정성태7/20/200618800Team Foundation Server: 12. 사용자 계정 재생성에 따른 Version Control 영향
302정성태8/21/200620509Team Foundation Server: 11. TFS Team Build와 VC++ Project 설정
299정성태7/23/200620561개발 환경 구성: 5. VMWare - VM 생성 화면 캡쳐
300정성태7/15/200625469    답변글 개발 환경 구성: 5.1. VMWare 오류 유형 - The handle is invalid.
301정성태7/18/200620426    답변글 개발 환경 구성: 5.2. VMWare - 사용 후기.
298정성태7/14/200620433개발 환경 구성: 4. VMWare Server를 64bit 운영체제에 설치 시 주의 사항 [2]
296정성태7/10/200629800.NET Framework: 73. [ASP.NET] HTC(DHTML Control Behavior)를 WebResource.axd로 제공하는 방법 [3]
295정성태7/1/200623035VC++: 25. Microsoft National Language Support Downlevel APIs 1.0 사용 방법파일 다운로드1
294정성태6/30/200619567.NET Framework: 72. XSDObjectGen.EXE 기능 개선
293정성태6/29/200621445Team Foundation Server: 10. TFS 버전 컨트롤(TFVC)에 참여시킨 프로젝트의 로컬 경로를 옮기는 방법
290정성태6/26/200618772Team Foundation Server: 9. HTTPS를 통한 Team Server 접근 - 두 번째 이야기 [1]
291정성태6/26/200620043    답변글 Team Foundation Server: 9.1. [선행 작업] HTTPS 를 통한 Team Server 접근 - 두번째 이야기 [1]
292정성태6/26/200619828    답변글 Team Foundation Server: 9.2. TF30177 오류 발생
307정성태8/3/200621717    답변글 Team Foundation Server: 9.3. Team Server 접근 이름을 바꾸는 방법 [1]파일 다운로드1
308정성태2/18/200721364        답변글 Team Foundation Server: 9.4. Team Server HTTPS 접근 완료. ^^ [1]
288정성태6/26/200635979오류 유형: 10. error MSB6006: "aspnet_merge.exe" exited with code 1
286정성태6/23/200624097웹: 4. 웹 사이트 식별자(Identifier) 값 변경
285정성태6/20/200624401오류 유형: 9. [TFS] Report 관련 서비스를 조회할 때 rsErrorImpersonatingUser 오류 메시지 발생 [1]
284정성태6/19/200622206VS.NET IDE: 40. FxCop - IDE 에서 제공해 주는 SuppressMessage 코드
283정성태1/19/200723311Team Foundation Server: 8. 소스 세이프에서 TFS SourceControl 로 마이그레이션 [2]
279정성태12/27/200628819개발 환경 구성: 3. VS.NET 원격 디버깅 [1]
280정성태6/12/200628071    답변글 개발 환경 구성: 3.1. VS.NET 2003 원격 디버깅 설정
281정성태8/11/200629570    답변글 개발 환경 구성: 3.2. VS.NET 2005 원격 디버깅 설정
315정성태8/11/200630191        답변글 개발 환경 구성: 3.3. VS.NET 2005 원격 디버깅 설정 - ASP.NET F5 디버깅
278정성태6/11/200626599오류 유형: 8. [Outlook] 0x8004011D 에러 - "Exchange over the Internet" 환경
276정성태6/7/200619966Team Foundation Server: 7. 외부 빌드 머신 구성
... 181  182  183  184  185  186  187  [188]  189  190  191  192  193  194  195  ...