Thank you for investigating where this error comes from. Yes, according to the PyMC example mentioned above, or the documentation of LKJCholeskyCov (see the sd_dist parameter), shape should be equal to n. Indeed, it seems there is an issue with the initialization based on model.initial_point(), which, unless I am mistaken, shouldn’t happen.