jBillou edited Estimating phases from signal.tex  about 9 years ago

Commit id: 8c2ba81945440da5bfe01f9d281ae43dc3045333

deletions | additions      

       

\end{equation}  where $\xi$ is a normally distributed random variable. The model for the cell cycle is similar but links the cell cycle area to the cell cycle phase. In order to distinguish between amplitude or baseline changes and phase deformations, we imposed that the amplitude and the baseline can only change on long timescale (about a day). More precisely we modeled the phase as a diffusion-drift process and the amplitude and background as Ornstein-Uhlenbeck processes.  Using that model we calibrated HMM's capable of inferring the two phases from the signals. signals (Methods \ref{sec:HMM_methods}).  Once the estimated phases are known for a given cell trace, the trajectory can be plotted in the $(\theta, \phi)$ plane and the velocity in both direction can be estimated by computing the discrete derivative of the trajectory (Figure \ref{fig:inference_of_circadian_phase}B). Given a large collection of such traces it is possible to estimate the velocity on the whole $(\theta, \phi)$ plane \cite{Rosenblum2001} and thus to reconstruct the dynamics of the system. Inferred phase in the data are shown in Figure \ref{fig:inference_of_circadian_phase}C-D.