code issue in "bart_par4"

[vincent5290]

In BART_par4 model, the formula is l - vk, the optimal number of pumps on trial k (paper, titled "Development ofa novel computational model for the Balloon Analogue Risk Task: The exponential-weight mean–variance model"). However, in stanfile of "BART_par4", line 73 d[j, k, l] ~ bernoulli_logit(tau[j] * (omega - l)); is omega - l. The code seems contrast to the formula.