my $lfit = Statistics::LineFit->new();
$lfit->setData($x_ref, $y_ref);
  
my ($intercept, $slope) = $lfit->coefficients();
my $rSquared = $lfit->rSquared();
my $pValue = $lfit->pValue();  #  I need a line like this