I’m writing “toy” Gibbs sampler and my idea was to use LKJCholeskyCov for covariance matrix prior in
“stand alone mode” (without pymc model). But seems it has no “random” method like others distributions. So my question is why and is it possible to use LKJ distribution outside of a model.
Thanks a lot.