def bg():
bg = pr.load_dataset("chipseq_background")
return bg
import pyranges as pr
gr = pr.load_dataset("epigenome_roadmap")
rle = gr["chr1"].coverage()
print(list(rle.runs)[:20])
print(list(rle.values)[:20])
raise
import numpy as np
from pyrle import Rle
import pandas as pd
r = pd.Series([1, 2, 3, 4], dtype=np.int16)
# v = pd.Series([-1, 2.3, 3, 4.976], dtype=np.float)
r1 = Rle(r, r)
r2 = Rle(r * 2, r * 2)
# > r2
# numeric-Rle of length 20 with 4 runs
# Lengths: 2 4 6 8
# Values : 2 4 6 8
# > r4
# numeric-Rle of length 20 with 5 runs
# Lengths: 1 2 3 4 10
# Values : 1 2 3 4 0
# > r2 + r4
# numeric-Rle of length 20 with 7 runs
# Lengths: 1 1 1 3 4 2 8
# Values : 3 4 6 7 10 6 8
r3 = r1 + r2
def cs():
cs = pr.load_dataset("chipseq")
return cs
def chipseq_dataset():
gr = pr.load_dataset("chipseq")
return gr
