sub correlation
{
  . . .
  my( $x_ref, $y_ref ) = @_;     
  my $lfit = Statistics::LineFit->new();
  $lfit->setData( $x_ref, $y_ref );
  . . .
}