微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

使用谱系 = phylo 对象在 R 中运行 inverseA 函数时出错

如何解决使用谱系 = phylo 对象在 R 中运行 inverseA 函数时出错

我正在尝试计算系统发育相关性求和矩阵的逆(将输入 MCMCglmm)。我有一个系统发育树,最初运行如下:

inv.phylo <- inverseA(pedigree=tree,nodes="ALL",scale=TRUE,reduced = FALSE)

这引发了与我的系统发生树格式相关的错误,因此我对树进行了超度量并使用代码将树扎根:

tree.UM <- chronos(tree)
tree.root <- root(tree.UM,"Ovibos_moschatus",resolve.root=TRUE)

现在当我跑步时:

inverseA(pedigree=tree.root,reduced = FALSE)

我收到错误

Error in pedigree[,2] : incorrect number of dimensions

我注意到 treeUM 比树多一条边,但无法确定这是 1) 导致问题,还是 2) 如果是,该怎么办?到目前为止,我已经尝试反过来运行 chronosroot,并使用 chronoMPL 而不是 chronos。谁能帮我让 inverseA 运行?

提前致谢:)

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。