Code to compare different models for branching processes with negative binomial offspring distribution