We know the generalized linear models (GLMs) are a broad class of models. When fitting GLMs in R, we need to specify which family function to use from a bunch of options like gaussian, poisson, binomial, quasi, etc. I would like to have your advice regarding how to determine the optional family function used for GLM fitting in R. Thanks!