mr_radialegger_stan {mrbayes}

Bayesian inverse variance weighted model with a choice of prior distributions fitted using RStan.

Description

Bayesian inverse variance weighted model with a choice of prior distributions fitted using RStan

Usage

mr_radialegger_stan(
  data,
  prior = 1,
  n.chains = 3,
  n.burn = 1000,
  n.iter = 5000,
  rho = 0.5,
  seed = 12345,
  ...
)

Arguments

data

A data of class mr_format.

prior

An integer for selecting the prior distributions;

• 1 selects a non-informative set of priors;

• 2 selects weakly informative priors;

• 3 selects a pseudo-horseshoe prior on the causal effect;

• 4 selects joint prior of the intercept and causal effect estimate.

n.chains

Numeric indicating the number of chains used in the HMC estimation in rstan, the default is 3 chains.

n.burn

Numeric indicating the burn-in period of the Bayesian HMC estimation. The default is 1000 samples.

n.iter

Numeric indicating the number of iterations in the Bayesian HMC estimation. The default is 5000 iterations.

rho

Numeric indicating the correlation coefficient input into the joint prior distribution. The default is 0.5.

seed

Numeric indicating the random number seed. The default is 12345.

...

Additional arguments passed through to rstan::sampling().

Value

An object of class stanfit.

References

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.

Stan Development Team (2020). "RStan: the R interface to Stan." R package version 2.19.3, https://mc-stan.org/.

Examples

if (requireNamespace("rstan", quietly = TRUE)) {
# Note we recommend setting n.burn and n.iter to larger values
radegger_fit <- mr_radialegger_stan(bmi_insulin, n.burn = 500, n.iter = 1000)
print(radegger_fit)
}