;   vector[nSubj] bSubj;   vector[nGroup] bGroup;   vector[nBev] bBev;   vector[nBev] bGxB[nGroup];   real<lower=0> sigma_s;   real<lower=0> sigma_g;   real<lower=0> sigma_b;   real<lower=0> sigma_gb; 