function [] = wilcoxonSRT(one, two)

warning ("off", "all");
significance = 0;

x = load(one);
y = load(two);
# Define Original X1 - Y1
z = x - y;
# Define Absolute |X1 - Y1|
zAbs = abs(z);
# Perform remaining operations
[zAbs, indices] = sort(zAbs);
z = z(indices);
sizeZ = size(zAbs);
rankings = [];
rCount = 1;
for i = 1 : sizeZ
	if (zAbs(i) > 0)
			rankings(1,rCount) = zAbs(i);
			if(z(i) > 0)
				rankings(2,rCount) = 1;
			end
			if(z(i) < 0)
				rankings(2,rCount) = -1;
			end
			rCount = rCount + 1;
	end
end
if (size(rankings) > 1)
	[rankings(1,:)] = ranks(rankings(1,:));
	# Calculate W value	
	[junk, rSize] = size(rankings);
	w = 0;
	for i = 1 : rSize
		w = w + (rankings(2,i) * rankings(1,i));
	end
	# calculate Ow
	oW = sqrt((rSize * (rSize + 1) * ((2 * rSize)+1))/6);
	zScore = (w - 0.5) / oW;
	# Test significance at 0.95 level
	if (zScore > 1.645 & mean(x) > mean(y))
		significance = 1;
	end
	if (zScore < 1.645 & mean(x) < mean(y))
		significance = -1;
	end
end
significance