Peak caller based on the Zero Inflated Negative Binomial (ZINB) distribution Cusco &al 2016 Bioinformatics.
The algorithm works in two phases:
- Fit the distribution parameters on the entire vector by using multiple random initialization.
- Fit the distribution parameters on widows of size WS using the optimization on the entire vector as initialization.
To get a binary vector with a threshold at a p-value of
vectB = Peaks_Caller(vect)
To get the p-value vector
vectB = Peaks_Caller(vect,getPval=True)
Note: By default, the distribution parameters inside windows are computed in parallel,
using all the cores of the machine (Ncore=-1).