Hola Juan.
findMRCA takes the tip labels - not the tip numbers (if I remember
correctly). That means that the vector tips should be a vector of class
'character'. So, even if your tip labels are 16 & 17, this should be
tips=c("16","17"). I'm not sure that this is your problem - but it might be.
All the best, Liam
Liam J. Revell, Assistant Professor of Biology
University of Massachusetts Boston
web: http://faculty.umb.edu/liam.revell/
email: liam.rev...@umb.edu
blog: http://blog.phytools.org
On 7/30/2013 9:47 AM, Juan Antonio Balbuena wrote:
Hello
I am using the function findMRCA of phytools. I tried it with a
simulated tree and it worked as it should. Now I want to use it with a
real tree but get an odd error message:
findMRCA(TreeP, tips= c(16, 17))
Error in if (anc == 0) break : missing value where TRUE/FALSE needed
This is the info about TreeP:
summary(TreeP)
Phylogenetic tree: TreeP
Number of tips: 17
Number of nodes: 16
Branch lengths:
mean: 0.05285941
variance: 0.001793694
distribution summary:
Min. 1st Qu. Median 3rd Qu. Max.
0.000001 0.013200 0.036010 0.084470 0.143000
No root edge.
First ten tip labels: Tbar
Tmin
Gtho
Gact
Gpero
Gche
Gcos
Gpan
Gset
Gcha
No node labels.
Any help will be much appreciated.
Juan Antonio Balbuena
--
Dr. Juan A. Balbuena
Marine Zoology Unit
Cavanilles Institute of Biodiversity and Evolutionary Biology
University of Valencia http://www.uv.es/~balbuena
<http://www.uv.es/%7Ebalbuena>
P.O. Box 22085 http://www.uv.es/cavanilles/zoomarin/index.htm
46071 Valencia, Spain http://cetus.uv.es/mullpardb/index.html
e-mail: j.a.balbu...@uv.es <mailto:j.a.balbu...@uv.es> tel. +34 963
543 658 fax +34 963 543 733
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
*NOTE!*For shipments by EXPRESS COURIER use the following street address:
C/ Catedrático José Beltrán 2, 46980 Paterna (Valencia), Spain.
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
_______________________________________________
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 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/