You can download the R source code for the Jive package here.
This package includes functions to analyse phenotypic data using a hierarchical Bayesian model that accounts for intraspecific variation. The publication is in preparation for Systematic Biology (Kostikova et al. in prep). An example of the application is in press in Proc R Soc B (Litsios et al. in press).