import scipy.stats
f = open('beauty.csv', 'r')
head = f.readline()
data = f.readlines()
f.close()
for k in xrange(2, 8) :
wage = []
x = []
for row in data :
p = row.split(';')
wage.append(float(p[0]))
x.append(int(p[k]))
wage, x = vector(wage), vector(x)
rb = scipy.stats.pointbiserialr(x, wage)
print k, rb.correlation, rb.pvalue