theta_plot function

Plot theta values of negative binomial models versus non-zero count for count data