标签:文件 python RECList transList 行与行 int rho SNPPosList line
在处理文本文件时,很多时候需要我们处理跨行的数据,但是用for循环处理不是很方便,想了一个歪招来处理不是很大的数据。
核心思想就是将上一行的东西存在一个列表里,到下一行用完这个数据在循环体里将列表初始化,再将这行数据存到列表里,这样循环就可以把两行数据联系起来。
另外这样做的还有另外一个好处,就是可以初始化第一行,我们可以通过判断这个“传递列表”空与否来区分标题行的输入和其他循环体的输入情况。
from __future__ import division from sys import argv from itertools import islice Mutation = open('DivFile_sort_chr' + argv[1] +'.txt',"r") REC = open('GeneticMap_Chr' + argv[1],"r") RECList = open('REC_Chr_' + argv[1], "a") transList = [] SNPPosList = [] rholist = [] for lines in islice(Mutation, 1, None): line = lines.strip('\n').split("\t") POS = line[0] SNPPosList.append(POS) print "SNPPos done" for lines in islice(REC, 1, None): line = lines.strip('\n').split("\t") if transList == []: transList.append((line[0],line[1])) RECList.write("position\trate\n") RECList.write('%s\t%s\n' %(line[0],0)) print "headline" else: rho = (float(line[1]) - float(transList[0][1])) / (float(line[0]) - float(transList[0][0])) for i in SNPPosList: if int(i) < int(line[0]) and int(i) > int(transList[0][0]): index = int(SNPPosList.index(i) - 1) site_rho = (int(i) - int(SNPPosList[index])) * rho RECList.write('%s\t%s\n' %(i,site_rho)) else: continue transList = [] transList.append((line[0],line[1],rho)) rholist.append(rho) print rholist for i in SNPPosList: if i > int(transList[0][0]): index = int(SNPPosList.index(i) - 1) site_rho = (int(i) - int(SNPPosList[index])) * rho RECList.write('%s\t%s\n' %(i,site_rho)) else: continue Mutation.close() REC.close() RECList.close()
标签:文件,python,RECList,transList,行与行,int,rho,SNPPosList,line 来源: https://www.cnblogs.com/muuyouzhi/p/10494093.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。