应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Shell&Perl » perl程序求助~~
21/1返回列表
查看: 429  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

philip890220

铜虫 (初入文坛)

[求助] perl程序求助~~

正如附件及图片所示，注意下一啊，gene那一列中的gene是有重复的，我的目的就是将重复出现的基因每列分别取平均值，最后输出一个文件！
就拿BRO1来说，他出现两次，第二列中的两个数字分别是0.0375155和0.0691334，九曲他俩的平均值，第三列亦然~~，最后的输出结果应该是这样的：
GENE        GSM630089        GSM630090
BRO1               
COX11               
MET16               
SGO1               
YOR121C               
数字我省略了，反正就是有重复的取平均值，不重复的直接输出！，求助~~~       
谢谢喽~~~PS:实际处理的文件较大，附件中仅是例子~~~

1.gif

[ Last edited by philip890220 on 2012-9-25 at 16:48 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : test.txt
    • 2012-09-25 16:32:22, 256 bytes

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-09-25 16:46:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jackie1179

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
philip890220: 金币+10, ★★★很有帮助, good~~ 2012-09-29 00:34:55
这个应该可以解决你的问题，输出文件为output.txt
CODE:
#!/usr/bin/perl

use warnings;
use strict;

my $in = "test.txt";

open IN,"<",$in || die "Cannot open this file!$!";
my %array;

;
while() {
  chomp;
  if(/^(\S+)\s+(\S+)\s+(\S+)/) {
    push @{$array{$1}{'col2'}},$2;
    push @{$array{$1}{'col3'}},$3;
  }
}

close IN;

foreach my $key(sort keys %array) {
  my $cnt = scalar(@{$array{$key}{'col2'}});
  foreach (@{$array{$key}{'col2'}}) {
    $array{$key}{'ave2'} += $_;     
  }
  
  $array{$key}{'ave2'} /= $cnt;

  foreach (@{$array{$key}{'col3'}}) {
    $array{$key}{'ave3'} += $_;
  }
  $array{$key}{'ave3'} /= $cnt;

}

my $out = "output.txt";
open OUT,">",$out || die "Cannot write to this file!$!";
foreach my $key(sort keys %array) {
  print OUT "$key\t$array{$key}{'ave2'}\t$array{$key}{'ave3'}\n";

}

close OUT;

exit(0);
赞一下(2人) 回复此楼
2楼2012-09-27 10:00:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 philip890220 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +9 majunge000 2026-08-11 11/550 2026-08-15 23:20 by botar
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 68/3400 2026-08-15 23:19 by 医学老男孩
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