powered by
The function initially intends to rework on GSMR outputs, but it would be appropriate for general use.
gsmr(data, X, Y, alpha = 0.05, other_plots = FALSE)
The result and plots.
Data to be used.
Exposure.
Outcome.
type I error rate for confidence intervals.
To add funnel and forest plots.
library(cowplot) library(ggplot2) library(gap) r <- gsmr(mr, "LIF.R", c("CAD","FEV1"))
Run the code above in your browser using DataLab