DESeq2/edgeR
DESeq2和edgeR是RNA-Seq數(shù)據(jù)分析中常用的差異表達分析軟件摇幻,這兩款軟件基于負二項分布模型剥汤,輸入count 數(shù)據(jù)進行差異分析挑社。其中DESeq2必須有重復(fù)才能運行最蕾,edgeR可以基于設(shè)置的離散度進行無重復(fù)比較分析碳蛋。
參考腳本
## 準備count表格
S.count.tsv
## 替換掉空格
sed 's/ /_/g' S.count.tsv > S.count.matrix
## 進行差異分析
perl run_DE_analysis.pl \
--matrix S.count.matrix \ # count表格
--method DESeq2 \ # 分析方法钧忽,可設(shè)置為edgeR
--samples_file sample.txt \ # 樣品分組信息
--output DESeq2_out # 輸出結(jié)果目錄
輸出結(jié)果文件為:S.count.matrix.A_vs_B.DESeq2.DE_results ;
可基于FDR 和 Foldchange 進行差異物種篩選毯炮。
常用的篩選標準:FDR < 0.05 && |log2(FC)|>1