Re: [R-sig-phylo] Estimating the posterior prob of the ancestral states of two nodes

2023-03-05 Thread Liam J. Revell
Dear Sishuo.

Yes, I think this would be a good way to do this.

All the best, Liam

Liam J. Revell
Professor of Biology, University of Massachusetts Boston
Web: http://faculty.umb.edu/liam.revell/
Book: Phylogenetic Comparative Methods in R 
 
(/Princeton University Press/, 2022)


On 3/5/2023 9:49 PM, Sishuo Wang wrote:
> CAUTION: EXTERNAL SENDER
> Dear all,
>
> I would like to ask for a little help from you.
>
> I wish to estimate the posterior probability that *two* ancestral 
> nodes take particular states, given the observations of discrete trait 
> at the tips of a tree. I was wondering if using a *SCM* approach 
> described below would work.
>
>
> Inline image
>
>
> If so, I think this can be easily done using /make.simmap/ from 
> phytools where I would like to assign an equal prior on the root to 
> all states.
>
> Does the above make sense? Any insights are highly appreciated!
>
> Thanks,
> Sishuo
> 
> Postdoctoral Fellow, Haiwei Luo Lab
> Room 103A, F. S. Li Marine Science Laboratory, The Chinese University 
> of Hong Kong
> Shatin, New Territories, Hong Kong
>
>
> ___
> R-sig-phylo mailing list -R-sig-phylo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
> Searchable archive athttp://www.mail-archive.com/r-sig-phylo@r-project.org/
[[alternative HTML version deleted]]

___
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/


[R-sig-phylo] Estimating the posterior prob of the ancestral states of two nodes

2023-03-05 Thread Sishuo Wang
Dear all,
I would like to ask for a little help from you.

I wish to estimate the posterior probability that two ancestral nodes take 
particular states, given the observations of discrete trait at the tips of a 
tree. I was wondering if using a SCM approach described below would work.




If so, I think this can be easily done using make.simmap from phytools where I 
would like to assign an equal prior on the root to all states.
Does the above make sense? Any insights are highly appreciated!
Thanks,SishuoPostdoctoral
 Fellow, Haiwei Luo Lab
Room 103A, F. S. Li Marine Science Laboratory, The Chinese University of Hong 
KongShatin, New Territories, Hong Kong
___
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/