A new flexible skewed bimodal distribution with multivariate extensions: Theory and application to traffic crash injury severity analysis