Power spectrum analysis with least-squares fitting: Amplitude bias and its elimination, with application to optical tweezers and atomic force microscope cantilevers
arXiv:0906.1708 · doi:10.1063/1.3455217
Abstract
Optical tweezers and AFM cantilevers are often calibrated by fitting their experimental powerspectra of Brownian motion. We demonstrate here that if this is done with typical weighted least-squares methods the result is a bias of relative size between -2/n and +1/n on the value of the fitted diffusion coefficient. Here n is the number of power-spectra averaged over, so typical calibrations contain 10-20% bias. Both the sign and the size of the bias depends on the weighting scheme applied. Hence, so do length-scale calibrations based on the diffusion coefficient. The fitted value for the characteristic frequency is not affected by this bias. For the AFM then, force measurements are not affected provided an independent length-scale calibration is available. For optical-tweezers there is no such luck, since the spring constant is found as the ratio of the characteristic frequency and the diffusion coefficient. We give analytical results for the weight-dependent bias for the wide class of systems whose dynamics is described by a linear (integro-)differential equation with additive noise, white or colored. Examples are optical tweezers with hydrodynamic self-interaction and aliasing, calibration of Ornstein-Uhlenbeck models in finance, models for cell-migration in biology, etc. Because the bias takes the form of a simple multiplicative factor on the fitted amplitude (e.g. the diffusion coefficient) it is straightforward to remove, and the user will need minimal modifications to his or her favorite least-square fitting programs. Results are demonstrated and illustrated using synthetic data, so we can compare fits with known true values. We also fit some commonly occurring power spectra once-and-for-all in the sense that we give their parameter values and associated error-bars as explicit functions of experimental power-spectral values.
20 pages, 10 figures
References in corpus (2)
Cited by in corpus (10)
- Harmonic Oscillator in Heat Bath: Exact simulation of time-lapse-recorded data, exact analytical benchmark statistics
- The mitotic spindle in the one-cell C. elegans embryo is positioned with high precision and stability
- Probing the screening of the Casimir interaction with optical tweezers
- What is superresolution microscopy?
- Measurement of probe displacement to the thermal resolution limit in photonic force microscopy using a miniature quadrant photodetector
- On the Estimation of Parameters from Time Traces originating from an Ornstein-Uhlenbeck Process
- Covariance Distributions in Single Particle Tracking
- Linking fluctuation and dissipation in spatially extended out-of-equilibrium systems
- Weighing an optically trapped microsphere in thermal equilibrium with air
- Robust and Efficient Parametric Spectral Estimation in Atomic Force Microscopy