1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
|
import pyranges as pr
import numpy as np
np.random.seed(0)
def test_stranded():
cpg = pr.data.cpg()
exons = pr.data.exons()
j = cpg.join(exons)
assert j.stranded
j.Strand = "."
assert not j.stranded
j.Strand = np.random.choice("+ -".split(), size=len(j))
assert j.stranded
for _, df in j:
assert len(df.Strand.drop_duplicates()) == 1
def test_unstrand():
exons = pr.data.exons()
cpg = pr.data.cpg()
print(exons)
print(exons.columns)
x = exons.unstrand()
print(x.stranded)
print(x)
print(x.columns)
for _, df in x:
assert not df.index.duplicated().sum()
# x = pr.PyRanges(x.df.reset_index(drop=True))
cpg.join(x)
|
