Support Normal and Gamma distributions

This uses the moment estimator of scale parameter as in R and obtains the same results. Because some of the tests (taken from R examples) use an offset, they will fail to converge without #58.