we

wc -l bed.v4.lnc v4.bed.rmCoding.cnci


 ls fpkm | perl -ne 'chomp ; if($. != 1){ print "\t";} print $_;' > header




for dir in `ls fpkm|head -1`; do  awk -F "\t" 'NR != 1' fpkm/$dir/genes.fpkm_tracking | sort -k 1 | cut -f 1 >fpkm.loci;   done
for dir in `ls fpkm`; do  awk -F "\t" 'NR != 1' fpkm/$dir/genes.fpkm_tracking | sort -k 1 | cut -f 10 |paste fpkm.loci - > t.fpkm.loci ; mv t.fpkm.loci  fpkm.loci;   done


for dir in `ls fpkm|head -1`; do  awk -F "\t" 'NR != 1' fpkm/$dir/isoforms.fpkm_tracking | sort -k 1 | cut -f 1 >fpkm.rna;   done
for dir in `ls fpkm`; do  awk -F "\t" 'NR != 1' fpkm/$dir/isoforms.fpkm_tracking | sort -k 1 | cut -f 10 |paste fpkm.rna - > t.fpkm.rna ; mv t.fpkm.rna  fpkm.rna;   done


tail fpkm.loci | less -S


表达值图示化




 wc -l v4.bed.rmCoding.cnci bed.v4.lnc  loci.tracking loci.loci