aertslab / scenicplus

SCENIC+ is a python package to build gene regulatory networks (GRNs) using combined or separate single-cell gene expression (scRNA-seq) and single-cell chromatin accessibility (scATAC-seq) data.
Other
178 stars 28 forks source link

No menr.pkl produced by `run_pycistarget` #113

Closed alexlenail closed 1 year ago

alexlenail commented 1 year ago

The scenic+ tutorial says to run

from scenicplus.wrappers.run_pycistarget import run_pycistarget
run_pycistarget(
    region_sets = region_sets,
    species = 'homo_sapiens',
    save_path = os.path.join(work_dir, 'motifs'),
    ctx_db_path = rankings_db,
    dem_db_path = scores_db,
    path_to_motif_annotations = motif_annotation,
    run_without_promoters = True,
    n_cpu = 8,
    _temp_dir = os.path.join(tmp_dir, 'ray_spill'),
    annotation_version = 'v10nr_clust',
    )

and then

import dill
menr = dill.load(open(os.path.join(work_dir, 'motifs/menr.pkl'), 'rb'))

However, there is no menr.pkl file when I run run_pycistarget as specified. Judging from the code,

https://github.com/aertslab/scenicplus/blob/61c17c842d25e537915f4caf7c604179e979603f/src/scenicplus/wrappers/run_pycistarget.py#L214-L216

It looks like this file might get created if the save_partial flag is set to true, but it's false by default:

https://github.com/aertslab/scenicplus/blob/61c17c842d25e537915f4caf7c604179e979603f/src/scenicplus/wrappers/run_pycistarget.py#L74-L75

alexlenail commented 1 year ago

Oops nevermind