fdr_msa: FDR-MSA correction

In nhejazi/methyvim: Targeted, Robust, and Model-free Differential Methylation Analysis

Description

Modified FDR Controlling Procedure for Multi-Stage Analyses (MJ van der Laan and C Tuglus, 2009, <doi:10.2202/1544-6115.1397>)

Usage

fdr_msa(pvals, total_obs)

Arguments

pvals	Numeric vector containing the p-values that result from any chosen statistical hypothesis testing procedure.
total_obs	Numeric indicating the total number of observations that would have been available for testing prior to the selection procedure employed in the multi-stage analysis performed.

Value

A numeric vector of corrected p-values, controlling the False Discovery Rate, using the method of Tuglus and van der Laan.

Examples

g <- 1e4
n <- 1e2
p <- abs(rnorm(n, mean = 1e-8, sd = 1e-2))
# treating the vector p as one of p-values, FDR-MSA may be applied
fdr_p <- fdr_msa(pvals = p, total_obs = g)

nhejazi/methyvim documentation built on April 30, 2020, 7:14 p.m.