单样本GSEA分析

单样本基因集富集分析（ssGSEA）是GSEA算法的单样本扩展。它为每个样本和基因集配对计算一个独立的富集分数[1]。

安装最新版本

pip install single_sample_gsea

用法

>>> from single_sample_gsea import ss_gsea

>>> gene_sets = {
    "gs1": {"gene2", "gene3"},
    "gs2": {"gene1", "gene4"},
    }

>>> data = {
    "gene": ["gene1", "gene2", "gene3", "gene4", "gene5"],
    "sample-1": [1, 3, 4, 7, 32],
    "sample-2": [25, 4, 6, 18, 1],
    }
>>> data = pd.DataFrame(data).set_index("gene")

>>> ss_gsea(data, gene_sets)
               gs1       gs2
sample-1 -1.333333 -0.962974
sample-2 -1.333333  2.543214

参考文献

[1] Barbie, D., Tamayo, P., Boehm, J. et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature 462, 108–112 (2009). https://doi.org/10.1038/nature08460