SummarizedExperiment

Container to represent genomic experiments, follows Bioconductor's SummarizedExperiment.

Install

Package is published to PyPI

pip install summarizedexperiment

Usage

Currently supports SummarizedExperiment & RangeSummarizedExperiment classes

First create necessary sample data

import pandas as pd
import numpy as np
from genomicranges import GenomicRanges

nrows = 200
ncols = 6
counts = np.random.rand(nrows, ncols)
df_gr = pd.DataFrame(
    {
        "seqnames": [
            "chr1",
            "chr2",
            "chr2",
            "chr2",
            "chr1",
            "chr1",
            "chr3",
            "chr3",
            "chr3",
            "chr3",
        ]
        * 20,
        "starts": range(100, 300),
        "ends": range(110, 310),
        "strand": ["-", "+", "+", "*", "*", "+", "+", "+", "-", "-"] * 20,
        "score": range(0, 200),
        "GC": [random() for _ in range(10)] * 20,
    }
)

gr = GenomicRanges.fromPandas(df_gr)

colData = pd.DataFrame(
    {
        "treatment": ["ChIP", "Input"] * 3,
    }
)

To create a SummarizedExperiment,

from summarizedexperiment import SummarizedExperiment

tse = SummarizedExperiment(
    assays={"counts": counts}, rowData=df_gr, colData=colData
)

To create a RangeSummarizedExperiment

from summarizedexperiment import RangeSummarizedExperiment

trse = RangeSummarizedExperiment(
    assays={"counts": counts}, rowRanges=gr, colData=colData
)

For more use cases including subset, checkout the documentation

Note

This project has been set up using PyScaffold 4.1.1. For details and usage information on PyScaffold see https://pyscaffold.org/.