import numpy as np
from scipy.stats.stats import pearsonr
ifile = open('/home/scottn/sixmercountinantisense.txt')
xmat = np.genfromtxt(ifile, delimiter='\t', names=None, dtype=np.float64)
#from matriximport import xmat
n = len(xmat)
#n = 5
output = np.zeros(shape=(n,n))
ofile = open('/home/scottn/kinshipsixmercorr.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]))
		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()