plot: Plot the stratified power assessment criteria

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

Plot the stratified power assessment criteria

Description

Plot the stratified Power, FDR (False Discovery Rate), TP (True Positive), FP (False Positive) based on prevalence or abundance.

Usage

plotStrata(assess.out, figure.type=c('power','fdr','type1error','tp','fp','otu'),stratify.by=c('prevalence','abundance'),is.errorbar=TRUE,is.legend=TRUE) 
plotStrataAll(assess.out, figure.type=c('power','fdr','type1error','tp','fp','otu'),stratify.by=c('prevalence','abundance'),is.errorbar=TRUE,is.legend=TRUE)

Arguments

assess.out	An object output by assess function
figure.type	Figure type could be Power("power"), FDR ("fdr"), TypeI error ("type1error"), TP("tp"), FP("tp"), #OTU ("otu")
stratify.type	Taxon is stratified by either "prevalence" or "abundance"
is.errorbar	Standard error of all simulations is plot
is.legend	Legend of sample sizes is plot

Examples

## Not run: 
assess.out = assess(powmeta.out, alpha.type="fdr", alpha.level=0.1,stratify.type='prevalence')
plotStrataAll(assess.out,stratify.by='prevalence')
plotStrata(assess.out,figure.type='power',stratify.by='prevalence')
 
## End(Not run)

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