The existing studies on drivers’ injury severity include numerous statistical models that assess potential factors affecting the level of injury. These models should address specific concerns tailored to different crash characteristics. For rear-end crashes, potential correlation in injury severity may present between the two drivers involved in the same crash. Moreover, there may exist unobserved heterogeneity considering parameter effects, which may vary across both crashes and individuals. To address these concerns, a random parameters bivariate ordered probit model has been developed to examine factors affecting injury sustained by two drivers involved in the same rear-end crash between passenger cars. Taking both the within-crash correlation and unobserved heterogeneity into consideration, the proposed model outperforms the two separate ordered probit models with fixed parameters. The value of the correlation parameter demonstrates that there indeed exists significant correlation between two drivers’ injuries. Driver age, gender, vehicle, airbag or seat belt use, traffic flow, etc., are found to affect injury severity for both the two drivers. Some differences can also be found between the two drivers, such as the effect of light condition, crash season, crash position, etc. The approach utilized provides a possible use for dealing with similar injury severity analysis in future work.