-
Notifications
You must be signed in to change notification settings - Fork 208
Open
Labels
bugSomething isn't workingSomething isn't working
Description
Hi there,
I encountered an error when running pyscenic ctx with the following code, and I would like to know the cause of this error. Thank you!
ncore=20
disease_name="MM_with_HD_eRNA_count"
loom_file="/users/cuiyiran/RData/DEG_enhancer/Fig3/loom_file/${disease_name}.loom"
output_base="/home/cuiyiran/data/caixp/SCENIC/output/${disease_name}"
short_name="${disease_name%_eRNA_count}"
step1_out="${output_base}/step1_${short_name}_grn.tsv"
step2_out="${output_base}/step2_${short_name}_ctx.tsv"
auc_out="${output_base}/${short_name}_SCENIC.loom"
##step2 CTX
pyscenic ctx \
${step1_out} \
/home/cuiyiran/software/pySCENIC/hg38_500bp_up_100bp_down_full_tx_v10_clust.genes_vs_motifs.rankings.feather \
/home/cuiyiran/software/pySCENIC/hg38_10kbp_up_10kbp_down_full_tx_v10_clust.genes_vs_motifs.rankings.feather \
--annotations_fname /home/cuiyiran/software/pySCENIC/motifs-v10nr_clust-nr.hgnc-m0.001-o0.0.tbl \
--expression_mtx_fname ${loom_file} \
--mode "dask_multiprocessing" \
--output ${step2_out} \
--num_workers ${ncore}
Traceback (most recent call last):
File "/home/cuiyiran/miniconda3/bin/pyscenic", line 8, in <module>
sys.exit(main())
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/pyscenic/cli/pyscenic.py", line 713, in main
args.func(args)
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/pyscenic/cli/pyscenic.py", line 247, in prune_targets_command
orthologuous_identity_threshold=args.min_orthologous_identity,
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/pyscenic/prune.py", line 434, in prune2df
module_chunksize,
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/pyscenic/prune.py", line 353, in _distributed_calc
num_workers=num_workers if num_workers else cpu_count(),
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/dask/base.py", line 290, in compute
(result,) = compute(self, traverse=False, **kwargs)
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/dask/base.py", line 573, in compute
results = schedule(dsk, keys, **kwargs)
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/dask/multiprocessing.py", line 231, in get
**kwargs,
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/dask/local.py", line 506, in get_async
raise_exception(exc, tb)
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/dask/local.py", line 314, in reraise
raise exc
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/dask/local.py", line 219, in execute_task
result = _execute_task(task, data)
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/dask/core.py", line 119, in _execute_task
return func(*(_execute_task(a, cache) for a in args))
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/dask/core.py", line 119, in <genexpr>
return func(*(_execute_task(a, cache) for a in args))
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/dask/core.py", line 119, in _execute_task
return func(*(_execute_task(a, cache) for a in args))
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/pyscenic/transform.py", line 377, in modules2df
for module in modules
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/pyscenic/transform.py", line 377, in <listcomp>
for module in modules
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/pyscenic/transform.py", line 288, in module2df
db, module, motif_annotations, weighted_recovery=weighted_recovery
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/pyscenic/transform.py", line 176, in module2features_auc1st_impl
df = db.load(module)
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/ctxcore/rnkdb.py", line 135, in load
regions_or_genes_type=self.ct_db.all_region_or_gene_ids.type,
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/ctxcore/ctdb.py", line 789, in subset_to_pandas
self.prefetch(region_or_gene_ids=region_or_gene_ids, engine=engine, sort=True)
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/ctxcore/ctdb.py", line 740, in prefetch
region_or_gene_ids=region_or_gene_ids, sort=sort
File "/home/cuiyiran/miniconda3/lib/python3.7/site-packages/ctxcore/ctdb.py", line 644, in _prefetch_as_pyarrow_table
"CisTargetDatabase object is in an inconsistent state: "
ValueError: CisTargetDatabase object is in an inconsistent state: "region_or_gene_ids_loaded" attribute is None, but "df_cached" is not.
Reactions are currently unavailable
Metadata
Metadata
Assignees
Labels
bugSomething isn't workingSomething isn't working