SummarizedExperiment 0.2.11

Container to represent data from genomic experiments

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/.