# 상수 설정 p_true = 0.05 # 실제로 이 값은 미지수다. N = 1500 # Ber(0.05)로부터 N개의 베르누이 확률변수를 표본추출한다. # 각 확률변수가 1일 확률은 0.05다. # 이것은 데이터 생성 단계다. occurrences = pm.rbernoulli(p_true, N) print(occurrences) # 파이썬에서 True == 1이고, False == 0이다. print(occurrences.sum())
결과
[False False False …, False False False] 85
관측된 빈도는 다음과 같다.
# Occurrences.mean은 n/N이다. print(“A 사이트에서 관측된 빈도는 얼마인가? %.4f“ % occurrences.mean()) print(“실제 빈도와 같은가? %s“ % (occurrences.mean() == p_true))
결과
A 사이트에서 관측된 빈도는 얼마인가? 0.0567 실제 빈도와 같은가? False