Still (and although easy to implement ) I don’t believe I’ll have to go to such lengths to marginalize since I’m working with likelihood here but it might worth a shot.
Before anything, thanks for your reply.
Yes, this is mainly a numerical problem since the physics here comes from the model in case (which I omitted). Nevertheless I would think this is standard in some areas of cosmology (as it is in many other areas of physics) therefore the question.
About...
Hello, I'm new to the world of "numerical" astronomy/cosmology and I've been studying maximum likelihood estimation for model tests. The (analytical) theory is somewhat easy to understand, but I've been struggling with the numerical aspect of the computations.
Suppose you have some model with a...