Finding Clocks in Genes: A Bayesian Approach to Estimate Periodicity

Figure 1

Flowchart of ABSR algorithm. is the observed profile of the th gene expression and is the standardized time profile from . is the detrended profile derived from and is the S-G filtered profile from . The AR models , , and are fitted for the detrended data by three methods of model parameter estimation: Yule-Walker method, Burg method, and MLE, respectively. The AR models , , and are fitted for the noise-reduced data by the three above-mentioned model parameter estimation methods, respectively. The frequency sets are obtained from the frequency spectra of the AR models , respectively. The harmonic models are constructed with and the candidate frequency sets . is the selected frequency set.