Bayesian Mendelian randomization testing of interval causal null hypotheses: ternary decision rules and loss function calibration

Abstract

We enhance the Bayesian Mendelian Randomization (MR) analysis framework of Berzuini et al. [13] by incorporating a novel Bayesian method for testing interval null causal hypotheses, in the interest of a healthier approach to causal discovery. A number of ideas are combined in our proposal. First, we replace the usual point null hypothesis for the causal effect with a region of practical equivalence (ROPE), in the spirit of Kruschke [1]. Second, we allow the hypothesis test decision to be taken on the basis of the Bayesian posterior odds of the causal effect falling within the ROPE, calculated via a new "Markov chain Monte Carlo cum Importance sampling" procedure. Third, we allow the test decision rule to accommodate an uncertain outcome, for those situations where the posterior odds is neither large nor small enough. Finally, we present an approach to calibration of the proposed method via loss function. We illustrate the method with the aid of a study of the causal effect of obesity on risk of juvenile myocardial infarction based on a unique prospective dataset.