Bayesian deconvolution analysis of pulsatile hormone concentration profiles

Many hormones are secreted into the circulatory system in a pulsatile manner and are cleared exponentially. The most common method of analyzing these systems is to deconvolve the hormone concentration into a secretion function and a clearance function. Accurate estimation of the model parameters depends on the number and location of the secretion pulses. To date, deconvolution analysis assumes the number and approximate location of these pulses are known a priori. In this article, we present a novel Bayesian approach to deconvolution that jointly models the number of pulses along with all other model parameters. Our method stochastically searches for the secretion pulses. This is accomplished by viewing the set of parameters that define the pulses as a point process. Pulses are determined by a birth-death process which is embedded in Markov chain Monte Carlo algorithm. This idea originated with Stephens (2000, Annals of Statistics 28, 40-74) in the context of finite mixture model density estimation, where the number of mixture components is unknown. There are several advantages that our model enjoys over the traditional frequentist approaches. These advantages are highlighted with four datasets consisting of serum concentration levels of luteinizing hormone obtained from ovariectomized ewes.