pints icon indicating copy to clipboard operation
pints copied to clipboard

Parameter transformation for Nested Sampling

Open chonlei opened this issue 4 years ago • 0 comments

Making a separate issue for using Transformation class in Nested Sampling, as mentioned in PR #1165:

Because we cannot sample from the transformed prior; we can sample only from the original prior and then transform the samples.

The evidence Z = \int L(x) p(x) dx, where L(.) is the likelihood, p(.) is the prior, and x is the model parameter. If we apply a transformation of y = g(x), then Z = \int L(x(y)) p(x(y)) (dx/dy) dy.

My understanding of nested sampling, from a quick read here, we need to first sample from the prior. But we when do the transformation, we can no longer sample parameters y from the transformed prior p(y)dy = p(x(y)) (dx/dy) dy (and then calculate the likelihood with L(x(y))). What we can do is sample y from the original prior p(x(y))dy, and the calculate the transformed likelihood L(y) = L(x(y)) dx/dy...

So that we can merge that PR and close #1070 first.

@ben18785 @martinjrobins @MichaelClerx comments/thoughts welcome!

chonlei avatar Jul 06 '20 13:07 chonlei