est_power_curve: est_power_curve

In slzhao/RnaSeqSampleSize: RnaSeqSampleSize

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 (Treatment/Control).
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))

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