How do I test the goodness of fit for multinomial logistic regression, having used the subpop option? So far I have tried using the mlogitgof and fitstat commands, but the test results were incomplete. The linktest command worked pretty well for the binary logistic regression models in the same analysis.