SummarizedExperiment 0.3.1

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 & RangedSummarizedExperiment 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.from_pandas(df_gr)

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

To create a SummarizedExperiment,

from summarizedexperiment import SummarizedExperiment

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

To create a RangedSummarizedExperiment

from summarizedexperiment import RangedSummarizedExperiment

trse = RangedSummarizedExperiment(
    assays={"counts": counts}, row_ranges=gr, col_data=colData
)

For more examples, checkout the documentation.

Note

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