for a k \times k matrix \Sigma, for \mathrm{Wishart}(\Sigma, \nu) is only defined for \nu > k. I had been playing around with the dimensions so I put an offset in to ensure I wouldn’t inadvertently sample from undefined likelihoods. The “right” way to do this would be with pm.Bound, I think.