figsize(12.5, 5) simulations = mcmc.trace(“bernoulli_sim”)[:].astype(int) print(“시뮬레이션 배열의 형태:”, simulations.shape) plt.title(“사후확률분포의 모수를 사용한 시뮬레이션 데이터셋”) figsize(12.5, 6) for i in range(4): ax = plt.subplot(4, 1, i + 1) plt.scatter(temperature, simulations[1000 * i, :], color=“k”, s=50, alpha=0.6) plt.ylabel(“사고?”) plt.xlabel(“외부 기온(화씨)”)
결과
시뮬레이션 배열의 형태: (10000, 23)