pm.HalfStudentT.dist() failing to give correct shape

continued… σ=pm.HalfNormal(“σ”,2000,dims=‘species’) also works.