Estimating second order system parameters from noise power spectra using nonlinear regression
… the safe use of regression requires a good deal of thought and a good dose of skepticism
Overview
In the limits of detection lab, you will use regression to estimate the parameters of second order system models of the optical trap and AFM. MATLAB includes a powerful nonlinear regression function called nlinfit. Like all sharp tools, nlinfit must be used with extreme care. Regression is a delicate business.
This page illustrates several flawed regressions with the hope of helping you navigate the concealed hazards of fitting power spectra. Many common regression missteps spawn subtle errors rather than catastrophic failures. Because of this, bad regressions can masquerade as good ones. — hiding their insidious flaws behind a lovely, plausible pair of overlapping curves.
Remain vigilant at all times when you regress.
Regression review
Regression assumes that you are trying to find the relationship between a response variable $ y_n $ and an explanatory varialbe $ x_n $. The two quantities are related by a model function $ f(x_n,\Beta) $, where $ \Beta $ is a vector of model parameters. $ y_n $ is measured in the presence of noise, which is represented mathematically by a random variable $ \epsilon_n $. The basic regression model is:
$ y_n=f(x_n,\Beta)+\epsilon_n $,
where $ f $ is the model function, $ \Beta $ is a vector of model parameters and $ \epsilon $ is a random error term.