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.