Hi!
In my programme’s parameters block I have this alpha_0 defined as:
row_vector[2] alpha_0[num_prod];
Later on, in the model block, I assign a prior to those parameters, vectorized over num_prod:
for (k in 1:2) {
alpha_0[k] ~ normal(0, 1);
}
This has been working fine until I faced an unusual case of num_prod = 1. In this case, I believe alpha_0 should be a row vector of length two, so the above approach should still work, but it does not. How can I solve this issue? Thanks in advance!