import numpy as np
from scipy.stats.stats import pearsonr
# from weavile import xmat
ifile = open('/home/scottn/fourmercountantisense.txt')
xmat = np.genfromtxt(ifile, delimiter='\t', dtype=None, names=None)
n = len(xmat)
#n = 50
output = np.zeros(shape=(n,n))
ofile = open('/home/scottn/pythoutfourmerantisense.txt', 'w+')
for i in range(0, n):
	output[i][i] = 1
	for j in range(i+1, n):
		c = pearsonr(xmat[i], xmat[j])
	#	print(abs(c[0]))
	#	if c[0] < 0:
	#		output[i][j] = 0
	#	else:
		output[i][j] = round(c[0], 5)
		output[j][i] = output[i][j]
for row in output:
	ofile.write('\t'.join([str(x) for x in row])+'\n')
ofile.close()