power.seqDepth {PROPER} | R Documentation

## Estimate the marginal power under different sequencing depth

### Description

The function helps study the power under different sequencing depth and sample sizes. It estimates the marginal powers based on existing simulation results, given new sequencing depth.

### Usage

    power.seqDepth(simResult, powerOutput, depth.factor = c(0.2, 0.5, 1, 2, 5, 10))

### Arguments

`simResult` | Result object from "runSims" function.
`powerOutput` | Result object from "comparePower" function.
`depth.factor` | A vector of numbers specifying the *relative* sequencing depth, comparing to the depth used in the simulation. 1 means using the same number of total reads as the simulation.

### Details

The powers under different sequencing depth and sample sizes provides important guidence in experimental design. Under the same total number of sequence reads, investigator can choose to use more replicates and shallower coverage for each, or less replicates and deeper coverage.

This function provides estimated marginal power holding all experimental variables fixed (biological variation, effect sizes, sample sizes, etc.) except the sequencing depth. Changing sequencing depth will only alter the marginal distribution of average counts. Since the stratified power (by average counts) won't change, those numbers are used in estimating the powers under different depth. This approaches allows skipping new simulations, which saves computation.

### Value

A matrix for marginal powers. Each row is for a sequencing depth, each columns is for a sample size.

### Author(s)

Hao Wu <hao.wu@emory.edu>

### See Also

comparePower, summary.power

### Examples

    ## Not run:
    simOptions = RNAseq.SimOptions.2grp()
    simRes = runSims(Nreps=c(3,5,7), sim.opts=simOptions, nsims=5,
                     DEmethod="edgeR")
    powers = comparePower(simRes)
    power.seqDepth(simRes, powers)

    ## End(Not run)

[Package PROPER version 1.18.0 Index]