fit.NGS {bda} | R Documentation |
Fit a finite normal mixture model for pre-binned NGS data
fit.NGS(x, iter.max=30,reflect=TRUE)
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. |
We fit a mixture model to take care the zero values using component-0.
B. Wang bwang@southalabama.edu
AS 254, ...
# To be updated. x = rexp(100,1) x = c(rep(0,20),x) out = fit.NGS(x) plot(out)