est_power {RnaSeqSampleSize} | R Documentation |
A function to estitamete the power for differential expression analysis of RNA-seq data.
est_power( n, w = 1, rho = 2, lambda0 = 5, phi0 = 1, alpha = 0.05, f, m = 20000, m1 = 200 )
n |
Numer of samples. |
w |
Ratio of normalization factors between two groups. |
rho |
minimum fold changes for prognostic genes between two groups. |
lambda0 |
Average read counts for prognostic genes. |
phi0 |
Dispersion for prognostic genes. |
alpha |
alpha level. |
f |
FDR level |
m |
Total number of genes for testing. |
m1 |
Expected number of prognostic genes. |
Estimate power
n<-63;rho<-2;lambda0<-5;phi0<-0.5;f<-0.01 est_power(n=n, rho=rho, lambda0=lambda0, phi0=phi0,f=f)