With , the fitted logistic model is
Risk category 1 is the reference level, so its effect is included in the intercept; adding a separate coefficient for both levels would make the design matrix linearly dependent.
Holding risk category fixed, increasing pollution by one unit multiplies the disease odds by
Thus the fitted odds increase by about twenty percent per unit of pollution.
There are observations. The null model estimates one intercept, so its residual degrees of freedom are . The fitted model estimates three coefficients, so its residual degrees of freedom are .
Fit the intercept-only model with
glm(dis ~ 1, family = binomial, data = disease)
Its deviance is . Since it has one fitted coefficient, its AIC on the same convention is
The first likelihood-ratio test compares
with the alternative that at least one is nonzero. The deviance drop is compared with , giving ; there is strong evidence that the covariates improve the model.
The second test compares the additive model with
against a model with a pollution-by-risk interaction. Its likelihood-ratio p-value is , so there is no evidence for an interaction.

Articles by others on the same topic (0)

There are currently no matching articles.