Sampling rate when running independent regressions in vectorized form

Maybe try putting a hyper prior on coeff_pal and initialise the sampler at the MLE?