deletions | additions       diff --git a/Reconstruction coupling functions.tex b/Reconstruction coupling functions.tex  index 29deb5f..e61dfc8 100644  --- a/Reconstruction coupling functions.tex  +++ b/Reconstruction coupling functions.tex 
...
The velocities are computed by taking the finite difference of the phases $v_{\theta}(\phi,\theta) = (\theta(t+\Delta t)-\theta(t))/\Delta t$, with $\Delta t$ begin equal to $0.5h$. We estimated the functions $F1$ and $F_2$ by doing a least square fit of a 2D Fourier serie with 10 harmonics on the collection of velocities $\{v_{\theta}\}$ and $\{v_{\phi} \}$. Before doing so we verified that the distribution of velocities at every position in the phase plane was mainly unimodal. The functions were estimated on a $40$ by $40$ grid, given by the discrete hidden states of our HMMs.  In total we selected $2753$ time traces with a minimum length of $24h$ and at least two divisions. These gave $208'762$ velocities for each function. The data density (number of data points per discrete position in the phase plane) is shown in \fref{fig:dataDensity}. Figure \ref{fig:dataDensity}.