mr_radialegger_rjags: Bayesian radial MR-Egger model with a choice of prior...
In mrbayes: Bayesian Summary Data Models for Mendelian Randomization Studies

Description Usage Arguments Value References Examples

Bayesian radial MR-Egger model with a choice of prior distributions fitted using JAGS.

mr_radialegger_rjags(
  object,
  prior = "default",
  betaprior = "",
  sigmaprior = "",
  n.chains = 3,
  n.burn = 1000,
  n.iter = 5000,
  seed = NULL,
  rho = 0.5,
  ...
)

`object`	A data object of class `mr_format`.
`prior`	A character string for selecting the prior distributions; `"default"` selects a non-informative set of priors; `"weak"` selects weakly informative priors; `"pseudo"` selects a pseudo-horseshoe prior on the causal effect; `"joint"` selects a joint prior on the intercept and slope.
`betaprior`	A character string in JAGS syntax to allow a user defined prior for the causal effect.
`sigmaprior`	A character string in JAGS syntax to allow a user defined prior for the residual standard deviation.
`n.chains`	Numeric indicating the number of chains used in the MCMC estimation, the default is `3` chains.
`n.burn`	Numeric indicating the burn-in period of the Bayesian MCMC estimation. The default is `1000` samples.
`n.iter`	Numeric indicating the number of iterations in the Bayesian MCMC estimation. The default is `5000` iterations.
`seed`	Numeric indicating the random number seed. The default is the rjags default.
`rho`	Numeric indicating the correlation coefficient input into the joint prior distribution. The default is `0.5`.
`...`	Additional arguments passed through to `rjags::jags.model()`.

An object of class radialeggerjags containing the following components:

AvgPleio: The mean of the simulated pleiotropic effect
CausalEffect: The mean of the simulated causal effect
StandardError: Standard deviation of the simulated causal effect
sigma: The mean of the simaulted residual standard deviation
CredibleInterval: The credible interval for the causal effect, which includes the lower (2.5%), median (50%) and upper intervals (97.5%)
samples: Output of the Bayesian MCMC samples
Prior: The specified priors

Bowden, J., et al., Improving the visualization, interpretation and analysis of two-sample summary data Mendelian randomization via the Radial plot and Radial regression. International Journal of Epidemiology, 2018. 47(4): p. 1264-1278. doi: 10.1093/ije/dyy101.

if (requireNamespace("rjags", quietly = TRUE)) {
fit <- mr_radialegger_rjags(bmi_insulin)
summary(fit)
plot(fit$samples)
# 90% credible interval
fitdf <- do.call(rbind.data.frame, fit$samples)
cri90 <- quantile(fitdf$Estimate, probs = c(0.05,0.95))
print(cri90)
}