前兩列是位點信息,希望統(tǒng)計同一個位點下,第三列x有多少個,X有多少個,在python中應(yīng)該用什么方法能比較快的實現(xiàn)呢
Chr1 1001 x
Chr1 1001 X
Chr1 1001 x
Chr1 10018 X
Chr1 1002 x
Chr1 1002 X
目前我想的是根據(jù)前兩列為鍵,最后一列是X和x的頻率為值,建立兩個字典,然后分別輸出成兩個文件,怎么輸出成一個列表類似:
chr loc x X
Chr1 1001 2 1
chr1 1002 1 1
chr1 10018 0 1
思考了一下,好像可以通過列表解決,tlist統(tǒng)計所有行,mlist統(tǒng)計只有X的行,uniqlist統(tǒng)計不重復(fù)行
通過遍歷uniqlist統(tǒng)計X行和x+X行的個數(shù)
tlist=[]
mlist=[]
uniqlist=[]
for l in f.readlines():
line = l.strip().split("\t")
loc="\t".join(line[0:2])
tlist.append(loc)
uniqlist=list(set(tlist))
if line[2]=="X":
mlist.append(loc)
for i in uniqlist:
out.write(i+"\t"+str(uniqlist.count(i))+"\t"+str(mlist.count(i))+"\n")北大青鳥APTECH成立于1999年。依托北京大學優(yōu)質(zhì)雄厚的教育資源和背景,秉承“教育改變生活”的發(fā)展理念,致力于培養(yǎng)中國IT技能型緊缺人才,是大數(shù)據(jù)專業(yè)的國家
達內(nèi)教育集團成立于2002年,是一家由留學海歸創(chuàng)辦的高端職業(yè)教育培訓(xùn)機構(gòu),是中國一站式人才培養(yǎng)平臺、一站式人才輸送平臺。2014年4月3日在美國成功上市,融資1
北大課工場是北京大學校辦產(chǎn)業(yè)為響應(yīng)國家深化產(chǎn)教融合/校企合作的政策,積極推進“中國制造2025”,實現(xiàn)中華民族偉大復(fù)興的升級產(chǎn)業(yè)鏈。利用北京大學優(yōu)質(zhì)教育資源及背
博為峰,中國職業(yè)人才培訓(xùn)領(lǐng)域的先行者
曾工作于聯(lián)想擔任系統(tǒng)開發(fā)工程師,曾在博彥科技股份有限公司擔任項目經(jīng)理從事移動互聯(lián)網(wǎng)管理及研發(fā)工作,曾創(chuàng)辦藍懿科技有限責任公司從事總經(jīng)理職務(wù)負責iOS教學及管理工作。
浪潮集團項目經(jīng)理。精通Java與.NET 技術(shù), 熟練的跨平臺面向?qū)ο箝_發(fā)經(jīng)驗,技術(shù)功底深厚。 授課風格 授課風格清新自然、條理清晰、主次分明、重點難點突出、引人入勝。
精通HTML5和CSS3;Javascript及主流js庫,具有快速界面開發(fā)的能力,對瀏覽器兼容性、前端性能優(yōu)化等有深入理解。精通網(wǎng)頁制作和網(wǎng)頁游戲開發(fā)。
具有10 年的Java 企業(yè)應(yīng)用開發(fā)經(jīng)驗。曾經(jīng)歷任德國Software AG 技術(shù)顧問,美國Dachieve 系統(tǒng)架構(gòu)師,美國AngelEngineers Inc. 系統(tǒng)架構(gòu)師。