symmetric positive definite domain about transformvariables.jl HOT 10 OPEN

I will add one (perhaps in the examples on how to use CustomTransform). But I am curious what you need it for, since in practice its Cholesky factor is much more useful for inference.

from transformvariables.jl.

I have a model where the prior is inverse Wishart.

from transformvariables.jl.

I added it as a test/example.

I will think about adding it as a regular transform, but I don't think Wishart is a good prior, so I would rather encourage using half-Cauchy for standard deviations and LKJ for the correlation matrix (apologies if you know about this and decided you still need Wishart).

from transformvariables.jl.

I agree, but i'm replicating someone else's results.

from transformvariables.jl.

and thanks for the example!

from transformvariables.jl.

@simonbyrne: just to let you know, I am merging an API change now. I have updated the example in the tests.

If you are using this library in production, you may want to pin to df4234b, or update your code.

As the library is currently in development and I am experimenting, there will be no deprecations or releases until it stabilizes. So you may want to use Github's watch functionality for notifications.

from transformvariables.jl.

from transformvariables.jl.

+1 to please add this. I don't think preference against the Inverse Wishart prior is a strong reason not to include this transform - appropriateness of the Wishart/IW prior depends on the user's likelihood. I'm using a Wishart it for estimating a symmetric positive definite matrix that isn't a covariance (i.e. my likelihood is not multivariate Gaussian with unknown Sigma or anything remotely near), where LKJ is not necessarily appropriate, and Wishart with identity hyperparameter is a natural choice due to its group invariance properties.

from transformvariables.jl.

So would you need a transformation from an n⋅(n+1)/2 vector of reals to an n × n positive definite matrix?

from transformvariables.jl.

Yes - just as in Stan - this would be fantastic!

from transformvariables.jl.