Generalized linear model with log link and Poisson family for non-negative integer outcomes (counts).
Poisson regression models log(E[Y]) as a linear function of predictors. The exponentiated β is a rate ratio — exp(β) = factor by which the expected count multiplies per unit change in the predictor. For rates with varying exposure (events per person-year), add log(exposure) as an offset via the GLM-general handler.
The Poisson distribution assumes mean = variance. Real-world counts are usually overdispersed (variance > mean) because of unmeasured heterogeneity — the residual deviance / df ratio diagnoses this. > 1.5 ⇒ try negative binomial; > 2 ⇒ definitely use negbin.
Integer-valued response + numeric predictors + intercept toggle.
β + SE + Wald z + p per coefficient + rate ratios = exp(β) + CI + AIC + McFadden R² + overdispersion = residual deviance / df.
Rate ratio is the canonical effect-size. RR = 1.5 means a 50% increase in the expected count per unit of the predictor.