Dear Leonardo.
This is covered in Chapter 7, pp. 215-219 of my book with Luke Harmon
(https://press.princeton.edu/books/phylogenetic-comparative-methods-in-r),
but don't worry -- I also just posted a solution on my blog:
Dear group,
I am performing analyses using threshBayes from the phytools package. I was
wondering if someone has ever extracted the highest posterior density
intervals (HPD) values using this function. How could I make this?
Thanks in advance
--
Leonardo Matheus Servino
Pós-graduação em