est_power: est_power

In slzhao/RnaSeqSampleSize: RnaSeqSampleSize

est_power

Description

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

Usage

est_power(
  n,
  w = 1,
  k = 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.
k	Ratio of sample size between two groups (Treatment/Control).
rho	minimum fold changes for prognostic genes between two groups (Treatment/Control).
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)

slzhao/RnaSeqSampleSize documentation built on March 24, 2022, 2:21 a.m.