[R-sig-phylo] generate liability/threshold plot from ancThresh

Hi folks, I'm trying to generate the thresholds and liabilities plot, including the posterior support values, from Revell (2014) Evolution 68(3), Figure 6B. Does anybody have a script for this? I see the liabilities output from ancThresh(), but no thresholds. Thanks! ~John John S. S. Denton,

Re: [R-sig-phylo] generate liability/threshold plot from ancThresh

Hi John. The posterior sample of the relative positions of the thresholds (if there is more than one) are output in the matrix "par", i.e.: mcmc<-ancThresh(...) mcmc$par ## contains the posterior sample of the thresholds To get a visualization of the posterior density you can, for instance,