#!/usr/bin/perl -w
use strict;
use Statistics::Distributions;
$|=1;

open my $OUTPUT, ">", "mh_dist_outliers1d.tsv" or die $!;

my %PCA=();
open my $PRCOMPFILE, "<", "strain_PCA.csv" or die $!;
while (my $line=<$PRCOMPFILE>){
    chomp $line;
    my @line=split /,/, $line;
    $PCA{$line[0]}=$line[1];
}
close $PRCOMPFILE;

my %results=();
my @sizes=(1,25,50,100);
for my $size (@sizes){
    print "\n###Starting pop$size###\n";
    open my $POPLIST, "<", "pops.$size.tsv" or die $!;
    while (my $line=<$POPLIST>){
        print ".";
        chomp $line;
        my @line=split /\t/, $line;
        my $pop=shift @line;
        my @pos= map {$PCA{$_};} @line;