Volcano plot of the log-changes in the target causal paramter against the log raw p-values from the moderated t-test.

volcano_ic(biotmle, ate_bound = 1, pval_bound = 0.2)

biotmle | object of class |
---|---|

ate_bound | A |

pval_bound | A |

object of class `ggplot`

containing a standard volcano plot of
the log-fold change in the causal target parameter against the raw log
p-value computed from the moderated tests in `modtest_ic`

.

if (FALSE) { library(dplyr) library(biotmleData) library(SuperLearner) library(SummarizedExperiment) data(illuminaData) colData(illuminaData) <- colData(illuminaData) %>% data.frame() %>% mutate(age = as.numeric(age > median(age))) %>% DataFrame() benz_idx <- which(names(colData(illuminaData)) %in% "benzene") biomarkerTMLEout <- biomarkertmle( se = illuminaData, varInt = benz_idx, bppar_type = BiocParallel::SerialParam(), g_lib = c("SL.mean", "SL.glm"), Q_lib = c("SL.mean", "SL.glm") ) limmaTMLEout <- modtest_ic(biotmle = biomarkerTMLEout) volcano_ic(biotmle = limmaTMLEout) }