Hi Sam, Santiago,
To be clear, this procedure yields the "squared-change parsimony"
ancestral states for a continuous character (see Maddison, 1991; Syst.
Zool.). This will not provide MP ancestral state estimates for a
discrete character, if this is what you are looking for.
Note that Maddison (1991) also proposed "weighted squared-change
parsimony" which will provide exactly the same ancestral state estimates
as ML under Brownian motion.
- Liam
Liam J. Revell
NESCent, Duke University
web: http://anolis.oeb.harvard.edu/~liam/
NEW email: lrev...@nescent.org
Samantha Price wrote:
Hi Santiago,
Look at the ace function in the Ape package - to get a max parsimony
reconstruction set branch lengths =1 and run a maximum likelihood
estimation. You can find further information here:
http://www.r-phylo.org/wiki/HowTo/Ancestral_State_Reconstruction
Sam
On Jan 23, 2010, at 8:32 PM, Santiago Claramunt wrote:
Hi all,
Does anybody know a function in R for performing maximum parsimony
ancestral reconstructions?
Best,
Santiago
Santiago Claramunt
Museum of Natural Science,
119 Foster Hall,
Louisiana State University,
Baton Rouge, LA70803
scla...@tigers.lsu.edu
http://www.museum.lsu.edu/Claramunt/Home.html
_______________________________________________
R-sig-phylo mailing list
R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
_______________________________________________
R-sig-phylo mailing list
R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
_______________________________________________
R-sig-phylo mailing list
R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo