Published online by Cambridge University Press: 11 September 2020
Bonus-malus systems (BMSs) are widely used in actuarial sciences. These systems are applied by insurance companies to distinguish the policyholders by their risks. The most known application of BMS is in automobile third-party liability insurance. In BMS, there are several classes, and the premium of a policyholder depends on the class he/she is assigned to. The classification of policyholders over the periods of the insurance depends on the transition rules. In general, optimization of these systems involves the calculation of an appropriate premium scale considering the number of classes and transition rules as external parameters. Usually, the stationary distribution is used in the optimization process. In this article, we present a mixed integer linear programming (MILP) formulation for determining the premium scale and the transition rules. We present two versions of the model, one with the calculation of stationary probabilities and another with the consideration of multiple periods of the insurance. Furthermore, numerical examples will also be given to demonstrate that the MILP technique is suitable for handling existing BMSs.