問い：ビーチパーティーを開くべきか？

RによるベイズデータAnalysisの基礎

Rasmus Bååth

Data Scientist

問い

7月20日の平均水温はどのくらいか？
_次の_ 20日に水温が18度以上になる確率はどのくらいか？

事後分布

pars

  mu sigma probability
17.5   1.9      0.0001
18.0   1.9      0.0003
18.5   1.9      0.0014
19.0   1.9      0.0043
19.5   1.9      0.0094
20.0   1.9      0.0142
20.5   1.9      0.0151
21.0   1.9      0.0112
21.5   1.9      0.0058
22.0   1.9      0.0021
...    ...       ...

sample_indices <- sample(1:nrow(pars), size = 10000, 
                         replace = TRUE, prob = pars$probability)

sample_indices <- sample(1:nrow(pars), size = 10000, 
                         replace = TRUE, prob = pars$probability)

head(sample_indices)

430  428 1010  383  343  385

pars_sample <- pars[sample_indices, c("mu", "sigma")]
head(pars_sample)

     mu sigma
1  20.0   2.8
2  19.0   2.8
3  17.5   6.7
4  19.0   2.5
5  21.5   2.2
6  20.0   2.5
7  20.0   2.8
8  20.5   1.6
9  19.0   2.5
10 17.0   4.0

平均気温の確率分布

hist(pars_sample$mu, 30)

平均気温の確率分布

quantile(pars_sample$mu, c(0.05, 0.95))

  5%  95% 
17.5 22.5

20日の気温は18度以上か？

pred_temp <- rnorm(10000, mean = , sd = )

20日の気温は18度以上か？

pred_temp <- rnorm(10000, mean = pars_sample$mu, sd = pars_sample$sigma)

20日の気温は18度以上か？

pred_temp <- rnorm(10000, mean = pars_sample$mu, sd = pars_sample$sigma)
hist(pred_temp, 30)

20日の気温は18度以上か？

pred_temp <- rnorm(10000, mean = pars_sample$mu, sd = pars_sample$sigma)
hist(pred_temp, 30)
sum(pred_temp >= 18) / length(pred_temp )

0.73

ゾンビのIQについて考えてみよう

RによるベイズデータAnalysisの基礎