add loglike to family, use deviance based on loglike in gamma (now correct deviance), add dof correction for aic for known/estimated scale (unclear for gamma)
add estimated scale count in aic of regression.py