Is there a way to get the result of a vectorized probability functions * not* as the sum of all evaluations? I want to avoid a

`for`

loop over cumulative normal distributions and write something like this:`lt += bernoulli_lpmf(choice[:Tr] | normal_cdf(0, to_vector(delta[:Tr]), weber * to_vector(Total[:Tr])));`

instead of:

```
for (t in Tr) {
lt += bernoulli_lpmf(choice[t] | normal_cdf(0,delta[t], weber *(Total[t])));
}
```

Thanks!!