… if the test rejects, then one concludes that there is either a problem of model misspecification or a problem of a best not being a global maximum, but one cannot distinguish between these two. However, in many cases other techniques are available for model checking. … to the best of our knowledge no statistical method is available for checking whether a given root corresponds to a global maximum, even if the the model has been determined correct. This is the main reason for writing this article.

Thanks to physics, astrophysicists know evidently the chosen model is correct but the model could have a complex parameter space providing multiple solutions. I thought this test for a global maximum would be useful under this circumstance. It’s cheaper to run the fitting process with different initial values instead of searching full parameter space. On the other hand, the circumstance that the global maximum in unphysical part of the parameter space does not satisfy the regularity conditions.

]]>The global maximum could be located in “unphysical” part of parameter the space. Do these test work in such situation, i.e. when the parameter space is constrained within the physically acceptable solutions while the global max is not?

When you work with observations you typically think about how well the global max is constrained given the data. Within the “confidence” regions in the parameter space you can find a range of acceptable solutions, so even if you determine that yes, this is a global max, your goal is to collect such data that this range is small. so if within the constrained 90% intervals you have many “local” maxima does it really matter that you landed in one of them? I think you need to be sure that you can constrain the 90% intervals well.

]]>