A Bayesian Survival Analysis of a Historical Dataset: How Long Do Popes Live?

University courses in statistical modeling often place great emphasis on methodological theory, illustrating it only briefly by means of limited and repeatedly used standard examples. Unfortunately, this approach often fails to actively engage and motivate students in their learning process. The teaching of statistical topics such as Bayesian survival analysis can be enhanced by focusing on innovative applications. Here, we discuss the visualization and modeling of a dataset of historical events comprising the post-election survival times of popes. Inference, prediction, and model checking are performed in the Bayesian framework, with comparisons being made with the frequentist approach. Further opportunities for similar statistical investigations are outlined. Supplementary materials for this article are available online.