bda (version 10.1.9)

fit.NGS: Fit a Normal Mixture Model to binned NGS data

Description

Fit a finite normal mixture model for pre-binned NGS data

Usage

fit.NGS(x, iter.max=30,reflect=TRUE)

Arguments

x

Next-Generation Sequencing data or similar. Could have a lot of zeros but negative values are not allowed

iter.max

maximum iteration to find the optimal fit.

reflect

Use reflect and replicate method to adjust for boundary effect.

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