1

我正在尝试根据我使用包中phyloCorrelogram的数据创建系统发育相关图phylosignal，以测试系统发育信号的存在。我的数据是所谓的phylo4d格式，叫做tree.

现在，当我运行时phyloCorrelogram(tree)，我返回以下错误：

library("phylobase")
library("ape")
library("phylosignal") # contains phyloCorrelogram()

> phyloCorrelogram(tree, ci.bs = 10, n.points = 10)
Error in boot::boot(X, function(x, z) moranTest(xr = x[z], Wr = prop.table(Wi[z,  : 
  no data in call to 'boot'

我已经在互联网上进行了详尽的搜索以寻找解决此问题的方法，但没有成功。

由于我使用的数据维度太大，无法在此处发布，因此我已将数据以.rda格式发布在dropbox上。

有谁知道缺陷吗？

4

1 回答 1

2

您没有与树关联的数据。您的树是一个 phylo4d 对象，它具有“树”信息但没有附加数据。你需要这样的东西

library(phylobase)
g1 <- as(geospiza_raw$tree, "phylo4")
geodata <- geospiza_raw$data
g2 <- phylo4d(g1, geodata)
pc <- phyloCorrelogram(g2)
于 2022-01-14T10:23:15.407 回答