assess: Perform power assessment

In lichen-lab/powmic: powmic: an R package for power assessment in microbiome case-control

Perform power assessment

Description

Assess the simulation results from powmic by providing overall and stratified FDR (False Discovery Rate), TP (True Positive), FP (False Positive), TPR (True Positive Rate), FP (False Positive Rate)

Usage

assess(powmic.out, alpha.type=c("fdr","pval"), alpha.level=0.1,stratify.type=c("prevalence","abundance"))

Arguments

powmic.out	An powmic object
alpha.type	Power assessment using FDR ("fdr") or pvalue ("pval")
alpha.level	Significant level. Default is 0.1
stratify.type	OTU is stratified by either "prevalence" or "abundance"

See Also

summaryAssess, plotStrata, plotStrataAll

Examples

 powmic.out=powmic(n1s=c(20,40,60), n2s=c(20,40,60), nsims=10, 
                 sim.params=sim.params,distrib=distrib,DAmethod='edgeR')
 assess.out = assess(powmic.out, alpha.type='fdr',alpha.level=0.1,stratify.type='prevalence')

lichen-lab/powmic documentation built on April 7, 2023, 4:40 p.m.