2020-10-09-cufflinks转录组分析结果FPKM值合并

cd ~/RNAseq/result/out_file/  #进入cufflinks结果输出目录

vim merge.sh

find ./ -name "gene*" |xargs -I {} awk '{print $1}' {} >track_id.list

for i in *_transcripts; do awk '{print "\t"$10}' $i/genes.fpkm_tracking >$i.txt && sed -i s/FPKM/$i/g $i.txt && paste track_id.list *.txt >all.fpkm; done

rm *.txt

rm track_id.list

保存并退出。

运行：sh merge.sh