您好,登錄后才能下訂單哦!
小編給大家分享一下如何使用awk篩選差異基因,相信大部分人都還不怎么了解,因此分享這篇文章給大家參考一下,希望大家閱讀完這篇文章后大有收獲,下面讓我們一起去了解一下吧!
命令行篩選差異基因
將圖片拖文件格式:
ID log2fc FDR gene1 1 0.001 gene2 -1 0.001 gene3 1 0.001 gene4 2 0.0001
篩選條件:log2fc的絕對值大于1,FDR<0.05
命令:
awk 'BEGIN{OFS=FS="\t"}{if(FNR==1) print $0; else {abs_log2fc=($2<0?$2*(-1):$2);if(abs_log2fc>=1 && $3<0.05) print $0;}}' de_gene.txt
$2代表的是第二列的log2fc,你的文件的log2fc在哪一列就用$第幾列,$3代表FDR,你的文件第幾列是也是寫$第幾列
篩選差異基因存儲到不同的文件
命令:
awk 'BEGIN{OFS=FS="\t"; up="up"; dw="dw";}{if(FNR==1) {print $0 >up; print $0 >dw;} else if ($3<0.05) {if ($2>=1) print $0 >up; else if($2<=-1) print $0 >dw;}}' de_gene.txt
以上是“如何使用awk篩選差異基因”這篇文章的所有內容,感謝各位的閱讀!相信大家都有了一定的了解,希望分享的內容對大家有所幫助,如果還想學習更多知識,歡迎關注億速云行業資訊頻道!
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。