function [r,test]=corrsign(x,y,ndof)
test=0;
dd=cov(x,y)/std(x)/std(y);
r=dd(1,2);
if (~exist('ndof'))
   ndof=fix(sqrt(dof(x)*dof(y)));
end
t=abs(r)*sqrt(ndof)/sqrt(1-r^2);
alpha=0.025;
ts=tq(1-alpha/2,ndof);
test=tp(t,ndof);
