Closed seabbs closed 3 months ago
This doesn't look like it works when bp_sd is length 0 so won't be feasible for us
Could we not use something like
array[bp_n] real<multiplier = (length(bp_sd > 0) ? bp_sd[1] : 0)> bp_effects; // Rt breakpoint effects
Maybe though I have no idea of what a multiplier of 0 does?
This is a small PR that tests making the random walk non-centred.