skutil.linear_model will provide the basis for Poisson and Tweedy regressions and gradient computations.