Question

我正在尝试使用R包MiRKAT创建一个树，但由于文档很难理解，我希望有人在这里可以帮助我。

我的数据在结构中：

         Acaryochloris Achromobacter Acidiphilium Acidisphaera
Sample1             23             1            0            4                 
Sample2              0             0            2            0                 
Sample3              4             0            4            0                 
Sample4              0            30            0            5                 
Sample5             11             0            5            0                 
Sample6              0             0            0            0

包裹在小插图中采用了结构树：

List of 4
 $ edge       : int [1:1710, 1:2] 857 858 859 860 861 862 863 864 865 866 ...
 $ Nnode      : int 855
 $ tip.label  : chr [1:856] "1883" "3114" "1483" "2576" ...
 $ edge.length: num [1:1710] 0.00846 0.09451 0.01012 0.01208 0.03501 ...
 - attr(*, "class")= chr "phylo"

可以在http://research.fhcrc.org/content/dam/stripe/wu/files/MiRKAT/MiRKAT_Vignette.pdf

上找到插图

有没有办法将我的数据转换为这样的树结构？如果是这样，你能推荐任何套餐吗？

Answer 1

对不起伙计们，我发现了自己。

mydata_dist <- dist(mydata)
mydata_hc <- hclust(mydata_dist)
mydata_phyl <- as.phylo(mydata_hc)

R帮助为MiRKAT创建树

1 个答案: