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)