0

嗨,我正在尝试重新排序 multiphylo 对象中的所有树,这应该只是在列表上应用函数的一种形式,我以前做过。但是,我似乎无法从包 APE 中获取“重新排序”命令来处理跨多门对象的 lapply 或 treeapply 。我的函数确实使用 lapply 跨越了 multiphylo 对象,但实际上并没有重新排序树。它似乎只是返回原始的无序树。任何帮助将非常感激!

我的代码示例

编辑:它似乎也不只是我的数据,我添加了一行代码来模拟 100 棵树,并且该函数也没有正确地重新排序模拟的树。

#I normally read my multiphylo object in from a file, but simulate one 
#here as an example
#trees<-read.trees("multiphyloobject.tre")

#Simulate 100 trees
trees<-pbtree(n=100,scale=10,nsim=100)

myfun <- function(x) {  
reorder(x,order="pruningwise")}

ordtrees<-lapply(trees,myfun)
4

0 回答 0