bda (version 13.4.13)

deg.NGS: A test of DEG detection

Description

To test whether a gene differentially expressed in two experiments.

Usage

deg.NGS(x,y, gene,iter=1001)
  deg.test(normal.x, normal.y,counts.x, counts.y, gene,weights.x, weights.y)

Arguments

x,y

Two data files

counts.x, counts.y

Counts data files

normal.x, normal.y

Normalized data files

iter

maximum iteration to find the optimal fit.

gene

name of gene to be tested

weights.x, weights.y

weights

Details

We fit a mixture model to take care the zero values using component-0.

References

AS 254, ...

Examples

# NOT RUN {
 # To be updated.
 x = rexp(100,1)
 x = c(rep(0,20),x)
 out = fit.NGS(x)
 plot(out)
# }