sub correlation {

  print "I am in the correlation subroutine";

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

  print "\nWe are in the correlation subroutine\n\n";


 # print Dumper ($x_ref, $y_ref);

  print Dumper length($x_ref);
  print Dumper length($y_ref);


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

  my ($intercept, $slope) = $lfit->coefficients();
  my $rSquared = $lfit->rSquared();