Thanks for the suggestion. This helped. I just had to made slight changes with the way I define intercept level random effect estimates: prior B: ~ normal(0, var=1e6); prior S2 ~ igamma(0.01, scale = 0.01); prior S2g ~ general(0, lower=0); random b0 ~ normal(0,var=S2g) subject=tier monitor=(b0) ; random b1~ normal(0,var=S2g) subject=tier monitor=(b1); random b2~ normal(0,var=S2g) subject=tier monitor=(b2); random b2~ normal(0,var=S2g) subject=tier monitor=(b3); Mu = (B0+b0) + (B1+b1)*P1 + (B2+b2)*P2 +(B3+b3)*P3 ;
... View more
