Bayesian adaptive regression splines for hierarchical data.

Journal Article (Journal Article)

This article considers methodology for hierarchical functional data analysis, motivated by studies of reproductive hormone profiles in the menstrual cycle. Current methods standardize the cycle lengths and ignore the timing of ovulation within the cycle, both of which are biologically informative. Methods are needed that avoid standardization, while flexibly incorporating information on covariates and the timing of reference events, such as ovulation and onset of menses. In addition, it is necessary to account for within-woman dependency when data are collected for multiple cycles. We propose an approach based on a hierarchical generalization of Bayesian multivariate adaptive regression splines. Our formulation allows for an unknown set of basis functions characterizing the population-averaged and woman-specific trajectories in relation to covariates. A reversible jump Markov chain Monte Carlo algorithm is developed for posterior computation. Applying the methods to data from the North Carolina Early Pregnancy Study, we investigate differences in urinary progesterone profiles between conception and nonconception cycles.

Full Text

Duke Authors

Cited Authors

  • Bigelow, JL; Dunson, DB

Published Date

  • September 2007

Published In

Volume / Issue

  • 63 / 3

Start / End Page

  • 724 - 732

PubMed ID

  • 17403106

Electronic International Standard Serial Number (EISSN)

  • 1541-0420

International Standard Serial Number (ISSN)

  • 0006-341X

Digital Object Identifier (DOI)

  • 10.1111/j.1541-0420.2007.00761.x


  • eng