Fit framework

This page contains only a general description of FRETBursts fitting functionalities. The content of this page is:

For the reference documentation for fitting multi-channel histograms see:

Overview

FRETBursts uses of the powerful lmfit library for most fittings (like E or S histogram fitting). Lmfit should be automatically installed when installing FRETBursts, but in any case it is easily installable via pip install lmfit. For more installation info see FRETBursts Installation.

FRETBursts requires lmfit version 0.8 or higher.

Fitting E or S histograms

The module fretbursts.mfit provides a class fretbursts.mfit.MultiFitter that allow to build histograms and KDE on a multi-channel sample population (typically E or S values for each burst). The MultiFitter class can find the max peak position of a KDE or fit the histogram with an arbitrary model. A set of predefined models is provided to handle common cases. Sensible defaults are applied but the user can control every detail of the fit by setting initial values, parameter bounds (min, max), algebraic constrains and so on. New models can be created by composing simpler models (by using + operator). See the lmfit documentation for more info on how to define models and composite models.

A convenience function fretbursts.burstlib_ext.burst_fitter() can be used to create a MultiFitter object to fit either E or S. As an example let suppose having a measurement loaded in the variable d. To create a fitter object and compute the FRET histogram we execute:

bext.burst_fitter(d) # Creates d.E_fitter
d.E_fitter.histogram() # Compute the histogram for all the channels

Now we fit the E histogram with a 2-Gaussians model:

d.E_fitter.fit_histogram(mfit.factory_two_gaussians)

And plot the histogram and the fitted model:

dplot(d, hist_fret, show_model=True)

More detailed example can be found in the tutorials in notebooks on us-ALEX analysis.

Lmfit introduction

Lmfit provides a simple and flexible interface for non-linear least squares and other minimization methods. All the model parameters can be fixed/varied, have bounds (min, max) or constrained to an algebraic expression.

Moreover lmfit provides a Model class and a set of built-in models that allows to express curve-fitting problems in an compact and expressive form. Basic models (such as a Gaussian peak) and be composed allowing an easy definitions of a variety of models (2 or 3 Gaussians).

For more information refer to the official lmfit documentation.

Legacy Fit functions

A set of legacy functions used in versions of FRETBursts < 0.4 are defined in fretbursts/fit. This function are retained for backward compatibility but should not be used in new analysis.

These are low-level (i.e. generic) fit functions to fit gaussian or exponential models.