foreach my $dataset (sort keys %datasets) {

    my $aov = Statistics::ANOVA->new();

    $aov->load( "$dataset\1", \@{$datasets{$dataset}[0]} );
    $aov->add( "$dataset\2", \@{$datasets{$dataset}[1]} );

    $str = $aov->anova(independent => 1, parametric => 1, ordinal => 0);

    print $str->{_stat}->{p_value} . "\n";
}