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