Based on the algorithm discussed in the materials and methods of this paper,
I created a cython
script that takes a time trace as input and then:
- Creates a window over a given time interval
- Uses linear regression (using the
cython_gsl
package) to fit the deterministic motion of data over the window. - Determines the average mean squared displacement from the linear regression over the window.
- Loop over all possible windows and apply steps 1-3.
See the IPython Notebook illustrating how to use the package in the doc folder.