This code and supporting Jupyter notebooks document our implentation of the O-C method for
RR Lyrae variables using a modified Hertzsprung method, as well as their Markov chain Monte Carlo
fitting using the emcee
package. We fit the O-C diagrams with a mixed model containing terms for
both the light-travel time effect, as well as linear period change.
Further details of the techniques are available in the original paper of Hajdu et al. (2021).
These routines were tested on a machine running Linux Mint, using the following versions of the required packages:
Python
3.7.10Numpy
1.20.1emcee
3.0.2Matplotlib
3.3.4Scipy
1.6.2scikit-learn
0.24.1tqdm
4.59.0Numba
0.53.1ChainConsumer
0.31.0