fdr_msa: FDR-MSA correction

Description Usage Arguments Value Examples

View source: R/fdr_msa.R

Description

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

Usage

1
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

1
2
3
4
5
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)

methyvim documentation built on Nov. 8, 2020, 11:11 p.m.