Add support for other families of generalized linear models

[kschmaus]

R glm functions supports the following set: binomial(link = "logit"), gaussian(link = "identity"), Gamma(link = "inverse"), inverse.gaussian(link = "1/mu^2"), poisson(link = "log"), quasi(link = "identity", variance = "constant"), quasibinomial(link = "logit"), quasipoisson(link = "log"). We should try to support as many as possible.