Abstract:
|
We present an example of a Bayesian hierarchical model, for an advanced undergraduate Bayesian Statistics class, with a focus on explaining the variability in coronavirus testing. This demonstration was applied to data from various Trump events throughout 2020, where testing and social distancing were not enforced. Using R and the JAGS software for Gibbs Sampling, we developed predictive, posterior distributions for the total number of people who would have tested positive for coronavirus at each (and every) event. Additionally, we developed posterior distributions for the coronavirus test sensitivity and specificity. The MCMC chains appeared to properly converge to the posterior distributions, and was further suggested so, by the Gelman-Rubin statistic.
|