2016/08/09 by Carlo Berzuini, Berzuini, Carlo, Hui Guo +5 · 1 citation
Biochemistry, Genetics and Molecular Biology · Mathematics · #FOS: Computer and information sciences #FOS: Mathematics #Genetic Associations and Epidemiology #Genetic and phenotypic traits in livestock #Methodology (stat.ME) #Statistical Methods and Inference #Statistics Theory (math.ST) #math.ST #stat.ME #stat.TH
paper · pdf · doi:10.48550/arxiv.1608.02990
21 pages, 6 figures and 2 tables
openalex publication_date 2016/08/09 · arxiv created 2017/01/31 · arxiv updated 2017/02/01 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
Our Bayesian approach to Mendelian Randomisation uses multiple instruments to assess the putative causal effect of an exposure on an outcome. The approach is robust to violations of the (untestable) Exclusion Restriction condition, and hence it does not require instruments to be independent of the outcome conditional on the exposure and on the confounders of the exposure-outcome relationship. The Bayesian approach offers a rigorous handling of the uncertainty (e.g. about the estimated instrument-exposure associations), freedom from asymptotic approximations of the null distribution and the possibility to elaborate the model in any direction of scientific relevance. We illustrate the last feature with the aid of a study of the metabolic mediators of the disease-inducing effects of obesity, where we elaborate the model to investigate whether the causal effect of interest interacts with a covariate. The proposed model contains a vector of unidentifiable parameters, β, whose jth element represents the pleiotropic (i.e., not mediated by the exposure) component of the association of instrument j with the outcome. We deal with the incomplete identifiability by assuming that the pleiotropic effect of some instruments is null, or nearly so, formally by imposing on β Carvalho's horseshoe shrinkage prior, in such a way that different components of β are subjected to different degrees of shrinking, adaptively and in accord with the compatibility of each individual instrument with the hypothesis of no pleiotropy. This prior requires a minimal input from the user. We present the results of a simulation study into the performance of the proposed method under different types of pleiotropy and sample sizes. Comparisons with the performance of the weighted median estimator are made. Choice of the prior and inference via Markov chain Monte Carlo are discussed.