Change possible scale range of alpha and delta

JAGS/cpt_hierarchical_prelec-98.txt

@@ -36,14 +36,14 @@ model
 
   ## retransfrom group level means and extend scales to a reasonable range (cf. Scheibehenne & Pachur, 2015)
 
-  mu.alpha <- 2*phi(mu.alpha.prob)
+  mu.alpha <- phi(mu.alpha.prob)
   mu.gamma <- 2*phi(mu.gamma.prob)
-  mu.delta <- 5*phi(mu.delta.prob)
+  mu.delta <- 10*phi(mu.delta.prob)
   mu.rho <- exp(mu.rho.log)
 
   # agent level
 
-  for (j in n_agents)
+  for (j in 1:n_agents)
   {
 
   ## transformed agent-level parameters are drawn from group-level normal distributions on the respective scale
@@ -62,7 +62,7 @@ model
 
     # define CPT model
 
-    for (i in start:stop)
+    for (i in 1:n_gambles)
 	  {
 
     # value function
@@ -76,8 +76,8 @@ model
 
 	  w.a.p2[i, j] <- exp(-delta[j]*(pow((-1*(log(a_p2_exp[i, j]))), gamma[j])))
 	  w.a.p1[i, j] <- 1-w.a.p2[i, j]
-	  w.b.p1[i, j] <- exp(-delta[j]*(pow((-1*(log(b_p1[i, j]))), gamma[j])))
-	  w.b.p2[j, i] <- 1-w.b.p1[j, i]
+	  w.b.p1[i, j] <- exp(-delta[j]*(pow((-1*(log(b_p1[i]))), gamma[j])))
+	  w.b.p2[i, j] <- 1-w.b.p1[i, j]
 
 	  Vf.a[i, j]  <- w.a.p1[i, j] * v.a.o1[i, j] + w.a.p2[i, j] * v.a.o2[i, j]
 	  Vf.b[i, j]  <- w.b.p1[i, j] * v.b.o1[i, j] + w.b.p2[i, j] * v.b.o2[i, j]