Addressing spatial heterogeneity of injury severity using Bayesian multilevel ordered probit model. (May 2020)