Error message: Random variable[1] is nan, but must not be nan!

A tilde statement does not assign a value but only computes the probability of a preassigned value. If you want to draw random samples of alpha it needs to be declared in the parameters block.

parameters {
  ...
  real alpha[R,J-1];
}
model {
  ...
  for (r in 1:R) {
    ...
    alpha[r] ~ smvn(..);
    ...
  }
}
3 Likes