Inference of changes in effective HIV population size in Roka
We jointly inferred phylogenetic relationships, evolutionary rate, and demographic parameters among Roka HIV-Env sequences using three demographic models acting as coalescent priors on phylogeny structure: a Bayesian SkyGrid non-parametric model of population size change, an exponential growth model, and a model allowing a period of constant population size followed by exponential growth (constant + exponential).
We additionally conducted a sensitivity analysis to look at the robustness of population size inference under an exponential growth demographic model to changes in the clock rate prior, which is an informed prior.
The following types of files are available in this directory:
.xmlfiles are the input files for the joint inference performed in BEAST
.logfiles contain the logged draws from the posterior distribution for various parameter values
.tsvfiles contain the mean, median, upper HPD, and lower HPD values of Ne used to make the skyline plot
This other R script was used to plot the posterior distributions of NeTau under three different clock models.