如何将引导值包含在ape包中的系统发育树中

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)

我想找到引导值，所以我使用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结果。