OTUtable標準化方法是什么

這篇文章主要介紹“OTUtable標準化方法是什么”，在日常操作中，相信很多人在OTUtable標準化方法是什么問題上存在疑惑，小編查閱了各式資料，整理出簡單好用的操作方法，希望對大家解答”O(jiān)TUtable標準化方法是什么”的疑惑有所幫助！接下來，請跟著小編一起來學習吧！

OTUtable標準方法

1.抽平方法：

single_rarefaction.py -i pick_de_novo_otus/otu_table_clean.biom -o pick_de_novo_otus/otu_table_clean_rare.biom -d 2032

2.css方法或者deseq2方法

normalize_table.py -i otu_table_clean.biom -a CSS -o otu_table_clean_css.biom
normalize_table.py -i otu_table_clean.biom -a DESeq2 -z -o otu_table_clean_deseq2_norm_no_negatives.biom

標準化之后數(shù)據(jù)處理：

###############################################################################
#對OTU豐度進行標準化CSS方法 ：https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4010126/  http://qiime.org/scripts/normalize_table.html
mkdir normalized_otu_css/
normalize_table.py -i $workdir/5.pick_otu/pick_de_novo_otus/otu_table_clean.biom -a CSS -o normalized_otu_css/CSS_normalized_otu_table.biom
biom convert -i  normalized_otu_css/CSS_normalized_otu_table.biom -o  normalized_otu_css/CSS_normalized_otu_table1.txt --to-tsv --header-key taxonomy
#去掉只有界注釋的OTU
grep -v 'k; _; _; B; a; c; t; e; r; i; a'  normalized_otu_css/CSS_normalized_otu_table1.txt >normalized_otu_css/CSS_normalized_otu_table.txt
#注意biom格式轉(zhuǎn)變 json to hdf5  ：https://www.億速云.com/article/1314
biom convert -i normalized_otu_css/CSS_normalized_otu_table.txt -o  normalized_otu_css/CSS_normalized_otu_table_hdf5.biom --table-type="OTU table" --to-hdf5  --process-obs-metadata taxonomy
biom summarize-table -i normalized_otu_css/CSS_normalized_otu_table_hdf5.biom -o normalized_otu_css/CSS_normalized_otu_table_summary.txt
#對OTU豐度進行標準化 相對豐度
echo Summarize taxa
echo summarize_taxa:level 2,3,4,5,6,7 >taxa_summary_parm.txt
echo plot_taxa_summary:x_width 8 >>taxa_summary_parm.txt
echo plot_taxa_summary:bar_width 0.5 >>taxa_summary_parm.txt
echo plot_taxa_summary:chart_type bar >>taxa_summary_parm.txt
summarize_taxa.py -i normalized_otu_css/CSS_normalized_otu_table_hdf5.biom -o normalized_otu_css/absolute_abundance -a --level 2,3,4,5,6,7
summarize_taxa.py -i normalized_otu_css/CSS_normalized_otu_table_hdf5.biom -o normalized_otu_css/relative_abundance  --level 2,3,4,5,6,7
mkdir normalized_otu_css/group_otus_css
for i in  city loc country;do 
collapse_samples.py -b normalized_otu_css/CSS_normalized_otu_table_hdf5.biom -m $fastmap --collapse_fields $i --output_mapping_fp normalized_otu_css/group_otus_css/group.$i.txt --output_biom_fp normalized_otu_css/group_otus_css/otu_table_$i.biom
biom convert -i normalized_otu_css/group_otus_css/otu_table_$i.biom -o normalized_otu_css/group_otus_css/otu_table_${i}.txt --to-tsv --header-key taxonomy
summarize_taxa.py -i normalized_otu_css/group_otus_css/otu_table_$i.biom -o normalized_otu_css/group_otus_css/taxa_summary_relative.${i} -L 2,3,4,5,6,7
summarize_taxa.py -i normalized_otu_css/group_otus_css/otu_table_$i.biom -o normalized_otu_css/group_otus_css/taxa_summary_absolute.${i} -a -L 2,3,4,5,6,7
done

到此，關于“OTUtable標準化方法是什么”的學習就結束了，希望能夠解決大家的疑惑。理論與實踐的搭配能更好的幫助大家學習，快去試試吧！若想繼續(xù)學習更多相關知識，請繼續(xù)關注億速云網(wǎng)站，小編會繼續(xù)努力為大家?guī)砀鄬嵱玫奈恼拢?/p>