use warnings; use strict; my $runs = 1; # for testing code # Program vars my $chr_X_input = "bootstrap_data.txt"; my $range = 1552; # total number of array elements my $pi_sum; my $L_sum; my $differences_sum; my $coverage_sum; open my $CHR_X_INPUT, '<', $chr_X_input or die "Can't open chromosome X input: $!"; chomp( my @X_info = <$CHR_X_INPUT> ); close $CHR_X_INPUT; for ( 1 .. $runs ) { for ( 1 .. 26 ) { my ( $pi, $L, $differences, $coverage ) = split /\t/, $X_info[ int( rand($range) ) ]; $pi_sum += $pi; $L_sum += $L; $differences_sum += $differences; $coverage_sum += $coverage; my $PAR1_diversity = ( ( $pi_sum / $L_sum ) / ( $differences_sum / $coverage_sum ) ); } }