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
网友评论