免费注册 查看新帖 |

Chinaunix

  平台 论坛 博客 文库
最近访问板块 发新帖
查看: 1700 | 回复: 5
打印 上一主题 下一主题

求助:这段程序哪里出了问题呢? [复制链接]

论坛徽章:
0
跳转到指定楼层
1 [收藏(0)] [报告]
发表于 2009-12-15 10:58 |只看该作者 |倒序浏览
DATA
相关基因.txt
HLA-DRB1        3123
NFKBIL1        4795
PADI4        23569
SLC22A4        6583
RUNX1        861
STAT4        6775
PTPN22        26191
TRAF1        7185
C5        727   
IRF5        3663
CD40        958  
pathway__all_gene.txt
hsa00010 GENE 2645 3098 3099 3101 2821 5211 5213 5214 2203 8789 226 229 230 7167 729708 2597 5230 5232 441531 5223 5224 728188 2023 2026 2027 5313 5315 5160 5161 5162 1737 1738 160287 3939 3945 3948 92483 124 125 126 127 128 130 131 10327 217 218 219 223 224 501 220 221 222 55902 84532 130589 5236 55276 2538 57818 669 5105 5106
hsa00020 GENE 1431 47 48 50 3417 3418 3419 3420 3421 4967 55753 1743 1738 8802 283398 8801 8803 6389 6390 6391 642502 6392 651820 2271 4190 4191 5091 5105 5106 5160 5161 5162 1737
hsa00030 GENE 2821 2539 25796 9563 5226 6120 7086 8277 84076 6888 22934 51071 64080 5236 55276 221823 5631 5634 226 229 230 5211 5213 5214 2203 8789
hsa00040 GENE 2990 10720 10941 54490 54575 54576 54577 54578 54579 54600 54657 54658 54659 7363 7364 7365 7366 7367 79799 7358 7360 6120 9942 231 51181
hsa00051 GENE 4351 5372 5373 29925 29926 2762 7264 8790 197258 2645 3098 3099 3101 5211 5213 5214 2203 8789 5207 5208 5209 5210 29085 8776 8898 9107 57016 80270 3795 6652 231 226 229 230 7167 729708
hsa00052 GENE 2584 2585 2582 7360 5236 55276 2645 3098 3099 3101 2538 57818 3938 2683 3906 8704 2592 2717 231 5211 5213 5214 8972 2548 2595
hsa00053 GENE 7358 10720 10941 54490 54575 54576 54577 54578 54579 54600 54657 54658 54659 7363 7364 7365 7366 7367 79799 55586 217 218 219 223 224 501
hsa00061 GENE 31 32 27349 2194 54995 55301
hsa00062 GENE 10449 3032 3028 3033 3295 3030 1892 51102 5538 9374
hsa00071 GENE 38 39 30 10449 3032 3028 3033 3295 3030 1962 1892 51 8310 35 34 33 36 37 2639 2180 2181 2182 23305 51703 126129 1374 1375 1376 10455 1632 1579 284541 124 125 126 127 128 130 131 217 218 219 223 224 501
hsa00072 GENE 3157 3158 3155 5019 64064 38 39 56898 622
hsa00100 GENE 3156 4598 10654 4597 3422 91734 2224 9453 2222 6713 4047 1595 7108 6307 50814 51478 10682 1718 6309 1717 1594 1728 2677 79001
hsa00120 GENE 3988 1056 6646 8435 1581 6715 6716 6718 1109 1593 124 125 126 127 128 130 131 217 218 219 223 224 501 10998 28976 57016 80270 30 10449 3032 570
hsa00130 GENE 27235 51805 51004 84274 10229
hsa00140 GENE 1583 3283 3284 1589 1586 1584 1585 3290 3291 6718 1109
hsa00150 GENE 3283 3284 1588 3292 3294 51144 51478 7923 10720 10941 54490 54575 54576 54577 54578 54579 54600 54657 54658 54659 7363 7364 7365 7366 7367 79799 1584 1585 3290 3291 6715 6716 1109 6718 6820 412 6783 414 415 3293 114049 131965 51409 51451 55798 60487 9836



#!/usr/bin/perl
open(in1,"E:\\相关基因.txt");
open(in2,"E:\\pathway__all_gene.txt");
open(stdout,">E:\\RA_result1.txt");
while($line1 = <in1>)
{
    chomp $line1;
    @array1 = split(/\t/,$line1);
    $y{$array1[1]} = $array1[0];
}
while($line2 = <in2>)
{
    chomp $line2;
    @array2 = split(/ /,$line2);
    $count = @array2;
    for($i=2;$i<$count;$i++)
    {
      if(exists($y{$array2[$i]}) == 1)
      {
         print "$array2[$i]\t$array2[0]\n";
      }
    }
}

结果存在问题,看了好久都没发现问题出在哪?高手指点一下啊
我只列了部分的数据,期望得到相关基因出现在哪一个通路的结果,可发现对于有的基因在实际通路中存在的,结果中却没有列出来,问题出在哪呢?谢谢大家的指点了。。

[ 本帖最后由 yn521yn 于 2009-12-15 16:05 编辑 ]

论坛徽章:
0
2 [报告]
发表于 2009-12-15 11:00 |只看该作者
flw曰:不会缩进的程序员不是好程序员。
LZ把code排下版吧~

论坛徽章:
78
双子座
日期:2013-10-15 08:50:09天秤座
日期:2013-10-16 18:02:08白羊座
日期:2013-10-18 13:35:33天蝎座
日期:2013-10-18 13:37:06狮子座
日期:2013-10-18 13:40:31双子座
日期:2013-10-22 13:58:42戌狗
日期:2013-10-22 18:50:04CU十二周年纪念徽章
日期:2013-10-24 15:41:34巨蟹座
日期:2013-10-24 17:14:56处女座
日期:2013-10-24 17:15:30双子座
日期:2013-10-25 13:49:39午马
日期:2013-10-28 15:02:15
3 [报告]
发表于 2009-12-15 11:15 |只看该作者

回复 #1 yn521yn 的帖子

什么问题呢?描述下吧

论坛徽章:
0
4 [报告]
发表于 2009-12-15 11:22 |只看该作者

回复 #1 yn521yn 的帖子

换行符或者分割符问题吧。楼主最好查下文件格式。打印下两个变量的内容是否正确。

[ 本帖最后由 toniz 于 2009-12-15 11:25 编辑 ]

论坛徽章:
0
5 [报告]
发表于 2009-12-15 11:22 |只看该作者
至少也要说说  结果哪里出问题了啊  要什么样的结果  
否则让大家在那里猜啊?

论坛徽章:
0
6 [报告]
发表于 2009-12-15 14:03 |只看该作者
代码最好排版一下···好看点····

楼主也没说是虾米问题··不好解决··
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

  

北京盛拓优讯信息技术有限公司. 版权所有 京ICP备16024965号-6 北京市公安局海淀分局网监中心备案编号:11010802020122 niuxiaotong@pcpop.com 17352615567
未成年举报专区
中国互联网协会会员  联系我们:huangweiwei@itpub.net
感谢所有关心和支持过ChinaUnix的朋友们 转载本站内容请注明原作者名及出处

清除 Cookies - ChinaUnix - Archiver - WAP - TOP