                                   (bGroup[gIndex[i]] * sigma_g) +                                   (bBev[bIndex[i]] * sigma_b) +                                   (bSubj[sIndex[i]] * sigma_s) +                                   (bGxB[gIndex[i]][bIndex[i]] * sigma_gb),                                 cutpoints);   } }  generated quantities {   real y_draw[N];   real 