Description Usage Arguments Details Value References Examples
View source: R/fdr-sampsize-v1.0e.R
Use the formula of Li et al (2013) to compute power for comparing RNA-seq expression across two groups assuming the Poisson distribution.
1 | power.li (n, alpha, rho, mu0, w = 1, type = "w")
|
n |
per-group sample size |
alpha |
p-value threshold |
rho |
fold-change, usual null hypothesis is that rho=1 |
mu0 |
average count in control group |
w |
ratio of total number of |
type |
type of test: "w" for Wald, "s" for score, "lw" for log-transformed Wald, "ls" for log-transformed score. |
This function computes the power for each of a series of two-sided tests defined by the input parameters. The power is based on the sample size formulas in equations 10-13 of Li et al (2013). Also, note that the null.effect is set to 1 in the examples because the usual null hypothesis is that the fold-change = 1.
vector of power estimates for two-sided tests
C-I Li, P-F Su, Y Guo, and Y Shyr (2013). Sample size calculation for differential expression analysis of RNA-seq data under Poisson distribution. Int J Comput Biol Drug Des 6(4). doi:10.1504/IJCBDD.2013.056830
1 2 3 4 5 6 7 8 9 |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.