modal.plot | R Documentation |
Modal plot to detect heterogeneity
modal.plot( b_exp = NULL, b_out = NULL, se_exp = NULL, se_out = NULL, data = NULL, k = 1.5, weight.option = c("MLE", "shrinkage"), beta.range = NULL )
data |
Alternatively, dataset can be passed by the argument |
k |
Locality of the robust likelihood (smaller |
beta.range |
range of beta in the plot |
data(lipid.cad) data <- subset(lipid.cad, lipid == "hdl" & restrict & gwas.selection == "teslovich_2010" & gwas.outcome == "cardiogramplusc4d_1000genome" & pval.selection < 1e-5) modal.plot(data$beta.exposure, data$beta.outcome, data$se.exposure, data$se.outcome, k = 1) data <- subset(lipid.cad, lipid == "ldl" & restrict & gwas.selection == "teslovich_2010" & gwas.outcome == "cardiogramplusc4d_1000genome" & pval.selection < 1e-5) modal.plot(data$beta.exposure, data$beta.outcome, data$se.exposure, data$se.outcome)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.