ncRNA num

最新推荐文章于 2014-02-17 09:06:49 发布

NeoBilly

最新推荐文章于 2014-02-17 09:06:49 发布

阅读量1k

点赞数

分类专栏：随手记

本文链接：https://blog.csdn.net/NeoBilly/article/details/14047439

版权

随手记专栏收录该内容

11 篇文章 0 订阅

订阅专栏

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