import numpy as np
# from scipy.stats.stats import pearsonr
# from weavile import xmat
ifile = open('/home/scottn/sixmerbin.csv')
xmat = np.genfromtxt(ifile, delimiter=',', dtype=None, names=None).transpose()
n = len(xmat)
# n = 50
output = np.zeros(shape=(n,n))
ofile = open('/home/scottn/pythoutsixmer.csv', 'w+')
for i in range(0, n):
	output[i][i] = 1
	for j in range(i+1, n):
		output[i][j] = round(np.mean(xmat[i]*xmat[j])/4, 3)
		output[j][i] = output[i][j]
for row in output:
	ofile.write(','.join([str(x) for x in row])+'\n')
ofile.close()