Diff of /man/plot.bioTMLE.Rd [000000] .. [efa494]

Switch to side-by-side view

--- a
+++ b/man/plot.bioTMLE.Rd
@@ -0,0 +1,51 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/plots.R
+\name{plot.bioTMLE}
+\alias{plot.bioTMLE}
+\title{Plot p-values from moderated statistical tests for class biotmle}
+\usage{
+\method{plot}{bioTMLE}(x, ..., type = "pvals_adj")
+}
+\arguments{
+\item{x}{object of class \code{biotmle} as produced by an appropriate call
+to \code{biomarkertmle}}
+
+\item{...}{additional arguments passed \code{plot} as necessary}
+
+\item{type}{character describing whether to provide a plot of unadjusted or
+adjusted p-values (adjustment performed via Benjamini-Hochberg)}
+}
+\value{
+object of class \code{ggplot} containing a histogram of the raw or
+ Benjamini-Hochberg corrected p-values (depending on user input).
+}
+\description{
+Histogram of raw or FDR-adjusted p-values from the moderated t-test.
+}
+\examples{
+\dontrun{
+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)
+
+plot(x = limmaTMLEout, type = "pvals_adj")
+}
+}