Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

nrun for denovo SV signatures #51

Open
AndreVidas opened this issue Aug 26, 2020 · 0 comments
Open

nrun for denovo SV signatures #51

AndreVidas opened this issue Aug 26, 2020 · 0 comments

Comments

@AndreVidas
Copy link

AndreVidas commented Aug 26, 2020

Dear Palimpsest team.

When reading through your Introduction to Palimpsest, I can see that you recommend nrun = 20 for large data set for SBS/DBS/indel mutational signatures analysis. Is this setting also the recommendation for SV mutational signature analysis? I have used Palimpsest v1.0.0 for SV mutational signature extraction so far and here I have usually used nrun = 100 for large data set (1000 samples), since the recommendation at that time was nrun > 20 for large data set if I remember correctly. Nonetheless, If I try with the newest Palimpsest version and setting nrun = 100, I get an error:

There were 50 or more warnings (use warnings() to see the first 50)
Error in [<-.data.frame(*tmp*, , Sig.max.prob.col, value = list(1, :
replacement element 38961 is a matrix/data frame of 1 row, need 78479
Calls: signature_origins -> [<- -> [<-.data.frame
In addition: Warning message:
In [<-.data.frame(*tmp*, , Sig.max.col, value = list("SV_denovo_2.prob", :
provided 78479 variables to replace 1 variables
Execution halted

Because of this error I'm not able to get the output file generated from the signature_origins() function (containing the probability of each event being due to each process). Furthermore the signature_content_plot.pdf looks pretty weird (attached). Apart from that I get all the other plots and output files correctly generated. Here I get 3 SV signatures in total.

image

If I set the nrun = 20 for the same data, I instead get 10 signatures, so I'm wondering if nrun = 20 is a good choice. Also, the signature bars in the plot exceeds the figure (attached) and the signature_content_plot.pdf also looks weird as before.

image

Do you have an idea of what might be causing these problems as well as the best suited nrun?
Best regards,
André

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant