私は肺と口からの 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オブジェクトの作り方がわからない