AlexandrovLab / SigProfilerAssignment

Assignment of known mutational signatures to individual samples and individual somatic mutations
BSD 2-Clause "Simplified" License
46 stars 10 forks source link

SBS96_Decomposition_Plots.pdf is empty #120

Closed odoveestrella closed 6 months ago

odoveestrella commented 6 months ago

Description

Hello! I have run the Analyze.decompose_fit function from SigProfilerAssignment to decompose several SBS96 de novo signatures using a signatures database. There were no problems during the execution, however when exploring the output I realised that the SBS96_Decomposition_Plots.pdf is empty.

Expected behaviour

SBS96_Decomposition_Plots.pdf file containing the decomposition plots for each of the de novo signatures.

Observed behaviour

SBS96_Decomposition_Plots.pdf file is empty.

Steps to reproduce

Please, find below the code I have executed:

from SigProfilerAssignment import Analyzer as Analyze

Analyze.decompose_fit(
    samples = .../snv_filtered,
    output = .../out_sigassign_owndb,
    input_type = "vcf",
    signatures = ../SBS1536toSBS96_De-Novo_Signatures.tsv,
    genome_build = "GRCh38",
    signature_database = ../v3.3_ALL.tsv

And here the contents of JOB_METADATA_SPA.txt

THIS FILE CONTAINS THE METADATA ABOUT SYSTEM AND RUNTIME

-------System Info-------
Operating System Name: Linux
Nodename: 006
Release: 3.10.0-1062.9.1.el7.x86_64
Version: #1 SMP Fri Dec 6 15:49:49 UTC 2019

-------Python and Package Versions------- 
Python Version: 3.9.18
SigProfilerPlotting Version: 1.3.21
SigProfilerMatrixGenerator Version: 1.2.25
SigProfilerAssignment Version: 0.1.4
Pandas version: 1.5.3
Numpy version: 1.26.4

--------------EXECUTION PARAMETERS--------------
INPUT DATA
    input_type: vcf
    output: .../out_sigassign_owndb
    samples: .../snv_filtered
    reference_genome: GRCh38
    context_types: 96
    exome: False
COSMIC MATCH
    cosmic_version: 3.4
    nnls_add_penalty: 0.05
    nnls_remove_penalty: 0.01
    initial_remove_penalty: 0.05
    de_novo_fit_penalty: 0.02
    export_probabilities: True
    collapse_to_SBS96: True
    denovo_refit_option: False
    decompose_fit_option: True
    cosmic_fit_option: False

-------Analysis Progress------- 
[2024-04-04 12:54:02] Analysis started: 

 Decomposing De Novo Signatures  .....
 [2024-04-04 12:54:02] Assigning decomposed signature 

 Finished Decompose fitting! 
Execution time:0:02:26.314548

[2024-04-04 12:56:29] Analysis ended: 

-------Job Status------- 
Assignment of mutational signatures completed successfully! 
Total execution time: 0:02:26 
Results can be found in:  .../out_sigassign_owndb  folder

Thank you!

mdbarnesUCSD commented 6 months ago

Hi @odoveestrella,

Thanks for reaching out with detailed information. Decomposition plots are not currently supported for custom signature database files. We plan to update this behavior in a future release.