estimate cpt over experienced probabilities

JAGS/cpt_trial_level.txt

@@ -27,7 +27,7 @@ model
 
 	# weighting function
 
-	w.a.p2[i] <- (delta * (pow(a_p2[i], gamma))) / (delta * (pow(a_p2[i], gamma)) + pow(a_p1[i], gamma))
+	w.a.p2[i] <- (delta * (pow(a_p2_exp[i], gamma))) / (delta * (pow(a_p2_exp[i], gamma)) + pow(a_p1_exp[i], gamma))
 	w.a.p1[i] <- 1-w.a.p2[i]
 	w.b.p1[i] <- (delta * (pow(b_p1[i], gamma))) / (delta * (pow(b_p1[i], gamma)) + pow(b_p2[i], gamma))
 	w.b.p2[i] <- 1-w.b.p1[i]

R/cpt.R

@@ -47,8 +47,8 @@ for(set in seq_len(nrow(params_sim))){
                  a_o2 = choices_grouped[[set]]$a_o2, # higher risky outcome (x)
                  b_o1 = choices_grouped[[set]]$b_o1, # safe outcome (x)
                  b_o2 = choices_grouped[[set]]$b_o2,
-                 a_p1 = choices_grouped[[set]]$a_p1,
-                 a_p2 = choices_grouped[[set]]$a_p2, # probability of higher risky outcome (x)
+                 a_p1_exp = choices_grouped[[set]]$a_p1_exp,
+                 a_p2_exp = choices_grouped[[set]]$a_p2_exp, # experienced probability of higher risky outcome (x)
                  b_p1 = choices_grouped[[set]]$b_p1,
                  b_p2 = choices_grouped[[set]]$b_p2, # probability of safe outcome
                  start = min(choices_grouped[[set]]$i),