Robust adjusted profile score

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.