Dear John,
can you please be a bit more specific with your error message. It is
always good to have a reproducible example, e.g. adding a tree which
where drop.tip fails, and to run traceback() just after the error to
get more information where the error occurred. It is also useful to
add information of the version of ape and your operating system.
Regards,
Klaus
On 2/28/13, john d <dobzhan...@gmail.com> wrote:
> Dear all,
>
> I'm trying to prune a set of 1000 post-burnin trees to include only a
> subset of taxa. Unfortunately the tree is too big to send to the list,
> but if it is really necessary I'll figure out a way to do it.
>
> tr is my tree and taxa is my list of selected terminals.
>
> for(i in 1:1000){
> write.tree(drop.tip(tr[[i]],tr[[i]]$tip.label[-match(taxa,
> tr[[i]]$tip.label)]), file="result.tre", append=TRUE)
> }
>
> If I run that code, it works for some trees, but not for others, for
> which I got the message "Error in kids[[parent[i]]] : subscript out of
> bounds".
>
> Any suggestions?
>
> John
>
> _______________________________________________
> 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/
>
--
Klaus Schliep
Phylogenomics Lab at the University of Vigo, 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/
