Abstract
This study proposes a Bayesian joint model with extended random effects structure that incorporates nested repeated measures and provides simultaneous inference on treatment effects over time and drop-out patterns. The proposed model includes flexible splines to characterize the circadian variation inherent in blood pressure sequences, and we assess the effectiveness of an intervention to resolve pediatric obstructive sleep apnea. We demonstrate that the proposed model and its conventional two-stage counterpart provide similar estimates of nighttime blood pressure but estimates on the mean evolution of daytime blood pressure are discrepant. Our simulation studies tailored to the motivating data suggest reasonable estimation and coverage probabilities for both fixed and random effects. Computational challenges of model implementation are discussed.