bda (version 10.1.9)

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)

Arguments

x,y

Two data files

iter

maximum iteration to find the optimal fit.

gene

name of gene to be tested

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