est_power: est_power

In RnaSeqSampleSize: RnaSeqSampleSize

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)

Example output

Loading required package: RnaSeqSampleSizeData
Loading required package: edgeR
Loading required package: limma
[1] 0.8

RnaSeqSampleSize documentation built on Nov. 8, 2020, 6:54 p.m.