http://www.perlmonks.org?node_id=969707


in reply to Population of HoAoA based on file contents

So, if I understand correctly, you want your results to look something like:
$data{"HG00553"}[162][16287215] = "0 0" $data{"HG00553"}[162][16287226] = "0 0" ... $data{"HG00638"}[162][16287851] = "1 1" ...
First, I'd suggest you do a little debug work to understand what your code is doing. Specifically:
while (<$in_file>) { print "$_"; chomp; my @snp_bins; if (/^SAMPLE/) { my ( $placeholder, @coords ) = split /,/; foreach my $coord (@coords) { push @snp_bins, int( $coord / 100_000 ); } } else { my ( $id, @snps ) = split /,/; push @individuals, $id; foreach my $individual (@individuals) { print "working on $individual\n"; foreach my $snp (@snps) { my $snp_bin = shift @snp_bins; print "snp_bin $snp_bin\n"; $data{$individual}[ [ $snp_bin ] ] = $snp; } } print "snp_bins: @snp_bins\n"; } }
I think that you'll find some things which you didn't expect. Also, I think that you might be better off using HoHoH instead. I believe that doing:
$data{$individual}[162] = "something"
will result in:
$data{$individual}[0..162]
even though you are only concerned with the former. So it would probably be better to:
$data{$individual}{162}{16287215} = "0 0"