plotPowerHist: Plot the histogram of power

In haowulab/PROPER: PROspective Power Evaluation for RNAseq

Description

Using histogram-alike graph to visualize the distribution of all genes and DE genes in all user specified strata.

Usage

plotPowerHist(powerOutput, simResult, main = "Histogram of power", return = FALSE)

Arguments

powerOutput	Result object from "comparePower" function.
simResult	Result object from "runSims" function.
main	Figure caption.
return	Return a matrix for average number of genes and true DE genes in each strata.

Value

A matrix of two rows. First row is for number of genes, and the second row is for number of DE genes. Columns are strata.

Author(s)

Hao Wu <[e-mail]>

Examples

## Not run: 
simOptions = RNAseq.SimOptions.2grp()
## run a few simulations
simRes = runSims(Nreps=c(3,5,7), sim.opts=simOptions, nsims=5,
                 DEmethod="edgeR")

## using FDR 0.1 to call DE, then look at power curves and summary
powers = comparePower(simRes)
plotPowerHist(powers, simRes)

## End(Not run)

haowulab/PROPER documentation built on Nov. 30, 2020, 2:22 a.m.