KissABC.jl icon indicating copy to clipboard operation
KissABC.jl copied to clipboard

Factored prior over multivariate distributions

Open zsteve opened this issue 3 years ago • 1 comments

Thanks for this package! At the moment, the Factored distribution only applies to UnivariateDistributions. Is there any reason currently the code does not support MultivariateDistribution?

The reason I am interested here is that I would like to estimate transition matrices using ABC. In this setting, the appropriate choice of prior would be a Dirichlet distribution for each row of the transition matrix. Then the prior over transition matrices should be something like

Factored([Dirichlet(zeros(N)) for _ in eachrow(P)]).

As a hack, I've just extended the Factored type to allow MultivariateDistribution. This seems to work as-is. Would be curious to see whether this is something that can be incorporated into the package. https://github.com/zsteve/KissABC.jl/tree/multivariatefactor

zsteve avatar Nov 10 '21 04:11 zsteve

Hi steve, open up a PR and add a few tests for the various sampling strategies :) I will be very glad to take a look at your improvements

francescoalemanno avatar Nov 23 '21 07:11 francescoalemanno