flist=list.files("./",pattern = ".csv")
lapply(flist,function(csvname){
tree1 <- read.csv(csvname,header = T)
tree1 <- phylo.maker(sp.list=tree1,tree=GBOTB.extended,nodes = nodes.info.1,scenarios = "S3")
write.tree(tree1$scenario.3,"tree.tre")
tree <- read.tree("tree.tre")
sigdat <- read.csv(csvname,row.names = 1)
name.check(tree,sigdat)
data1 <- sigdat[,4]
names(data1) <- rownames(sigdat)
sink("K.csv",append = TRUE)
K <- phylosig(tree,data1,method="K",test=T)
print(K)
sink()
})
Error in if (simK >= K) P <- P + 1/nsim :
需要TRUE/FALSE值的地方不可以用缺少值
上边是我写的一段代码,然后运行的时候中途报错了,之前运行的时候没有问题,这次出问题了,我检查了我的数据,格式什么的跟原来一样的,求教大佬们哪里错了
我点了rerun with debug他还是运行完了,这是什么情况