We assumed that the distribution of the number of secondary cases generated by a single primary case follows a negative binomial distribution with the basic reproduction number R0, i.e., the average number of secondary cases generated by a single primary case, and the dispersion parameter k. The probability of extinction π defined by the first generating moment [15] is then modeled as:(3) π=1(1+R0k(1−π))k.