3

After solving the inverse problem*, how to investigate whether the solution is unique or not? When I check theoretical values and published values during the convergence and search for the lowest value for chi-squared, is it sufficient?

*The task when the observed data is given (radial velocity, light curves for an eclipsing binary), and someone needs to determine the set of parameters for which the model yields the best match of the observed curves.

Anna-Kat
  • 505
  • 2
  • 6
  • I think you mean "when identifying the causes of a set of observations"... "how to I determine if my solution is unique. Is that right? – James K Aug 30 '21 at 07:20
  • 1
    I mean that the observed data is given (radial velocity, light curves for an eclipsing binary), and someone needs to determine the set of parameters for which the model yields the best match of the observed curves. – Anna-Kat Aug 30 '21 at 07:21
  • @JamesK yes it is right – Anna-Kat Aug 30 '21 at 07:24
  • 1
    In general you can't. An method based on convergence can hope at best to find a local minimum (of chi-squared or however you are judging the goodness of fit) There might be other, better minima that your method doesn't find. – James K Aug 30 '21 at 07:25
  • So the only inspection is to try whether slightly different initial values led to a similar solution and check that the resulting values make sence, right? – Anna-Kat Aug 30 '21 at 07:28

1 Answers1

2

In general you can't prove that your solution is unique. Methods for fitting a model to data that are based on some form of convergence can get caught in a local optimum and not find the global optimum.

For specific models with specific criteria the best model can be found analytically. For a simple example, fitting a linear model with least squares optimisation can be found by the usual linear regression line. No convergence is used and the method provably gives the best fit according to the given measure.

But for non-linear model fitting, there are no general proofs, and it is quite possible for a particular convergence to fail.

However "the best" model may not be required. You may only need a "good enough" model.

James K
  • 120,702
  • 5
  • 298
  • 423