BEGIN { n=0; } { if (group==1) { #group 634 x[n]=$2*1.0; y[n]=-$4*1.0; } else { # x[n]=$3*1.0; y[n]=-$4*1.0; } n++; } END { P=0; for (i=0; i x[j] && y[i] > y[j])) P++; if (x[i] == x[j] || y[i] == y[j]) P += 0.5; } } print ((4*P) / (n*(n-1))) - 1; } #sub kendall(\@\@) #{ my($x, $y) = @_; # my $P = 0; # # die "unequal number of x and y points when computing tau\n" if @$x != @$y; # # for(my $i = 0; $i < @$x; ++$i) # { for(my $j = $i + 1; $j < @$x; ++$j) # { ++$P if ($x->[$i] < $x->[$j] && $y->[$i] < $y->[$j]) # || ($x->[$i] > $x->[$j] && $y->[$i] > $y->[$j]); # $P += 0.5 if $x->[$i] == $x->[$j] || $y->[$i] == $y->[$j]; # } # } # # (4 * $P) / ( @$x * (@$x - 1)) - 1; #} #