随手记
NeoBilly
这个作者很懒,什么都没留下…
展开
-
we
wc -l bed.v4.lnc v4.bed.rmCoding.cnci ls fpkm | perl -ne 'chomp ; if($. != 1){ print "\t";} print $_;' > headerfor dir in `ls fpkm|head -1`; do awk -F "\t" 'NR != 1' fpkm/$dir/genes.f原创 2013-10-27 00:05:22 · 667 阅读 · 0 评论 -
中期报告数目统计
for org in celegans chicken cow dm human mouse tair yeast zebrafish; doecho $org;hgsql -e 'select count(id) from '${org}'_gene;' ncrnadonefor org in human mouse; doecho原创 2014-02-14 17:32:21 · 581 阅读 · 0 评论 -
NONCODE class code file
for org in dm celegans chicken cow human mouse tair yeast zebrafish;dohgsql -e 'select * from '$org'_gene_class into outfile "/data/tmp/NONCODEv4_GeneClass/'$org'_gene_class"' ncrna ;hgs原创 2014-01-12 23:22:17 · 697 阅读 · 0 评论 -
相识天数计算器
php -r '$timeNum0 = mktime(12,0,0,4,18,2011); $timeNum1 = time(); $days = ($timeNum1 - $timeNum0) / (24 * 60 * 60); echo intval($days) + 1; 'php -r '$timeNum0 = mktime(12,0,0,4,18,2011); $timeNum1 =原创 2014-01-02 09:25:00 · 4649 阅读 · 0 评论 -
创建gtf下载
搞错俩, 麻利的改回来cat id.v4_xloc | perl -ne '@t = split(/\t/, $_); chomp @t; $nId = "NONHSAG" . substr("000000", 0 , 6 - length($t[0]."") ) . $t[0]; print $t[1] . "\t" . $nId . "\n";' | perl /leo原创 2013-12-31 15:01:38 · 1273 阅读 · 0 评论 -
创建gtf下载数据
cat id.v4_xloc | perl -ne '@t = split(/\t/, $_); chomp @t; $nId = "NONHSAG" . substr("000000", 0 , 6 - length($t[0]."") ) . $t[0]; print $t[1] . "\t" . $nId . "\n";' | perl /leofs/noncode/NONCODEv4/原创 2013-12-12 11:12:38 · 1070 阅读 · 0 评论 -
NONCODE 入库
head bed.v4.ncrna head -2 v4.stats | tail -1 | perl -ne '@t= split(/\s/, $_); chomp @t; for($i = 6; $i for org in bosTau6 galGal3 tair10 dm3 ce10 danRer7 sacCer3docd $orghead -2原创 2013-11-05 20:42:21 · 1547 阅读 · 0 评论 -
ncRNA num
less draw.shfor((i=1; i < 20; i ++))do echo "ping 192.168.1."$i; donecd /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 | gr原创 2013-11-02 00:14:57 · 1081 阅读 · 0 评论 -
noncode
cat cnci.result.file | cut -f 1 | sed -e 's/>//' > cnci.1.tconscat cnci.result.file | perl -ne '@t = split(/\t/,$_); $t[0] =~ />(.*)/; $tid = $1; $t[2] =~ /score: (.*)/; $score = $1; print $t原创 2013-10-28 11:21:04 · 1540 阅读 · 0 评论 -
文件 级别的并行 画图
fpkmTable=$1pNum=$2cat $fpkmTable | awk '{print $0 >> "'$fpkmTable'"NR%'$pNum';}'####CAT THE HEADERfor (( i = 0; i cat header $fpkmTable""$i > $fpkmTable"."$idone####START THE PROGR原创 2013-10-27 02:27:20 · 558 阅读 · 0 评论 -
倒数日
php -r ' error_reporting(0); $nowSecInt = time(); $deadLineSecInt = mktime(12,0,0,7,12,2014);$days = ($deadLineSecInt - $nowSecInt) / (24*3600); print $days . "\n";'距离毕业还有145天原创 2014-02-17 09:06:49 · 1003 阅读 · 0 评论