Bayesian estimation of HIV-1 dynamics in vivo

Statistical analysis of viral dynamics in HIV-1 infected patients undergoing structured treatment interruptions were performed using a novel model that accounts for treatment efficiency as well as total CD8+ T cell counts. A brief review of parameter estimates obtained in other studies is given, pointing to a considerable variation in the estimated values. A Bayesian approach to parameter estimation was used with longitudinal measurements of CD4+ and CD8+ T cell counts and HIV RNA. We describe an estimation procedure which uses spline approximations of CD8+ T cells dynamics. This approach reduces the number of parameters that must be estimated and is especially helpful when the CD8+ T cells growth function has a delayed dependence on the past. Seven important parameters related to HIV-1 in-host dynamics were estimated, most of them treated as global parameters across the group of patients. The estimated values were mainly in keeping with the estimates obtained in other reports, but our paper also introduces the estimates of some new parameters which supplement the current knowledge. The method was also tested on a simulated data set.