data { int N; real values[N]; int cluster_num; } transformed data { real min_val = min(values); real max_val = max(values); real var_val = variance(values); } parameters { positive_ordered[cluster_num] mu; positive_ordered[cluster_num] k; simplex[cluster_num] ratio; } transformed parameters { real<lower=min_val/var_val> rate[cluster_num]; for (j in 1:cluster_num) { rate[j] = k[j] / mu[j]; } } mode