Block RW Metropolis Hastings

Hi

Is there a package implementing a block RWMH algorithm?

With block, I mean partitioning the parameter vector. Preferably one where I can choose the partition ex ante, but one with random blocks every draw works as well.

Thanks for any suggestions!

I put this in the Julia Slack, but I’ll repost here for posterity.

You can sort of do this with Gibbs sampling:

spl = Gibbs(
   MH(:a1, :a2),
   MH(:a3, :a4)
)
sample(model, spl, n)

if you want static MH (unconditional proposal distributions). If you have a proposal matrix and want to do random-walk MH you should be able to do

spl = Gibbs(
   MH(Sigma1, :a1, :a2),
   MH(Sigma2, :a3, :a4)
)
sample(model, spl, n)

for proposal matrices Sigma1 and Sigma2.