[i] = ordered_logistic_rng(a +                                        (bGroup[gIndex[i]] * sigma_g) +                                        (bBev[bIndex[i]] * sigma_b) +                                        (bSubj[sIndex[i]] * sigma_s) +                                        (bGxB[gIndex[i]][bIndex[i]] * sigma_gb),                                      cutpoints);   } 