Функция UniFrac
требует объекта phyloseq
:
physeq: (Required). ‘phyloseq-class’, containing at minimum a
phylogenetic tree (‘phylo-class’) and contingency table
(‘otu_table-class’). See examples below for coercions that
might be necessary
Не совсем уверен, имеет ли это смысл для вас, но, возможно, рассмотрим приведенный ниже пример. Ключ в том, чтобы правильно построить объект phyloseq, одна часть, в которой я не уверен, это как создать таблицу otu_table, я установил taxa_are_row=FALSE
:
#you get a tree
sample_tree = rtree(20,tip.label=letters[1:20])
# there are counts associated with each sample in tree
sample_counts = matrix(rnbinom(100*20,mu=20,size=1),
ncol=20,dimnames=list(1:100,letters[1:20]))
x1 = phyloseq(
otu_table(sample_counts,taxa_are_row=FALSE),
sample_tree)
UniFrac(x1)