R帮助为MiRKAT创建树

时间:2016-05-03 14:23:05

标签: r dataframe tree phylogeny

我正在尝试使用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

上找到插图

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

1 个答案:

答案 0 :(得分:0)

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

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