Svara på frågan: Ska jag ha en strandfest?

Grunderna i bayesiansk dataanalys i R

Rasmus Bååth

Data Scientist

Frågorna

Vad är troligtvis medelvattentemperaturen den 20 juli?
Hur stor är sannolikheten att vattentemperaturen är 18 grader eller mer nästa 20 juli?

Posteriorifördelningen

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

Sannolikhetsfördelningen över medeltemperaturen

hist(pars_sample$mu, 30)

Sannolikhetsfördelningen över medeltemperaturen

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

  5%  95% 
17.5 22.5

Är temperaturen 18 grader eller mer den 20:e?

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

Är temperaturen 18 grader eller mer den 20:e?

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

Är temperaturen 18 grader eller mer den 20:e?

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

Är temperaturen 18 grader eller mer den 20:e?

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

Vad sägs om zombiers IQ?

Grunderna i bayesiansk dataanalys i R