first working implementation of hierarchical bayes variant of Pareto/NBD