<p>用<a href="https://github.com/biocore-ntnu/pyranges" rel="nofollow noreferrer">pyranges</a>回答,这基本上是熊猫撒上了生物信息学的糖。在</p>
<p>设置:</p>
<pre><code>import numpy as np
np.random.seed(0)
import pyranges as pr
a = pr.random(int(1e6))
# + + -+ -+ +
# | Chromosome | Start | End | Strand |
# | (category) | (int32) | (int32) | (category) |
# | + -+ -+ |
# | chr1 | 8830650 | 8830750 | + |
# | chr1 | 9564361 | 9564461 | + |
# | chr1 | 44977425 | 44977525 | + |
# | chr1 | 239741543 | 239741643 | + |
# | ... | ... | ... | ... |
# | chrY | 29437476 | 29437576 | - |
# | chrY | 49995298 | 49995398 | - |
# | chrY | 50840129 | 50840229 | - |
# | chrY | 38069647 | 38069747 | - |
# + + -+ -+ +
# Stranded PyRanges object has 1,000,000 rows and 4 columns from 25 chromosomes.
# For printing, the PyRanges was sorted on Chromosome and Strand.
b = pr.random(int(1e6), length=1)
# + + -+ -+ +
# | Chromosome | Start | End | Strand |
# | (category) | (int32) | (int32) | (category) |
# | + -+ -+ |
# | chr1 | 52110394 | 52110395 | + |
# | chr1 | 122640219 | 122640220 | + |
# | chr1 | 162690565 | 162690566 | + |
# | chr1 | 117198743 | 117198744 | + |
# | ... | ... | ... | ... |
# | chrY | 45169886 | 45169887 | - |
# | chrY | 38863683 | 38863684 | - |
# | chrY | 28592193 | 28592194 | - |
# | chrY | 29441949 | 29441950 | - |
# + + -+ -+ +
# Stranded PyRanges object has 1,000,000 rows and 4 columns from 25 chromosomes.
# For printing, the PyRanges was sorted on Chromosome and Strand.
</code></pre>
<p>执行:</p>
^{pr2}$