SISnRiskGroups function

SIS model with multiple risk groups (P 3.2).