est_power_curve: est_power_curve

In RnaSeqSampleSize: RnaSeqSampleSize

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

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))

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