- 论坛徽章:
- 10
|
root@lp:~/lp/python/data# cat x.csv
chr1 100 300
chr3 60 110
chr4 40 200
chr5 25 90
root@lp:~/lp/python/data# cat get_data.py
- import pandas as pd
- a = pd.read_csv('a.txt',header=None,names=['id','start','end'],sep='\s+')
- b = pd.read_csv('b.txt',header=None,names=['id','start','end'],sep='\s+')
- x = []
- for i in range(b.shape[0]):
- bstart = b.iloc[i].start
- bend = b.iloc[i].end
- bid = b.iloc[i].id
- xx = a[a.id==bid]
- flag = False
- for j in range(xx.shape[0]):
- if xx.iloc[j].end >= bend and xx.iloc[j].start <= bstart:
- flag = True
- if flag == True:
- x.append(False)
- else:
- x.append(True)
- b[x].to_csv('x.csv',sep=' ',index=None,header=None)
复制代码 |
|