Gene Environment Interaction

Material related to the manuscript "Bayesian inference of gene-environment interaction from incomplete data: what happens when information on environment is disjoint from data on gene and disease," by Paul Gustafson and Igor Burstyn, to appear in Statistics in Medicine.

Here are the supplementary figures not appearing in the published manuscript.

Here is the R function to evaluate the limiting posterior distribution of the gene-environment interaction coefficient for given hyperparameters and true parameter values.

Here are R functions to obtain posterior distributions for retrospective analysis, using full (X,G|Y) data or reduced (G|Y) data plus a prior on exposure prevalence. Some test code illustrates their use.