0

私は肺と口からの 16s マイクロバイオーム データを分析しており、基本的には R を独学しています。そのため、qiime を使用して作成し、R に 2 つのファイルをアップロードしました。私の 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オブジェクトの作り方がわからない

4

1 に答える 1