import sys sys.path.append('/home/people/tc/svn/tc_sandbox/misc/') import statistics l1 = [48.2,54.6,58.3,47.8,51.4,52.0,55.2,49.1,49.9,52.6,] l2 = [52.3,57.4,55.6,53.2,61.3,58.0,59.8,54.8,] mean1,mean2,stderr,p = statistics.twosamplettest(l1,l2,) print p
for i in range(5): l[i] += [[]] fd = open('35C/%s/energies.txt' %(MD),'r') lines = fd.readlines() fd.close() for line in lines: for i in range(4): l[i][-1] += [float(line.split()[i+1])] Sum = float(line.split()[1])+float(line.split()[2])+float(line.split()[3])+float(line.split()[4]) l[4][-1] += [Sum] if MD == 'MD_2vb1_Glu35Asp52_01' and len(l[0][-1]) == 500: break Sum = 0 ##l = [[8.8,8.4,7.9,8.7,9.1,9.6,],[9.9,9.0,11.1,9.6,8.7,10.4,9.5],] mean1,mean2,stderr,p = statistics.twosamplettest(l[0][0],l[0][1],verbose=False,) Sum += mean1-mean2 print 'Asp52', mean1-mean2,1.962*stderr mean1,mean2,stderr,p = statistics.twosamplettest(l[1][0],l[1][1],verbose=False,) Sum += mean1-mean2 print 'protein', mean1-mean2,1.962*stderr mean1,mean2,stderr,p = statistics.twosamplettest(l[2][0],l[2][1],verbose=False,) Sum += mean1-mean2 print 'chloride', mean1-mean2,1.962*stderr mean1,mean2,stderr,p = statistics.twosamplettest(l[3][0],l[3][1],verbose=False,) Sum += mean1-mean2 print 'water', mean1-mean2,1.962*stderr mean1,mean2,stderr,p = statistics.twosamplettest(l[4][0],l[4][1],verbose=False,) print 'sum', mean1-mean2,1.962*stderr
fd.close() for line in lines: for i in range(4): l[i][-1] += [float(line.split()[i + 1])] Sum = float(line.split()[1]) + float(line.split()[2]) + float( line.split()[3]) + float(line.split()[4]) l[4][-1] += [Sum] if MD == 'MD_2vb1_Glu35Asp52_01' and len(l[0][-1]) == 500: break Sum = 0 ##l = [[8.8,8.4,7.9,8.7,9.1,9.6,],[9.9,9.0,11.1,9.6,8.7,10.4,9.5],] mean1, mean2, stderr, p = statistics.twosamplettest( l[0][0], l[0][1], verbose=False, ) Sum += mean1 - mean2 print 'Asp52', mean1 - mean2, 1.962 * stderr mean1, mean2, stderr, p = statistics.twosamplettest( l[1][0], l[1][1], verbose=False, ) Sum += mean1 - mean2 print 'protein', mean1 - mean2, 1.962 * stderr mean1, mean2, stderr, p = statistics.twosamplettest( l[2][0], l[2][1], verbose=False,
l1 = [ 48.2, 54.6, 58.3, 47.8, 51.4, 52.0, 55.2, 49.1, 49.9, 52.6, ] l2 = [ 52.3, 57.4, 55.6, 53.2, 61.3, 58.0, 59.8, 54.8, ] mean1, mean2, stderr, p = statistics.twosamplettest( l1, l2, ) print p