RnaSeqSampleSize (version 1.4.2)

est_power_curve: est_power_curve

Description

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

Usage

est_power_curve(n, w = 1, rho = 2, lambda0 = 5, phi0 = 1, alpha = 0.05, f = 0.05, ...)

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

...

other parameters for est_power function.

Value

A list including parameters, sample size and power.

Examples

## Not run: 
# result1<-est_power_curve(n=63, f=0.01, rho=2, lambda0=5, phi0=0.5)
# result2<-est_power_curve(n=63, f=0.05, rho=2, lambda0=5, phi0=0.5)
# plot_power_curve(list(result1,result2))
# ## End(Not run)