sub correlation {


  my ($x_ref, $y_ref) = @_;   
  

  my $lfit = Statistics::LineFit->new();
  $lfit->setData($x_ref, $y_ref);