## 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)
Run the code above in your browser using DataLab