Robust adjusted profile score

Usage

mr_raps(b_exp, b_out, se_exp, se_out, parameters = default_parameters())

Arguments

b_exp: Vector of genetic effects on exposure.
b_out: Vector of genetic effects on outcome.
se_exp: Standard errors of genetic effects on exposure.
se_out: Standard errors of genetic effects on outcome.
parameters: A list of parameters. Specifically, over.dispersion and loss.function. over.dispersion is a logical concerning should the model consider overdispersion (systematic pleiotropy). And loss.function allows using either the squared error loss ("l2") or robust loss functions/scores ("huber" or "tukey"). The default is parameters=list(overdispersion = TRUE, loss.function = "tukey").

Value

List with the following elements:

b: MR estimate
se: Standard error
pval: p-value
nsnp: Number of SNPs

Details

This function calls the mr.raps package. Please refer to the documentation of that package for more detail.

References

Qingyuan Zhao, Jingshu Wang, Jack Bowden, Dylan S. Small. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Forthcoming.