Comment on

Thanks!!!! It works fine now and, as you said, substantially faster also. This is what I did,

foreach my $pair (@lpairs) {
    my $n = 0;
    my $z = 0;
    my $hl;
    my $pvalue = 0;
    my $pvt = 0;
    my $fk;
    foreach $fk (%ipairs) {
        if( exists $ipairs{$fk} &&  exists $ipairs{$fk}{$pair} &&  exi
+sts $ipairs{$fk}{$pair}{'pvalue'}){
        #if($ipairs{$fk}{$pair}{'pvalue'}){
            $pvt = $ipairs{$fk}{$pair}{'pvalue'};
            if($pvt){
                unless($hl){
                    $hl = $ipairs{$fk}{$pair}{'head'};
                }
                $n++;
                $z+= qnorm($ipairs{$fk}{$pair}{'pvalue'});
            }
        }
    }
    if($n>2){
        $z = $z/sqrt($n);
        $pvalue = pnorm($z);
    }    
    if ($pvalue) {
        #printf "$pair -> %.4f\n", $pvalue; 
        printf OF "$hl %.4f $n\n", $pvalue;
    }

}
[download]

In reply to Re^2: memory issues by asqwerty
in thread memory issues by asqwerty

Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!

Titles consisting of a single word are discouraged, and in most cases are disallowed outright.

Read Where should I post X? if you're not absolutely sure you're posting in the right place.

Please read these before you post! —

Posts may use any of the Perl Monks Approved HTML tags:

a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr

You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)

	For:		Use:
	&		`&`
	<		`<`
	>		`>`
	[		`[`
	]		`]`

Link using PerlMonks shortcuts! What shortcuts can I use for linking?

See Writeup Formatting Tips and other pages linked from there for more info.


go ahead... be a heretic
	PerlMonks