normal_rng will return a scalar value when the arguments (mean, sd) are scalars.
Solutions include:
- using a loop to call
normal_rng(0, 2)for each element - passing vectors as arguments to
normal_rng
For more info check out this thread: Applying Vectorization in generated quantities section