24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1273  |  回复: 10

junruiwang

新虫 (小有名气)

[求助] 求各位大神,perl 对特定范围内的数据求平均值 已有1人参与

文件为:
Chr        Pos        Y
1        49820        -1.68855
1        54414        -1.55555
1        99585        -1.05338
1        103929        -0.523771
1        124119        -0.261121
1        200009        -0.1234
1        390189        2.2456
2        48625        0.0250024
2        88663        -0.659643
2        90128        -0.302134
2        122780        -0.366919
2        555583        -1.00859
想根据第一列和第二列的范围来求第三列的平均值,具体为:第二列的按照每100,000的长度来分组,把第三列求一个平均值,如果有些范围没有数值则写为0(如下)。想得到的结果是:
1  100,000    -1.432493333
1  200,000    -0.392446
1  300,000    -0.1234
1  400,000    2.2456
2  100,000    -0.3259234
2  200,000    0
2  300,000    0
2  400,000    0
2  500,000    0
2  600,000    -1.00859
其中100,000、200,000...代表0-100,000、100,000-200,000范围的意思。
求各位大神帮忙啦!
回复此楼

» 收录本帖的淘帖专辑推荐

程序

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

[perl code]
# Program in Perl
use Data:umper;
use POSIX;
my @sumY, @numY;
open my $fh, '<', "output.txt" or
     die "could not open file: $!";
