makemyprior: Intuitive Construction of Joint Priors for Variance Parameters in R

Priors allow us to robustify inference and to incorporate expert knowledge in Bayesian hierarchical models. This is particularly important when there are random effects that are hard to identify based on observed data. The challenge lies in understanding and controlling the joint influence of the priors for the variance parameters, and makemyprior is an R package that guides the formulation of joint prior distributions for variance parameters. A joint prior distribution is constructed based on a hierarchical decomposition of the total variance in the model along a tree, and takes the entire model structure into account. Users input their prior beliefs or express ignorance at each level of the tree. Prior beliefs can be general ideas about reasonable ranges of variance values and need not be detailed expert knowledge. The constructed priors lead to robust inference and guarantee proper posteriors. A graphical user interface facilitates construction and assessment of different choices of priors through visualization of the tree and joint prior. The package aims to expand the toolbox of applied researchers and make priors an active component in their Bayesian workflow.

[1]  Kevin Kröninger,et al.  BAT - The Bayesian Analysis Toolkit , 2008, Comput. Phys. Commun..

[2]  Haavard Rue,et al.  Bayesian Computing with INLA: A Review , 2016, 1604.00860.

[3]  M. Degroot,et al.  Information about Hyperparameters in Hierarchical Models , 1981 .

[4]  S. D. Winter,et al.  The Importance of Prior Sensitivity Analysis in Bayesian Statistics: Demonstrations Using an Interactive Shiny App , 2020, Frontiers in Psychology.

[5]  Joseph Hilbe,et al.  Data Analysis Using Regression and Multilevel/Hierarchical Models , 2009 .

[6]  Andrew Gelman,et al.  The Prior Can Often Only Be Understood in the Context of the Likelihood , 2017, Entropy.

[7]  Robust modeling of additive and nonadditive variation with intuitive inclusion of expert knowledge , 2021, Genetics.

[8]  Andrew Gelman,et al.  The No-U-turn sampler: adaptively setting path lengths in Hamiltonian Monte Carlo , 2011, J. Mach. Learn. Res..

[9]  H. Rue,et al.  Bayesian bivariate meta‐analysis of diagnostic test studies with interpretable priors , 2015, Statistics in medicine.

[10]  Sara Martino,et al.  Animal Models and Integrated Nested Laplace Approximations , 2013, G3: Genes, Genomes, Genetics.

[11]  Martyn Plummer,et al.  JAGS Version 3.3.0 user manual , 2012 .

[12]  Leonhard Held,et al.  Gaussian Markov Random Fields: Theory and Applications , 2005 .

[13]  J. Wakefield,et al.  Bayesian inference for generalized linear mixed models. , 2010, Biostatistics.

[14]  Aki Vehtari,et al.  Bayesian Workflow. , 2020, 2011.01808.

[15]  I. Steinsland,et al.  Flexible modelling of spatial variation in agricultural field trials with the R package INLA , 2019, Theoretical and Applied Genetics.

[16]  A. Bentley,et al.  A Two‐Part Strategy for Using Genomic Selection to Develop Inbred Lines , 2017 .

[17]  J. Besag,et al.  Bayesian image restoration, with two applications in spatial statistics , 1991 .

[18]  Sarah Depaoli,et al.  Where Do Priors Come From? Applying Guidelines to Construct Informative Priors in Small Sample Research , 2017 .

[19]  Haavard Rue,et al.  Penalised Complexity Priors for Stationary Autoregressive Processes , 2016, 1608.08941.

[20]  Jon Wakefield,et al.  Disease mapping and spatial regression with count data. , 2007, Biostatistics.

[21]  Jiqiang Guo,et al.  Stan: A Probabilistic Programming Language. , 2017, Journal of statistical software.

[22]  Gregor Gorjanc,et al.  AlphaSim: Software for Breeding Program Simulation , 2016, The plant genome.

[23]  Andrea Riebler,et al.  Robust modeling of additive and nonadditive variation with intuitive inclusion of expert knowledge , 2020, bioRxiv.

[24]  S. D. Winter,et al.  Dangers of the Defaults: A Tutorial on the Impact of Default Priors When Using Bayesian SEM With Small Samples , 2020, Frontiers in Psychology.

[25]  Andrea Riebler,et al.  An intuitive Bayesian spatial model for disease mapping that accounts for scaling , 2016, Statistical methods in medical research.

[26]  Thiago G. Martins,et al.  Penalising Model Component Complexity: A Principled, Practical Approach to Constructing Priors , 2014, 1403.4630.

[27]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .

[28]  Kevin Fiedler,et al.  Likelihood Bayesian And Mcmc Methods In Quantitative Genetics , 2016 .

[29]  Haavard Rue,et al.  Intuitive Joint Priors for Variance Parameters , 2019, Bayesian Analysis.

[30]  David R. Jones,et al.  How vague is vague? A simulation study of the impact of the use of vague prior distributions in MCMC using WinBUGS , 2005, Statistics in medicine.

[31]  J. Gabry,et al.  Bayesian Applied Regression Modeling via Stan , 2016 .

[32]  Howard D. Bondell,et al.  Bayesian Regression Using a Prior on the Model Fit: The R2-D2 Shrinkage Prior , 2016 .

[33]  Anders Nielsen,et al.  TMB: Automatic Differentiation and Laplace Approximation , 2015, 1509.00660.