Hello,

I came across MAESTRO(Model-based AnalysEs of Single-cell Transcriptome and RegulOme) and wanted to use it on my data.

I have used the instructions on github and created a conda environment for MAESTRO.

Github link for MAESTRO: github.com/liulab-dfci/MAESTRO

However, I am encountering an error. Such as:

[Tue Jun 15 16:24:03 2021]

Error in rule scrna_analysis:

jobid: 2

output: Result/Analysis/MAESTRO_DiffGenes.tsv, Result/Analysis/MAESTRO_cluster.png, Result/Analysis/MAESTRO_annotated.png, Result/Analysis/MAESTRO.PredictedTFTop10.

shell:

python /home/k/anaconda3/envs/MAESTRO/lib/python3.8/site-packages/MAESTRO/lisa_path.py --species GRCh38 --input /home/k/sam/RAID/chr/L_M/A_P/hg38_2.1.tar.gz; Rscript /home/k/anaconda3/envs/MAESTRO/lib/python3.8/site-packages/MAESTRO/R/scRNAseq_pipe.R --expression ../QC/MAESTRO_filtered_gene_count.h5 --species GRCh38 --prefix MAESTRO --signature human.immune.CIBERSORT --outdir Result/Analysis --thread 10

(one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)

I tried updating MAESTRO, but I still got the error. I was wondering if anyone encountered similar errors and if so, what steps did you follow so that the pipeline takes place without errors?

Extra question: For annotating the cells, MAESTRO uses built-in signatures like [human.immune.CIBERSORT]. But for my data, the immune cell types won't really fit well. However, I would still like to know the list of cell types that MAESTRO is using. How can I obtain that?

Thank you!
(This is my first time posting so if the format is not okay, apologies)



Source link