Question

我需要一个非常简单的循环，它从multiphylo对象中获取每个系统发生树，并计算每棵树的Ic并将其放入数据帧中。对不起如果这么简单，我是R的新手，我无法理解！

    library(apTreeshape)
    library(phytools)
   multi_trees<- pbtree(b=0.5, d=0, n=200, t=NULL, scale=NULL, 
  nsim=50, type="continuous", extant.only=TRUE) ### simulate 50 trees stored as multiphylo object
    converted_tree <- as.treeshape(multi_trees[[nsim=1]]) ## each tree need to be converted to class treeshapes using following function

 Ic <- colless(converted_tree, norm = NULL) #And finally calculate Ic for each tree

Answer 1

循环很简单，所以我假设它是自我解释的。

    library(apTreeshape)
    library(phytools)

    # simulate 50 trees stored as multiphylo object
    multi_trees<- pbtree(b=0.5, d=0, n=200, t=NULL, scale=NULL, nsim=50, type="continuous", extant.only=TRUE)

    # make data.frame for results
    Ic <- matrix(nrow=50, ncol=2)
    Ic <- data.frame(Ic)
    colnames(Ic) <- c("sim", "Ic")

    # loop
    for (i in 1:50) {
        converted_tree <- as.treeshape(multi_trees[[i]])
        Ic[i,1] <- i # get simulation number
        Ic[i,2] <- colless(converted_tree, norm = NULL) # get Ic
    }

    Ic

如何编写一个循环来重复计算模拟次数的分数？

1 个答案: