I agree to what you say, but also with diff_adapt = diff * adaptation[exp_idx] I get the error about non-matching shapes.
(of course also changing [exp_idx] for intercept and level)
I agree to what you say, but also with diff_adapt = diff * adaptation[exp_idx] I get the error about non-matching shapes.
(of course also changing [exp_idx] for intercept and level)