Beispiel #1
0
def enrichment_no_t_task(
    exp,
    block,
    T,
    gs,
    patterns,
    base_filename,
):

    if settings.CELERY_DEBUG:
        import sys
        sys.path.append(
            '/Migration/skola/phd/projects/miXGENE/mixgene_project/wrappers/pycharm-debug.egg'
        )
        import pydevd
        pydevd.settrace('localhost',
                        port=6901,
                        stdoutToServer=True,
                        stderrToServer=True)
    gene_set = gs.get_gs()
    patterns = patterns.get_gs()
    e = EnrichmentInGeneSets(patterns.genes)
    enrich = e.getModuleEnrichmentInGeneSets(patterns.genes,
                                             gene_set.genes,
                                             pval_threshold=T)
    enrich = dict(
        (mod, (genes,
               map(lambda x: (gene_set.description[x[0]], x[0], x[1]), terms)))
        for (mod, (genes, terms)) in enrich.items())
    ds = DictionarySet(exp.get_data_folder(), base_filename)
    ds.store_dict(enrich)
    return [ds], {}
def enrichment_no_t_task(exp, block,
                     T,
                     gs,
                     cs,
                     base_filename,
    ):

    if settings.CELERY_DEBUG:
        import sys
        sys.path.append('/Migration/skola/phd/projects/miXGENE/mixgene_project/wrappers/pycharm-debug.egg')
        import pydevd
        pydevd.settrace('localhost', port=6901, stdoutToServer=True, stderrToServer=True)
    gene_set = gs.get_gs()
    cs = cs.load_set()
    e = EnrichmentInGeneSets(cs)
    enrich = e.getModuleEnrichmentInGeneSets(cs, gene_set.genes, pval_threshold=T)
    enrich = dict((mod, (genes, map(lambda x: (gene_set.description[x[0]], x[0], x[1]), terms))) for (mod, (genes, terms)) in enrich.items())
    ds = DictionarySet(exp.get_data_folder(), base_filename)
    ds.store_dict(enrich)
    return [ds], {}