search

liamrevell / phytools

GNU General Public License v3.0
207 stars 56 forks source link

getting uncertainty/distributions of ancestral states from simmap #96

Closed HedvigS closed 3 years ago

HedvigS commented 3 years ago

Hi Liam. I'm sorry to bother you, I have yet another query. I really appreciate the help in #94, now I'm able to grab the ancestral states from the multiSimmap object by running make.simmap() with multiple iterations, passing it through summary() and inspecting the ace value.

I was wondering if there is some kind of way of easily deriving the distribution over the different lists from make.simmap(.., nsim>1)rather than the summary, and if so if this could be understood as the certainty of the ace reconstruction. I know that uncertainty can go into the SIMMAP calculation, I guess I'm wondering if it can also be derived from it. I hope I'm making some kind of sense and will be very happy with a "no that doesn't make sense"-answer.

HedvigS commented 3 years ago

I think i'm looking for something that doesn't exist, I'll give up.