## ----style, eval=TRUE, echo=FALSE, results='asis'---------------------------------------
    BiocStyle::latex(use.unsrturl=FALSE)

## ----echo=TRUE,eval=FALSE---------------------------------------------------------------
#    library(PROPER)
#    sim.opts.Cheung = RNAseq.SimOptions.2grp(ngenes = 20000, p.DE=0.05,
#      lOD="cheung", lBaselineExpr="cheung")

## ----echo=TRUE,eval=FALSE---------------------------------------------------------------
#    sim.opts.Bottomly = RNAseq.SimOptions.2grp(ngenes = 20000, p.DE=0.05,
#      lOD="bottomly", lBaselineExpr="bottomly")

## ----echo=TRUE,eval=FALSE,result=FALSE--------------------------------------------------
#    simres = runSims(Nreps = c(3, 5, 7, 10), sim.opts=sim.opts.Cheung,
#      DEmethod="edgeR", nsims=20)

## ----echo=TRUE,eval=FALSE---------------------------------------------------------------
#    powers = comparePower(simres, alpha.type="fdr", alpha.nominal=0.1,
#      stratify.by="expr", delta=0.5)

## ----echo=TRUE,eval=FALSE---------------------------------------------------------------
#    summaryPower(powers)

## ----eval=FALSE,echo=TRUE---------------------------------------------------------------
#    plotPower(powers)

## ----eval=FALSE,echo=TRUE---------------------------------------------------------------
#    plotPowerTD(powers)

## ----eval=FALSE,echo=TRUE---------------------------------------------------------------
#    plotFDcost(powers)

## ----echo=TRUE,eval=FALSE---------------------------------------------------------------
#    plotAll(powers)

## ----echo=TRUE,eval=FALSE---------------------------------------------------------------
#    power.seqDepth(simres, powers)

## ----echo=TRUE,eval=FALSE---------------------------------------------------------------
#    powers = comparePower(simres, alpha.type="fdr", alpha.nominal=0.1,
#      strata = c(0, 10, 2^(1:7)*10, Inf), filter.by="expr",
#      strata.filtered=1, stratify.by="expr", delta=0.5)

## ----echo=TRUE,eval=FALSE---------------------------------------------------------------
#      powers = comparePower(simres, alpha.type="fdr", alpha.nominal=0.1,
#        stratify.by="dispersion", target.by="effectsize", delta=1)

## ----echo=TRUE,eval=FALSE---------------------------------------------------------------
#      powers = comparePower(simres, alpha.type="pval", alpha.nominal=0.001,
#        stratify.by="dispersion", target.by="effectsize", delta=1)

## ----echo=TRUE, result=TRUE-------------------------------------------------------------
sessionInfo()