Hi Nkaimcaudle,
Thanks for the response. The numbers are closer to the MLE estimate now. However, a few of the chains are still diverging and the acceptance probability is above the target. The rhat values are close to 1.
In regards to y being (150,), any comments on why not (150, 1). I understand y as a (column) vector.