In our weekly meeting, Hou, Goodman, Foreman-Mackey, and I discussed Hou's results trying to evaluate the marginalized likelihood integral. Also in the room was Marshall, visiting for the week from Oxford. Hou's early attempts using a Gaussian trial function and importance sampling are not working very well, so I suggested we make a set of simple functional tests on easier problems and Goodman suggested that we plug the emcee ensemble sampling operations into some kind of nested sampling code, including possibly Brewer's Dnest.
For the rest of the day we worked on getting started on putting real gravitational lens models into the Tractor, but we got stuck on scipy and Astrometry.net compilation issues on Marshall's laptop. Argh! Open-source code is not always easy to use.
No comments:
Post a Comment