search

tritemio / FRETBursts

Burst analysis software for smFRET. **Moved to OpenSMFS organization**
https://github.com/OpenSMFS/FRETBursts
GNU General Public License v2.0
16 stars 17 forks source link

Lifetime Fitting #48

Closed ncodina closed 6 years ago

ncodina commented 7 years ago

I was wondering if there is already some function implemented to calculate lifetimes from the TCSPC histogram?

tritemio commented 7 years ago

Not out-of-the-box, but you can use straight lmfit to do the fitting. You can find an example of fitting with IRF deconvolution using least-squares or maximum likelihood in the following notebook (it uses simulated data):

https://github.com/tritemio/notebooks/blob/master/Lifetime_decay_fit.ipynb

Also you can fairly simply compute a per-burst "mean lifetime" with:

http://fretbursts.readthedocs.io/en/latest/plugins.html#fretbursts.burstlib_ext.calc_mean_lifetime

This can be used, for example, to build the plots like mean lifetime vs ratiometric E (see Seidel group MFD analysis).

tritemio commented 7 years ago

The notebooks recently published in this repository show how to do lifetime fitting on multispot data using both tail fitting and IRF reconvolution of the lifetime histogram with different models (1, 2 or 3 exponentials).

The notebooks are pretty long, for the part relative to lifetime fitting see this section:

http://nbviewer.jupyter.org/github/tritemio/polimi_tcspc_16ch/blob/master/16-CH%20POLIMI-TCSPC%20Bursts%20Analysis-2015-09-08_s5_FRET.ipynb#Population-Lifetime-Fit

tritemio commented 6 years ago

This issue was moved to OpenSMFS/FRETBursts#6