this looks wrong both, reg and unreg, have the extra parameter
L1 counts shape parameter in df_model and df_resid, NegativeBinomial does't
?