Thank you! It samples now…
I think that’s exactly connected to what I don’t understand about the coords / dimensions: why don’t I have to index adaptation (and the other participant specific parameters) in the code that follows with both indices like this
diff_adapt = pm.Deterministic("adaptation_diff", adaptation_diff * adaptation[participant_idx,exp_idx])