24小时热门版块排行榜    

查看: 378  |  回复: 1

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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楼2012-09-27 10:00:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 philip890220 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿吉大化学322求调剂 +3 17501029541 2026-03-23 4/200 2026-03-23 23:47 by Txy@872106
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +4 cccchenso 2026-03-23 4/200 2026-03-23 23:00 by 徐ckkk
[考研] 一志愿211 初试270分 求调剂 +4 谷雨上岸 2026-03-23 5/250 2026-03-23 21:18 by 不惑可乐
[考研] 070300化学求调剂 +8 苑豆豆 2026-03-20 8/400 2026-03-23 20:57 by baobaoye
[考研] 一志愿重庆大学085700资源与环境,总分308求调剂 +7 墨墨漠 2026-03-23 8/400 2026-03-23 20:36 by Creta
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +5 葵梓卫队 2026-03-18 7/350 2026-03-23 16:26 by lingjue
[考研] 316求调剂 +7 梁茜雯 2026-03-19 7/350 2026-03-23 16:21 by lingjue
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
信息提示
请填处理意见