如何在ape包中的系统发育树中包含bootstrap值

Question

我使用R包ape来分析存储在DNAbin对象中的一些序列：

library(ape)
my.seq  <- read.dna("sequences.txt", format = "clustal")
my.dist <- dist.dna(my.seq)
my.tree <- nj(my.dist)

我想找到bootstrap值，所以我使用boot.phylo：

boot <- boot.phylo(my.tree, my.seq, FUN = function(xx) nj(dist.dna(xx)), B = 100)

但是我收到一条错误消息：

Error in if (drop[j]) next : missing value where TRUE/FALSE needed

知道这意味着什么，以及如何解决它？我试过谷歌搜索错误信息，我找不到任何东西。

Answer 1

您的if条件导致NA。

必须有TRUE或FALSE结果。