Modelling two-vehicle crash severity by generalized estimating equations

•Crash severity levels of both parties in two-vehicle accidents are modeled.•The generalized estimating equations method is used to estimate ordered probit models.•Comparisons with univariate and bivariate models in a detailed prediction performance matrix were conducted.•Key risk factors are identified and corresponding countermeasures are proposed.