Fit GLMs with appropriate distributions and link functions
You have a response variable and several predictors and must fit one or more GLMs in R using glm().
Choose the response family (Poisson for counts, binomial for binary, Gamma for positive skewed amounts) and the link function (log for log-linear effects, logit for log-odds, identity for additive). Fit with glm(y ~ x1 + x2, family = ...). Iterate by adding interactions and basis splines suggested by EDA. Watch convergence warnings — they often signal separation or numerical issues.
glm(formula, family = family(link = link), data = df).