bda (version 10.1.9)

fit.lognormal: Fit a log-normal distribution to (binned) NGS data

Description

Fit a log-normal distribution to (binned) NGS data

Usage

fit.lognormal(x, x.limits)

Arguments

x

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

x.limits

a low threshold of a non-zero measure which depends on a specific experiment.

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