$line = <$fh>;
while (my $line = <$fh> {
    my($Chr, $Pos, $Y) =  split /\s+/, $line;
    $cPos = floor($Pos / 100000) ;
    $sumY[$Chr][$cPos]+=$Y;
    $numY[$Chr][$cPos]++;
}
for (my $j= 1; $j <= 2; $j++)
{
    for (my $i= 0; $i <= 5; $i++)
    {
        $range = ($i+1)*100000;
        if ( !defined $numY[$j][$i] )
           { $avg = 0; }
        else
           { $avg = $sumY[$j][$i] / $numY[$j][$i]; }
        print "$j $range $avg\n";
    }
}
close $fh;

[output]
1 100000 -1.43249333333333
1 200000 -0.392446
1 300000 -0.1234
1 400000 2.2456
1 500000 0
1 600000 0
2 100000 -0.3122582
2 200000 -0.366919
2 300000 0
2 400000 0
2 500000 0
2 600000 -1.00859
8楼2016-08-18 06:01:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

FMStation

至尊木虫 (知名作家)

看不懂

» 本帖已获得的红花(最新10朵)

2楼2016-08-17 17:52:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
#!/usr/bin/perl
use warnings;
use strict;

my $sum = 0;
my $n = 0;
while (<> {
    $sum += $_;
    $n++;
}
print $sum/$n, "\n";

https://stackoverflow.com/questi ... averages-using-perl
3楼2016-08-17 18:09:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

CODE:
#!/usr/bin/perl
use warnings;
use strict;

my $sum = 0;
my $n = 0;
while (<>) {
    $sum += $_;
    $n++;
}
print $sum/$n, "\n";

https://stackoverflow.com/questions/6906116/compute-averages-using-perl
4楼2016-08-17 18:10:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

James_ever

金虫 (文坛精英)

必须用这软件??
Take few risks early, the bonus following!科研MSN群:FutureChina@groups.live.com
5楼2016-08-17 19:07:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

junruiwang

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by FMStation at 2016-08-17 17:52:51
看不懂

就是按照0-100000,100001-200000,200001-300000......这些范围(每100000长度)对第三列进行求平均值。每个范围求一个平均值输出哦
6楼2016-08-17 19:20:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

junruiwang

新虫 (小有名气)

引用回帖:
5楼: Originally posted by James_ever at 2016-08-17 19:07:00
必须用这软件??

或者用别的语言写也成,只要能实现就行哈
7楼2016-08-17 19:23:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
junruiwang: 金币+100, ★★★★★最佳答案, 感谢~ 2016-08-20 14:33:13
CODE:
~~~~~~code~~~~~~~~~
# Program in Perl
use Data::Dumper;
use POSIX;
my @sumY, @numY;
open my $fh, '<', "output.txt" or
     die "could not open file: $!";
$line = <$fh>;
while (my $line = <$fh>) {
    my($Chr, $Pos, $Y) =  split /\s+/, $line;
    $cPos = floor($Pos / 100000) ;
    $sumY[$Chr][$cPos]+=$Y;
    $numY[$Chr][$cPos]++;
}
for (my $j= 1; $j <= 2; $j++)
{
    for (my $i= 0; $i <= 5; $i++)
    {
        $range = ($i+1)*100000;
        if ( !defined $numY[$j][$i] )
           { $avg = 0; }
        else
           { $avg = $sumY[$j][$i] / $numY[$j][$i]; }
        print "$j $range $avg\n";
    }
}
close $fh;

~~~~Output~~~~~~
1 100000 -1.43249333333333
1 200000 -0.392446
1 300000 -0.1234
1 400000 2.2456
1 500000 0
1 600000 0
2 100000 -0.3122582
2 200000 -0.366919
2 300000 0
2 400000 0
2 500000 0
2 600000 -1.00859

» 本帖已获得的红花(最新10朵)

9楼2016-08-18 06:02:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

junruiwang

新虫 (小有名气)

引用回帖:
9楼: Originally posted by FMStation at 2016-08-18 06:02:35
~~~~~~code~~~~~~~~~
# Program in Perl
use Data:umper;
use POSIX;
my @sumY, @numY;
open my $fh, '<', "output.txt" or
     die "could not open file: $!";
$line = < ...

太感谢啦!关于自学perl,大神能否给指条明路啊?实验室只有我一个人刚开始学,感觉好费劲。。。

发自小木虫Android客户端
10楼2016-08-18 11:48:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 junruiwang 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南昌大学324求调剂 +6 hanamiko 2026-03-29 6/300 2026-03-31 16:35 by hypershenger
[考研] 一志愿a区211,085601-307分求调剂 +9 党嘉豪 2026-03-31 18/900 2026-03-31 16:34 by sophie2180
[基金申请] 面上5B能上会吗? +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 一志愿华东师范大学有机化学专业,初试351分,复试被刷求调剂! +8 真名有冰 2026-03-29 9/450 2026-03-31 15:40 by whf1979
[考研] 085601 材料工程 313分 求调剂 +6 Ong3 2026-03-27 6/300 2026-03-31 15:18 by yanflower7133
[考研] 086000调剂 +4 7901117076 2026-03-26 4/200 2026-03-31 14:30 by 记事本2026
[考研] 生物学 296 求调剂 +7 朵朵- 2026-03-26 9/450 2026-03-31 14:26 by jp9609
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +6 星沉uu 2026-03-26 7/350 2026-03-31 10:19 by GdShizy
[考研] 276求调剂 +3 赵久华 2026-03-29 3/150 2026-03-31 10:06 by cal0306
[考研] 083000学硕274求调剂 +12 Li李鱼 2026-03-26 12/600 2026-03-31 10:01 by cal0306
[考研] 英一数一408,总分284,二战真诚求调剂 +3 12.27 2026-03-30 5/250 2026-03-31 00:47 by 不吃芒果17
[考研] 279求调剂 +12 j的立方 2026-03-29 12/600 2026-03-30 20:30 by dick_runner
[考研] 生物学学硕,一志愿湖南大学,初试成绩338 +7 YYYYYNNNNN 2026-03-26 9/450 2026-03-30 20:29 by YYYYYNNNNN
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
[考研] 332求调剂 +6 Lyy930824@ 2026-03-29 6/300 2026-03-30 10:53 by 1172367218
[考研] 求调剂 +10 张zz111 2026-03-27 11/550 2026-03-30 09:17 by 无际的草原
[考研] 0856材料化工调剂 总分330 +14 zhubinhao 2026-03-27 14/700 2026-03-29 10:01 by Sjndkwm
[硕博家园] 招收生物学/细胞生物学调剂 +4 IceGuo 2026-03-26 5/250 2026-03-29 01:25 by griffith2014
[考研] 调剂推荐 +5 清酒714 2026-03-26 6/300 2026-03-27 11:12 by 不吃魚的貓
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
信息提示
请填处理意见