less draw.sh
for((i=1; i < 20; i ++))do echo "ping 192.168.1."$i; done
cd /leofs/noncode/NONCODEv4
rm -rf CNCI_tmp/
nohup perl CNCI.pl -f cow.fa -p 50 -l libsvm-3.0 -o result_cow -b &
ps aux | grep CNCI | awk '{print $2;}'
cat cnci.result.file | perl -ne '@t = split(/\t/, $_); $t[0] =~ />(.*)/; $tid = $1; $t[2] =~ /score: (.*)/; $sc=$1; print $tid . "\t" . $sc . "\n"; ' > cnci.v4
awk '$2 < 0' cnci.v4 | wc -l
ln -s /leofs/noncode/NONCODEv4/A.Thaliana/cnci/result_tair/cnci.v4 cnci.v4
cat cnci.v4 | awk '$2 > 0 {print $1;}' | perl $CMD/rmByKey.pl - bed.v4.rmCoding 4 > bed.v4.rmCoding.cnci
cd /leofs/noncode/NONCODEv4/Cow
for org in hg19 mm9 bosTau6 galGal3 tair10 dm3 ce10 danRer7 sacCer3 Other
do
ls $org/bed.v4.ncrna $org/len.v3.noBed
done | xargs wc -l
for org in hg19 mm9 bosTau6 galGal3 tair10 dm3 ce10 danRer7 sacCer3 Other
do
cat $org/bed.v4.ncrna | getBedSeqLength.sh - | awk '$2 > 200{ print $1;}' | perl $CMD/leftJoin.pl - 1 $org/bed.v4.ncrna 4 | cut -f 2,3,4,5,6,7,8,9,10,11,12,13 >$org/bed.v4.lncrna
wc -l $org/bed.v4.lncrna ;
cat $org/len.v3.noBed | awk '$2 > 200' | wc -l ;
done | xargs wc -l
cat t | awk '{print $1}'| perl $CMD/sum.pl -
lncRNA#: 210831
ncRNA#: 595854
ncRNA num
最新推荐文章于 2022-08-10 17:49:09 发布