| ²é¿´: 1271 | »Ø¸´: 8 | ||
Àгæ (³õÈëÎÄ̳)
|
[ÇóÖú]
perl ´óÉñÇóÖú ÒÑÓÐ1È˲ÎÓë
|
|
Ç×ÃÇ£¬ÎÒµÄÁ½¸öÎļþÊÇÕâÑùµÄ 1Îļþchr1 - 9901.5 chr1 - 144612.5 chr1 - 227317.5 chr1 - 227321.5 chr1 + 267647.5 chr1 - 267648.5 chr1 - 267649.5 chr2 + 267651.5 chr2 - 267652.5 chr3 - 521268.5 chr3 - 569831.5 chr1 + 713294.5 chr1 - 713307.5 chr1 - 753396.5 chr1 + 762732.5 chr3 - 801094.5 chr1 + 808622.5 chr2 - 808633.5 chr1 - 819067.5 3Îļþ chr1 - 9802 144712 chr1 - 227222 227421 chr1 - 267552 267751 chr2 + 267552 267753 ´¦ÀíÉϱßÁ½¸öÎļþµÄÔ¤ÆÚÄ¿±êÊÇ chr1 - 9802 144712 3 chr1 - 227222 227421 2 chr1 - 267552 267751 2 chr2 + 267552 267753 0 ÎÒÏëÔÚÁ½¸öÎļþµÚÒ»Áк͵ڶþÁÐÏàͬµÄÇé¿öÏ£¬Èç¹û1ÎļþµÄµÚÈýÁÐÓÐÒ»¸ö´¦ÓÚ3ÎļþµÄµÚÈýÁк͵ÚËÄÁÐÖ®¼äµÄÇé¿öÏ£¬$n¾Í¼Ó1¡£ÓÐÁ½¸ö¾Í¼Ó2£¬ÒÔ´ËÀàÍÆ£¬$n³õʼֵΪ0.ÆäÖеÄ3,2£¬2,0ÊÇ$n. my %hash; my $n=0; open OUTFILE,'>','result1'; open FILE,'<','1'; while (<FILE> {chomp; ($a,$b,$c)=(split)[0,1,2]; $hash{$a}=1; $hash{$b}=2; open FILE,'<','3'; while (<FILE> {chomp; ($d,$e,$f,$g)=(split)[0,1,2,3]; if(exists $hash{$d} && exists $hash{$e}){ @array=sort{$a <=> $b}($f,$g); next if $c<$array[0] ; next if $c>$array[1]; $n++; }print OUTFILE"$d $e $f $g $n\n"; } } ÕâÊÇÎÒ±àдµÄ³ÌÐò£¬´ó¼ÒÄܲ»ÄܰïÎÒ¿´¿´ÊÇÄÄÀï³öÁËÎÊÌâ°¡»òÕß°ïÎÒÖØÐ±àдһ¸ö£¬ÉõÊǸÐл¸÷λ´óÅ£ÁË¡£ |
» ²ÂÄãϲ»¶
»¯Ñ§¹¤³Ìµ÷¼Á289
ÒѾÓÐ42È˻ظ´
282£¬Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
»¯¹¤µ÷¼ÁÇóµ¼Ê¦ÊÕÁô£¡Ò»Ö¾Ô¸Ê§Àû£¬Ì¤Êµ¿Ï¸É£¬ÓÐÖ²ÎïÌáÈ¡¿ÆÑоÀú
ÒѾÓÐ5È˻ظ´
277 ÊýÒ»104£¬Ñ§Ë¶£¬Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
087100³õÊÔ311Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
297Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
»¯Ñ§¹¤³ÌÓë¼¼Êõ324µ÷¼Á
ÒѾÓÐ10È˻ظ´
085404£¬285·ÖÇóµ÷¼Á
ÒѾÓÐ12È˻ظ´
0856ר˶Çóµ÷¼Á Ï£ÍûÊÇaÇøÔºÐ£
ÒѾÓÐ9È˻ظ´
Ò»Ö¾Ô¸211µç×ÓÐÅÏ¢347Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
ÂçºÏÌúÏà¹ØÎÊÌ⣬Çó´óÉñ´Í½Ì£¡
ÒѾÓÐ3È˻ظ´
¾Ûõ£Ñǰ·±¡Ä¤¼ì²âÏîÄ¿ÖеÄÄÍÈÈÐÔ±ê×¼ÎÊÌâ
ÒѾÓÐ18È˻ظ´
ÏëѯÎÊһϡ°Ä³»ùÒò¡±ºÍ¡°Ä³»ùÒò-like¡±µÄ²î±ð? ¹òÇó´óÉñÏêϸ½â´ð£¡
ÒѾÓÐ4È˻ظ´
¸÷λ´óÉñÇó´Í½Ì
ÒѾÓÐ4È˻ظ´
±à³Ì²ËÄñÇó´óÉñÒýµ¼ÈëÃÅ¡£
ÒѾÓÐ14È˻ظ´
¡¾ÇóÖú¡¿Çë´óÉñ°ïÎÒ·ÒëÒ»¾ä»°
ÒѾÓÐ3È˻ظ´
ÉúÎïÐÅϢѧÎÊÌ⣬ȾɫÌå»û±ä²éÕÒ
ÒѾÓÐ9È˻ظ´
ÐÞ¸Äͼ±êÎÊÌ⣬Çë¸÷λ´óÉñ°ï棡
ÒѾÓÐ5È˻ظ´
±à³ÌÓеãÎÊÌ⣬ÇóÖú
ÒѾÓÐ5È˻ظ´
Çó´óÉñÃÇÑ©ÖÐËÍÌ¿
ÒѾÓÐ12È˻ظ´
ÇóÖú´óÉñ°ïÎÒ¿´¿´Õâ¶Î³ÌÐòÊǸÉʲôµÄ£¿£¨¾ÍÊǶ¨ÒåÕâ¸öº¯ÊýËãʲô£©
ÒѾÓÐ6È˻ظ´
ÇóÖú Ò»¸öperl´úÂëÐÞ¸Ä~~¿´¿´°É...
ÒѾÓÐ6È˻ظ´
Çó´óÉñ×öÒ»¸ö»ã±àÓïÑÔ³ÌÐò
ÒѾÓÐ7È˻ظ´
perl³ÌÐòÇóÖú£¬Çë´óÉñÃǴͽÌ~~~
ÒѾÓÐ9È˻ظ´
Rg(x) vs time ,Rg(y) vs time,Rg(z) vs time
ÒѾÓÐ6È˻ظ´
Çë°ïæ±àдһ¸öperl³ÌÐò´¦ÀíÎı¾
ÒѾÓÐ26È˻ظ´
ÇóCÓïÑÔ´óÉñ°ïæ
ÒѾÓÐ6È˻ظ´
Çë½Ìperl½Å±¾µÄÒ»¸öifÓï¾ä½âÊÍ£¬Ð»Ð»
ÒѾÓÐ8È˻ظ´
Perl ±à³Ì Çë´óÉñÃǰïæ£¡£¡£¡ÉúÎïÐÅÏ¢µÄÅ£ÈËÇë½ø
ÒѾÓÐ11È˻ظ´
perl½Å±¾ÇóÖú£¬ÎÒÏ뽫һÍò¶à¸öÎļþ´Ófasta¸ñÐò£¬¸ÄΪphylip4¸ñʽ
ÒѾÓÐ8È˻ظ´
wizardfan
ÖÁ×ðľ³æ (ÖøÃûдÊÖ)
- BioEPI: 18
- Ó¦Öú: 599 (²©Ê¿)
- ¹ó±ö: 1.818
- ½ð±Ò: 24632.2
- É¢½ð: 197
- ºì»¨: 48
- ɳ·¢: 2
- Ìû×Ó: 2254
- ÔÚÏß: 400.7Сʱ
- ³æºÅ: 1879241
- ×¢²á: 2012-07-05
- ÐÔ±ð: GG
- רҵ: ÉúÎïÐÅϢѧ
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
À: ½ð±Ò+15, ¡ïÓаïÖú 2014-11-25 09:42:42
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
À: ½ð±Ò+15, ¡ïÓаïÖú 2014-11-25 09:42:42
|
Ê×ÏÈÏÈÈÃÎÒÀí½âÄãµÄËã·¨£º ¸ù¾ÝÄã˵µÄ chr1 - 9802 144712 µÚ3£¬4ÁÐÊÇ9802ºÍ144712£¬ÄÇÇø¼ä¾ÍÊÇ9802-144712 ÔÚ1ÎļþÖУ¬Ö»ÓеڶþÐÐchr1 - 144612.5 ÔÚÕâ¸öÇø¼ä£¬ÄÇÄãµÄ½á¹ûΪʲôÊÇ3ÄØ£¿ Æä´Î£¬ÇëÈ·¶¨1ÎļþºÍ3ÎļþµÄ·Ö¸î·ûÊǿոñ£¬»¹ÊÇtab ½Ó×Å ¿´ÉÏÈ¥µÚ¶þÁеÄÕý·´Á´ÐÅÏ¢ºÃÏñûÓÐÓÃµÄ ×îºó ÇëÔÚÄãµÄ³ÌÐòÉϼÓÉÏ×¢ÊÍ£¬ÎÒ¿ÉÒÔ°ïÄã¿´¿´£¬Äãµ½µ×´íÔÚÄÄÀ ÊÇËã·¨´íÎ󣬻¹ÊÇperlʵÏÖ´íÎó |
2Â¥2014-11-10 17:36:47
À
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 621.5
- É¢½ð: 4
- Ìû×Ó: 30
- ÔÚÏß: 15.2Сʱ
- ³æºÅ: 2652634
- ×¢²á: 2013-09-13
- רҵ: Ô×ӺͷÖ×ÓÎïÀí
|
²»ºÃÒâ˼£¬×òÌì°Ñ1ÎļþÉÏ´«´íÁË£¬Èç¹û°´Ò»ÎļþµÄ»°£¬Ô¤ÆÚµÄÄ¿±êÊÇ chr1 - 9802 144712 2 chr1 - 227222 227421 2 chr1 - 267552 267751 2 chr2 + 267552 267753 1 ÎÒ½âÊÍһϣ¬ÎÒÏë¸Éʲô(ÉϱßÎļþ1µÄ×ͷ»¹ÓÐÒ»ÐÐchr1 - 9901.5 )£¬ÎÒÁí$nµÄ³õʼֵΪ0£¬Èç¹ûÎļþ1µÄµÚÈýÁÐÓÐÂäÔÚÎļþÈýµÄµÚÈýÁк͵ÚËÄÁÐÖ®¼ä£¬ÄÇô$n¾Í¼Ó1£¬ÒÔ´ËÀàÍÆ£¬ chr1 - 9901.5 chr1 - 144612.5 ÈçÉϱßÁ½ÐеĵÚÈýÁж¼ÔÚ9802 144712 Ö®¼ä£¬ËùÒÔ$n¾ÍΪ2.ͬÑùϱßÒ²ÊÇÕâÑùµÄÒ²¾ÍÊÇÉúÎïÖÐÊýread¸öÊý¡£¿ÉÒÔ°ïÎÒ±àдһ¸öÒ²ÊÇ¿ÉÒԵģ¬ÒòΪ°ï±ðÈ˸ijÌÐòȷʵºÜÂé·³µÄ£¬·Ç³£¸Ðл |
3Â¥2014-11-11 08:34:54
wizardfan
ÖÁ×ðľ³æ (ÖøÃûдÊÖ)
- BioEPI: 18
- Ó¦Öú: 599 (²©Ê¿)
- ¹ó±ö: 1.818
- ½ð±Ò: 24632.2
- É¢½ð: 197
- ºì»¨: 48
- ɳ·¢: 2
- Ìû×Ó: 2254
- ÔÚÏß: 400.7Сʱ
- ³æºÅ: 1879241
- ×¢²á: 2012-07-05
- ÐÔ±ð: GG
- רҵ: ÉúÎïÐÅϢѧ
|
ÎÒ×Ô¼ºÐ´ºÜ¼òµ¥£¬µ«ÊÇÄãѧµ½µÄÉÙ¡£¸üÇãÏòÓÚÄã˵ÏÂÄãµÄËã·¨¡£ ±ÈÈç ¶ÁÈëÎļþ1£¬Ã¿ÐзֳÉ3ÁУ¬±£´æµ½hash±äÁ¿ÀkeyÊÇ***£¬valueÊÇ*** Èç¹ûÄã±à³Ì¾Ñé²»·á¸»µÄ»°£¬×öÕâôһ²½»áºÜÓаïÖú |
4Â¥2014-11-11 08:39:55
À
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 621.5
- É¢½ð: 4
- Ìû×Ó: 30
- ÔÚÏß: 15.2Сʱ
- ³æºÅ: 2652634
- ×¢²á: 2013-09-13
- רҵ: Ô×ӺͷÖ×ÓÎïÀí
|
#!/usr/bin/perl -w use strict; use warnings; open FILE, '<','3'; open AFILE, '<','1'; open OUT,'>','result'; my @array1 = <FILE>; my @array2 = <AFILE>; my @b; foreach my $array1(@array1){ chomp($array1); my @a1=split /\s+/,$array1; my $cnt1=0; foreach my $array2(@array2) { chomp($array2); my @a2 = split /\s+/,$array2; if ($a1[0] eq $a2[0] and $a1[1] eq $a2[1]){ @b=sort{$a<=>$b}($a1[2],$a1[3]); $cnt1++ if ($a2[2]>= $b[0] and $a2[2] <= $b[1]); } }print OUT "@a1 $cnt1\n" } Ç×£¬ÎÒ×Ô¼º¸Äд±ä³öÀ´ÁË£¬²»¹ýÎҸоõÔËÐкÜÂý£¬Äܲ»ÄܰïÎÒÐÞ¸ÄһϠ|
5Â¥2014-11-11 08:54:19
À
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 621.5
- É¢½ð: 4
- Ìû×Ó: 30
- ÔÚÏß: 15.2Сʱ
- ³æºÅ: 2652634
- ×¢²á: 2013-09-13
- רҵ: Ô×ӺͷÖ×ÓÎïÀí
|
my %hash; my $n=0; open OUTFILE,'>','result1'; open FILE,'<','1'; while (<FILE> { chomp; ($a,$b,$c)=(split)[0,1,2]; $hash{$a}=1; $hash{$b}=2; open FILE,'<','3'; while (<FILE>{ chomp; ($d,$e,$f,$g)=(split)[0,1,2,3]; if(exists $hash{$d} && exists $hash{$e})#ÕâÒ»ÐеÄÒâ˼ÊÇ£¬ÎÒÏëͬʱ´æÔÚÒ»ÎļþµÄµÚÒ»Áк͵ڶþÁÐ{ @array=sort{$a <=> $b}($f,$g);#ÒòΪÓÐÕý¸ºÁ´µÄÔÒò£¬ËùÒÔÎÒ½«3ÎļþµÄµÚÈýÁк͵ÚËÄÁнøÐÐÁËÒ»ÏÂÖØÐÂÅÅÐò next if $c<$array[0] ; next if $c>$array[1];#½«²»ÔÚÕâ¸öÇø¼äµÄËùÓÐ$c¸ø³ýÈ¥ $n++; #ÓжàÉÙ¸ö$c¾Í$n¾ÍΪ¶àÉÙ }print OUTFILE"$d $e $f $g $n\n"; } } |
6Â¥2014-11-11 09:03:03
wizardfan
ÖÁ×ðľ³æ (ÖøÃûдÊÖ)
- BioEPI: 18
- Ó¦Öú: 599 (²©Ê¿)
- ¹ó±ö: 1.818
- ½ð±Ò: 24632.2
- É¢½ð: 197
- ºì»¨: 48
- ɳ·¢: 2
- Ìû×Ó: 2254
- ÔÚÏß: 400.7Сʱ
- ³æºÅ: 1879241
- ×¢²á: 2012-07-05
- ÐÔ±ð: GG
- רҵ: ÉúÎïÐÅϢѧ
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
|
Ϊʲô²»ÄÜÌýÎҵĽ¨Ò飬ÏÈдËã·¨£¬ÔÙд³ÌÐòÄØ£¿£¡ Ëã·¨1 Õë¶ÔÎļþ1ÖеÄÿһÐУ¬¶ÁÒ»±éÎļþ3£¬¿´ÔÚ²»ÔÚÇø¼äÄÚ Ëã·¨2 Õë¶ÔÎļþ1ÖеÄÿһÐУ¬¸ù¾ÝȾɫÌå·Ö×é Õë¶ÔÎļþ3ÖеÄÿһÐУ¬¸ù¾ÝȾɫÌå·Ö×é Õë¶ÔÎļþ1ÖеÄÿһÌõȾɫÌ壬ÕÒµ½¶ÔÓ¦µÄÎļþ3ÖеÄȾɫÌå·Ö×飬Ȼºó½øÐзÖÎö Äã¿´Ï£¬Á½¸öËã·¨ÄĸöËã·¨ËÙ¶È»á¸ü¿ì£¿ |
7Â¥2014-11-11 09:34:48
À
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 621.5
- É¢½ð: 4
- Ìû×Ó: 30
- ÔÚÏß: 15.2Сʱ
- ³æºÅ: 2652634
- ×¢²á: 2013-09-13
- רҵ: Ô×ӺͷÖ×ÓÎïÀí
8Â¥2014-11-25 09:43:38
À
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 621.5
- É¢½ð: 4
- Ìû×Ó: 30
- ÔÚÏß: 15.2Сʱ
- ³æºÅ: 2652634
- ×¢²á: 2013-09-13
- רҵ: Ô×ӺͷÖ×ÓÎïÀí
9Â¥2014-11-25 09:44:01













{
»Ø¸´´ËÂ¥