0 
junctions = [2,9,15,20]

seq_1 = 'sauron'
seq_2 = 'corrupted'
seq_3 = 'numenor'
combined = 'sauroncorruptednumenor' #seq_1 + seq_2 + seq_3

count_1 = 1
count_2 = 1
count_3 = 2

我有一个包含 3 个字符串的列表（seq_1-3）。我将它们组合在一起以创建 1 个长字符串（组合）我有一个索引列表（连接点）。我为每个字符串设置了 3 个不同的计数器为零（count_1-3）

我要做的是在组合序列中找到每个结点 [2,9,15,20] 的位置。. . 如果是从 seq_1 --> count_1 += 1 如果是从 seq_2 --> count_2 += 1 从 seq_3 --> count_3 += 1

例子

junctions = [2,9,15,20]
count_1 = 0
count_2 = 0
count_3 = 0
combined = 'sauroncorruptednumenor'
seq_1 = 'sauron' #index 2 would be on 'u' in combined but originally from seq_1 so count_1 = count_1 + 1 
seq_2 = 'corrupted' #index 9 would be on 'r' in combined so count_2 += 1
seq_3 = 'numenor' #index 15 would be 'n' in combined so count_3 += 1, and 20 would be 'o' so count_3 += 1

让我知道我是否需要以不同的方式澄清

4

3 回答 3

1

你可以在这里使用collections.Counter和bisect.bisect_left：

>>> from collections import Counter
>>> import bisect
>>> junctions = [2,9,15,20]
>>> seq_1 = 'sauron'
>>> seq_2 = 'corrupted'
>>> seq_3 = 'numenor'
>>> lis  = [seq_1, seq_2, seq_3]

创建一个列表，其中包含每个seq_结尾处的索引：

>>> start = -1
>>> break_points = []
for item in lis:
    start += len(item) 
    break_points.append(start)
...     
>>> break_points
[5, 14, 21]

现在我们可以简单地循环并使用函数junctions查找列表中每个连接点的位置。break_pointsbisect.bisect_left

>>> Counter(bisect.bisect_left(break_points, jun)+1  for jun in junctions)
Counter({3: 2, 1: 1, 2: 1})

使用更好的输出collections.defaultdict：

>>> from collections import defaultdict
>>> dic = defaultdict(int)
for junc in junctions:
    ind = bisect.bisect_left(break_points, junc) +1
    dic['count_'+str(ind)] += 1
...     
>>> dic
defaultdict(<type 'int'>,
{'count_3': 2,
 'count_2': 1,
 'count_1': 1})

#accessing these counts
>>> dic['count_3']
2
于 2013-06-25T23:10:04.457 回答
1

你可以尝试一些基本的东西，比如

L_1 = len(seq_1)
L_2 = len(seq_2)
L_3 = len(seq_3)

junctions = [2, 9, 15, 20]
c_1, c_2, c_3 = (0, 0, 0)

for j in junctions:
    if j < L_1:
        c_1 += 1
    elif j < L_1 + L_2:
        c_2 += 1
    elif j < L_1 + L_2 + L_3:
        c_3 += 1
    else:
        Raise error
于 2013-06-25T23:11:03.827 回答
0

可以使用collections.Counter, and repeatand chainfrom itertools，例如：

from itertools import chain, repeat
from operator import itemgetter
from collections import Counter

junctions = [2,9,15,20]
seq_1 = 'sauron'
seq_2 = 'corrupted'
seq_3 = 'numenor'

indices = list(chain.from_iterable(repeat(i, len(j)) for i, j in enumerate([seq_1, seq_2, seq_3], start=1)))
print Counter(itemgetter(*junctions)(indices))
# Counter({3: 2, 1: 1, 2: 1})
于 2013-06-25T23:31:36.680 回答