Posterior probabilities

You've used RJAGS output to explore and quantify the posterior trend & uncertainty \(b\). You can also use RJAGS output to assess specific hypotheses. For example: What's the posterior probability that, on average, weight increases by more than 1.1 kg for every 1 cm increase in height? That is, what's the posterior probability that \(b > 1.1\)?

You will approximate this probability by the proportion of \(b\) Markov chain values that exceed 1.1. The weight_chains data frame with the 100,000 iteration Markov chain output is in your workspace.

Construct a density plot of the \(b\) Markov chain values and use geom_vline() to superimpose a vertical line at 1.1.
Use table() to summarize the number of \(b\) Markov chain values that exceed 1.1.
Use mean() to calculate the proportion of \(b\) Markov chain values that exceed 1.1.

Exercise

Posterior probabilities

Instructions