Approximate inference for the multinomial logit model

Higher-order asymptotic theory is used to derive p-values that achieve superior accuracy compared to the p-values obtained from traditional tests for inference about parameters of the multinomial logit model. Simulations are provided to assess the finite sample behavior of the test statistics considered and to demonstrate the superiority of the higher-order method. Stata code that outputs these p-values is available to facilitate the implementation of these methods for the end-user.