我想生成所有生成的排列并將 itertools.product() 保存到一個數(shù)組中。目前，我正在嘗試獲取字符“ACGT”的所有長度為 2 的排列。當(dāng)我嘗試使用 numpy.asarray() 時，僅保存最終排列 (['T', 'T'])，我認(rèn)為這是因為每次在數(shù)組中都會覆蓋該條目。我試過以下，import itertools as itimport numpy as npfor x in it.product('ACGT', repeat=2): array = np.asarray(x)print(array)['T', 'T']后來我想為更大的“單詞”長度做這個，但是當(dāng)我只期望 16 個結(jié)果時更容易測試。如果我使用的是 R，我會創(chuàng)建一個空向量并按順序添加到向量中……但是，我仍在嘗試掌握 Python 的竅門……請指教！
1 回答
胡說叔叔
TA貢獻(xiàn)1804條經(jīng)驗 獲得超8個贊
這不起作用的原因是因為您為該結(jié)果創(chuàng)建了一個it.product(..)數(shù)組。不是為了整個結(jié)果。
您可以使用以下方法創(chuàng)建這樣的矩陣：
np.array(list(it.product('ACGT', repeat=2)))
或與meshgrid：
dna = np.array(list('ACGT'))
np.transpose(np.meshgrid(dna, dna)).reshape(-1,2)
兩者都產(chǎn)生一個如下所示的數(shù)組：
array([['A', 'A'],
['A', 'C'],
['A', 'G'],
['A', 'T'],
['C', 'A'],
['C', 'C'],
['C', 'G'],
['C', 'T'],
['G', 'A'],
['G', 'C'],
['G', 'G'],
['G', 'T'],
['T', 'A'],
['T', 'C'],
['T', 'G'],
['T', 'T']], dtype='<U1')
添加回答
舉報
0/150
提交
取消