24小时热门版块排行榜    

查看: 1262  |  回复: 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的回帖

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的回帖
查看全部 11 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 312求调剂 +5 上岸吧ZJY 2026-03-22 7/350 2026-03-25 22:20 by 544594351
[考研] 材料求调剂 +4 .m.. 2026-03-25 4/200 2026-03-25 21:30 by peike
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 321求调剂 +3 璞玉~~ 2026-03-25 3/150 2026-03-25 19:07 by Zhanglab-TJU
[考研] 359求调剂 +3 王了个楠 2026-03-25 3/150 2026-03-25 12:50 by Dyhoer
[考研] 081700 调剂 267分 +10 迷人的哈哈 2026-03-23 10/500 2026-03-25 12:12 by userper
[考研] 0854电子信息求调剂 324 +4 Promise-jyl 2026-03-23 4/200 2026-03-25 11:36 by Sugarlight
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 11/550 2026-03-25 10:23 by userper
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 291求调剂 +8 hhhhxn.. 2026-03-23 8/400 2026-03-23 23:15 by peike
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
信息提示
请填处理意见