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,
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,