RnaSeqSampleSize (version 1.4.2)

est_power: est_power

Description

A function to estitamete the power for differential expression analysis of RNA-seq data.

Usage

est_power(n, w = 1, rho = 2, lambda0 = 5, phi0 = 1, alpha = 0.05, f, m = 20000, m1 = 200)

Arguments

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.

Value

Estimate power

Examples

n<-63;rho<-2;lambda0<-5;phi0<-0.5;f<-0.01
est_power(n=n, rho=rho, lambda0=lambda0, phi0=phi0,f=f)