= Solution
\b[This is not a valid ordinary R glm family specification.] A <multinomial distribution> is natural for the three damage counts conditional on the row total, but base R's <generalized linear model> interface does not supply a family named multinomial. Furthermore, the written response is one cell count, not a three-category response or grouped multinomial vector. A properly specified <multinomial logistic regression> could model category probabilities by group, but it is not the model fit by this command. The independent-count <log-linear model> in part (a), conditional on the margins, is the appropriate supplied alternative.
Back to article page