|
@@ -129,13 +129,13 @@ model {
|
|
|
//alphas_fp ~ gamma(2, 1);
|
|
|
for (i in 1:n_classes) {
|
|
|
//conf_sd[i,:] ~ normal(0, 1);
|
|
|
- omega[i,:] ~ pareto(1, 1.5);
|
|
|
+ omega[i,:] ~ pareto(1, 2);
|
|
|
//lambda_fp[:,i] ~ gamma(alphas_fp[i], alphas_fp[i]/mus_fp[i]);
|
|
|
for (j in 1:n_classes) {
|
|
|
mus[1,i,j] ~ exponential(1); // urban
|
|
|
mus[2,i,j] ~ exponential(1); // rural
|
|
|
- alphas[1,i,j] ~ lognormal(0,2); // urban
|
|
|
- alphas[2,i,j] ~ lognormal(0,2); // rural
|
|
|
+ alphas[1,i,j] ~ lognormal(0,1); // urban
|
|
|
+ alphas[2,i,j] ~ lognormal(0,1); // rural
|
|
|
for (c in 1:n_groups) {
|
|
|
lambda[c,i,j] ~ gamma(alphas[recording_rural[c]>0?2:1,i,j], alphas[recording_rural[c]>0?2:1,i,j]/(mus[recording_rural[c]>0?2:1,i,j]));
|
|
|
}
|