# concats current_file1 = open("500snp.summary.txt", 'w') current_file2 = open("500snp.simple.txt", 'w') wininfo={} inxs=open("window.key.txt") for line_idx, line in enumerate(inxs): cols = line.replace('\n', '').split('\t') # 1 scaffold_1 114 scaffold_1 104817 # 2 scaffold_1 104817 scaffold_1 192886 wininfo[cols[0]]=cols[1]+"."+cols[2]+"."+cols[3]+"."+cols[4] names={} good_bad={} nos=[0,0] src = open("samples.187.txt", "rU") # for line_idx, line in enumerate(src): cols = line.replace('\n', '').split('\t') # 9 109 0 in1=open(cols[1]+".minors.txt","rU") ccx=[0,0] for line_idx, line in enumerate(in1): colx = line.replace('\n', '').split('\t') # 1 1952 8 0.00408163265306 # 2 1941 4 0.00205655526992 ccx[0]+=int(colx[1]) ccx[1]+=int(colx[2]) current_file1.write(cols[1]+'\t'+wininfo[colx[0]]+'\t'+colx[1]+'\t'+colx[2]+'\t'+colx[3]+'\n') current_file2.write(colx[3]+'\n') nos[0]+=float(colx[1]) nos[1]+=float(colx[2]) print cols[1],ccx[0]+ccx[1],float(ccx[1])/float(ccx[0]+ccx[1]) print "Grand Mean frac rare ",nos[1]/(nos[1]+nos[0])