I Shaky model in least squares fit

AI Thread Summary
The discussion centers on a least squares fit of spectroscopic data that reveals significant deviations between the model and measurements, indicating a mediocre physical model. The user seeks to quantify the robustness of certain model parameters, which vary when fitting subsets of the data. Suggestions include using the F-value to assess systematic errors versus noise and applying the F-test to validate model assumptions. Additionally, exploring different models and fixing coefficients can help evaluate the fit's sensitivity. The Akaike information criterion is identified as a useful tool for model comparison.
Gigaz
Messages
110
Reaction score
37
I've come across a problem with my least squares fits and I think someone else must have analyzed this, but I don't know where to find it.

I have a converged least squares fit of my spectroscopic data. Unfortunately, the physical model, on which the fit is based, is mediocre. The deviations between measurement and model are much larger than the statistical errors at each data point. There is almost certainly nothing I can do about that. The fit reproduces the data reasonably well, but the model is incomplete.

I know that there are some parameters inside the model, which do not seem to be very robust. If I fit only half of my data (only s or only p polarization), they always come out differently. Other parameters remain totally unchanged.

I'm looking basically for an idea on how I could quantify this "robustness". It can probably been done based on some sort of artificial perturbation function, but I haven't seen anything like that.
 
Physics news on Phys.org
Hi Gigaz, welcome to PF,

Apparently your systematic error (error with respect to the proposed model) is much larger than the error due to noise.
To quantify that we can divide the variance in the systematic errors by the variance of the noise.
This is called the F-value.
We can use the F-test to verify how significant this is with the hypothesis that a proposed regression model fits the data well.

For the record, a least squares method assumes that the errors are independent, normally distributed, have equal variance everywhere, and have expectation zero.
The F-test can verify part of those assumptions.
 
  • Like
Likes Gigaz and berkeman
A case of data in search of a model, it seems. You can apply the same data to different models, and see how that changes the fit. For example, you can fix one coefficient at a time, and re-estimate the remaining parameters.
 
  • Like
Likes Gigaz
Many thanks for those suggestions. I will try it and see what comes out :)
 
Thanks, chiro. Apparently, what I need is listed in your link: The Akaike information criterion.
 
I was reading documentation about the soundness and completeness of logic formal systems. Consider the following $$\vdash_S \phi$$ where ##S## is the proof-system making part the formal system and ##\phi## is a wff (well formed formula) of the formal language. Note the blank on left of the turnstile symbol ##\vdash_S##, as far as I can tell it actually represents the empty set. So what does it mean ? I guess it actually means ##\phi## is a theorem of the formal system, i.e. there is a...
Back
Top