#!/usr/bin/perl

use strict;
use warnings;
use Data::Dump qw(dump);

my $ifile = "facehbi_dti_np_models_all.txt";
my $thresh = 0.3;
my %model;
open IDF, "<$ifile" or die "No such file\n";
while (<IDF>){
        if (/-------/ && $model{"r2"}>0.3) {
                print $model{"ni_var"}.",".$model{np_var}."\nr2 = ".$model{"r2"}.", p-value = ".$model{"pvalue"}."\npv_".$model{"ni_var"}." = ".$model{"pv_ni"}.", pv_SUVR = ".$model{pv_suvr}." \n";
        };
        if (/^NP:.*/) {($model{"np_var"}, $model{"ni_var"}) = /^NP:\s+(\w+)\s+NI:\s+(\w+)\s*$/};
        if (/^R2:.*/) {($model{"r2"}, $model{"pvalue"}) = /^R2:\s+(\S+)\s+p-value:\s+(\S+)/};
        if (/^p-value.*/) {($model{"pv_ni"}, $model{"pv_suvr"}) = /p-value.*:\s+(\S+)\s+p-value.*:\s+(\S+)/}
}
close IDF;