implementation of Abe's adaptation of Pareto/NBD, which assumes lognormal distribution for lambda and theta
removed deprecated parameters
added implementation of Gamma/Gompertz/NBD model (by Bemmaor/Glady)