fit.NGS {bda}R Documentation

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.

Author(s)

B. Wang bwang@southalabama.edu

References

AS 254, ...

Examples

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

[Package bda version 10.1.9 Index]