Also, for a little more context, I’m using a modified version of this to generate data for Bayesian linear regression models. Another similar shape bug I’ve encountered is the resulting data looks wonky when (coincidentally) samples = n. I haven’t had a chance to figure out exactly in what way the data is wonky so I haven’t posted anything about it yet, but the data plots look off and then my downstream methods produce incorrect results. My current workaround if samples = n is to first generate n-1 samples and then separately generate the last sample.