Trả lời câu hỏi: Có nên tổ chức tiệc bãi biển không?

Những nền tảng của Phân tích dữ liệu Bayesian trong R

Rasmus Bååth

Data Scientist

Các câu hỏi

Nhiệt độ nước trung bình vào ngày 20/7 thường là bao nhiêu?
Xác suất để nhiệt độ nước đạt từ 18 trở lên vào ngày 20 tới là bao nhiêu?

Phân phối hậu nghiệm

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

Phân phối xác suất của nhiệt độ trung bình

hist(pars_sample$mu, 30)

Phân phối xác suất của nhiệt độ trung bình

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

  5%  95% 
17.5 22.5

Nhiệt độ ngày 20 có đạt từ 18 trở lên không?

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

Nhiệt độ ngày 20 có đạt từ 18 trở lên không?

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

Nhiệt độ ngày 20 có đạt từ 18 trở lên không?

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

Nhiệt độ ngày 20 có đạt từ 18 trở lên không?

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

Còn IQ của zombie thì sao?

Những nền tảng của Phân tích dữ liệu Bayesian trong R