Shaky model in least squares fit

Click For Summary

Discussion Overview

The discussion revolves around challenges faced in performing least squares fits on spectroscopic data, particularly concerning the robustness of the physical model used for fitting. Participants explore methods to quantify the reliability of model parameters and the significance of systematic errors compared to statistical noise.

Discussion Character

  • Exploratory, Technical explanation, Debate/contested

Main Points Raised

  • One participant notes that the physical model for the least squares fit is inadequate, leading to significant deviations between the model and the data that exceed statistical errors.
  • Another participant suggests quantifying systematic errors by calculating the F-value, which compares the variance of systematic errors to the variance of noise, and mentions the F-test for verifying model assumptions.
  • A different viewpoint proposes applying various models to the same data to observe how fits change, including fixing one coefficient at a time to re-estimate others.
  • One participant expresses gratitude for the suggestions and indicates a willingness to explore them further.
  • Another participant points to the Akaike information criterion as a relevant metric for assessing model quality.

Areas of Agreement / Disagreement

Participants present multiple approaches to address the issues raised, indicating a lack of consensus on a single method or solution. The discussion remains open with various suggestions and perspectives offered.

Contextual Notes

Participants acknowledge the limitations of the current model and the assumptions underlying least squares fitting, such as the independence and normal distribution of errors, which remain unresolved in the context of the discussion.

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   Reactions: 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   Reactions: 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.
 

Similar threads

  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 13 ·
Replies
13
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 20 ·
Replies
20
Views
3K
Replies
24
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K