ตอบคำถาม: ควรจัดปาร์ตี้ริมชายหาดไหม?

พื้นฐานการวิเคราะห์ข้อมูลแบบ Bayesian ใน R

Rasmus Bååth

Data Scientist

คำถามที่ต้องการคำตอบ

อุณหภูมิน้ำเฉลี่ยในวันที่ 20 กรกฎาคมน่าจะเป็นเท่าไหร่?
ความน่าจะเป็นที่อุณหภูมิน้ำจะอยู่ที่ 18 องศาขึ้นไปในวันที่ 20 ครั้งหน้า คือเท่าไหร่?

การแจกแจงภายหลัง (Posterior Distribution)

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 ของซอมบี้ล่ะ?

พื้นฐานการวิเคราะห์ข้อมูลแบบ Bayesian ใน R