Phyloseq и тепловая карта

Я анализирую 16-секундные данные микробиома из легких и рта и в основном изучаю R. Поэтому я сделал это через qiime и загрузил два файла в R. Моя таблица OTU, используя:

otu=import_biom('C:\Users\OneDrive\otu_table.biom')

и мой файл «карты», который просто содержит SampleID, Location и Paired. Например:

 SampleID  Location   Paired
 1_L       Lung       1
 1_M       Mouth      1
 2_L       Lung       2
 2_M       Mouth      2

код:

map=read.csv('C:\Users\OneDrive\Mouth vs lung R map_adj.csv',
             header=T,row.name=1,stringsAsFactors=F)

Мои данные объединены, и я хочу сравнить легкое и рот, используя тепловую карту. У меня были проблемы с R. Во-первых, он не позволяет мне сделать это по местоположению. Я пытался:

plot_heatmap(otu, sample.label="SampleType")

И я получил эту ошибку:

Error in get_variable(physeq, sample.label) :
Your phyloseq data object does not have a sample-data component
Try ?sample_data for more details.

Я пошел в ?sample_data и понятия не имею, что мне здесь не хватает.

Во-вторых, я хотел бы сделать это по местоположению и в паре, если это имеет смысл.

Может ли кто-нибудь помочь мне с этим кодом и, возможно, объяснить, что мне здесь не хватает, потому что у меня также есть данные, в которых я смотрю на микробиом легких на исходном уровне, а затем снова через 2 месяца, так что это также парные данные.

Спасибо за помощь!

Вот весь мой код:

library("phyloseq")
packageVersion("phyloseq")

otu=import_biom('C:\Users\closed-ref-33031010\otu_table.biom')

map=read.csv('C:\Users\Qiime Maps\Mouth vs lung R map_adj.csv',
             header=T,row.name=1,stringsAsFactors=F)

sampledata1=sample_data(map)

otu2=merge_phyloseq(otu,sampledata1)
plot_heatmap(otu2)

Кроме того, я не уверен, как сделать объект phyloseq