Dear Monjes!!

I'm trying to generate a couple of things here. First I wrote a program that allows me to determine the generations that is needed for two lineages to meet back in time (evolving lineages). I'm assuming a mean of 4N/(k(k-1), with k= 2 and N=10000. My program is working and this is the script:

$k = 25; 
$N = 10000;

while ($k>=2)
    {
    if ($k==1)     
        {last;}
        $mean = (4*$N)/($k*($k-1));
        $time = (-log(rand)*$mean);
        
        push (@cls, $time);
        
        $k=$k-1;
        }
    open (GENE,'>coal.xls'); 
    for ($time=0; $time<24; $time ++) {
        print GENE "@cls[$time]\n";
        }

[download]

Now I need to write a simulation using an exponential distribution built from 5000 numbers, and afterwards get the mean, and the number of estimates above and below the mean. For this I need to use the same program as above. Although it is not working. So far this is my non-working scrypt:

#!usr/bin/perl
use warnings;

for ($sample=0; $sample<5000; $sample ++)
    {
    $k = 25; 
    $N = 10000;
    @cls=();
    $total=0;
    
while ($k>=2)
    {
    if ($k==1)     
        {last;}
        $mean = (4*$N)/($k*($k-1));
        $time = (-log(rand)*$mean);
        
        push (@cls, $time);
        $k=$k-1;
    }
        
        $ntotal+=$_ for @first;
        $nmean = $ntotal/5000;
        print "$nmean\n";
        foreach $first(@first)
        {
        if ($first<=$nmean) {push(@low, $first)}}
            $low= scalar @low;
            $hig= scalar (5000-@low);
                
            print "$low\n";
            print "$hig\n";
                
                
                }

[download]
Thanks in advance!

In reply to help finding mean, numbers above and below mean from an exponential distribution by kaupifalco

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post, it's "PerlMonks-approved HTML":


  • 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, details, 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, summary, 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:
    & &amp;
    < &lt;
    > &gt;
    [ &#91;
    ] &#93;
  • Link using PerlMonks shortcuts! What shortcuts can I use for linking?
  • See Writeup Formatting Tips and other pages linked from there for more info.