libpulsatile is a C++ library that unifies several Bayesian pulsatile hormone models. Over the last 20 years, the initial single-subject, single-hormone model has been extended to multiple subjects and driver-response pairs of hormones, as well as with different prior distributions, most importantly for the individual pulse locations.

This C++ library creates standardized sampling algorithms (Gibbs, Modified Metropolis Hastings, Birth-Death) and data structures allowing for a unified code base on which to finalize each of these models. When complete, an R package frontend will be developed and submitted to CRAN. The development repository is on GitHub.

[ed. add manuscript links for each model]