R Помощь в создании дерева для MiRKAT

Я пытаюсь создать дерево для использования с пакетом 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 ответ

Решение

Извините, ребята, я узнал сам.

mydata_dist <- dist(mydata)
mydata_hc <- hclust(mydata_dist)
mydata_phyl <- as.phylo(mydata_hc)
Другие вопросы по тегам