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
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