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

failed to obtain the final plot #110

Closed
Gon1976 opened this issue Oct 18, 2022 · 3 comments
Closed

failed to obtain the final plot #110

Gon1976 opened this issue Oct 18, 2022 · 3 comments

Comments

@Gon1976
Copy link

Gon1976 commented Oct 18, 2022

Hi
I install smudgeplot using conda. And also running the exec files directly.
I can obtain kmer-pairs_L_U_coverages.tsv, kmc_tools and smudgeplot.py hetkmers
But, finally I can run the last step to construct the graph (I think a library is missing in my R installation)

Error in library("argparse") : there is no package called ‘argparse’
Calls: suppressPackageStartupMessages -> withCallingHandlers -> library
Ejecución interrumpida

If I install the library in a personal library it will work? Or I need to ask the root to install the package?
thanks

@KamilSJaron
Copy link
Owner

Hi, I am a bit puzzled. I thought conda actually installs the package as well.

In any way, you should be able to install the package yourself and that should fix the problem.

@Gon1976
Copy link
Author

Gon1976 commented Oct 19, 2022

Thanks! I can make the plots!

I a bit confused about how to select low and upper coverage filter. I try with different values according to warning message of smudgeplot:
!! Careful, your coverage filter on the lower end (L = 100 ) is higher than half of the 1n coverage estimate ( 1n / 2 = 54.55. Using this L filter the plot suggest diplod genome, I try using L 20 (pentaplid) and L10 (multiploidy).

How can extract L and U values from genomescope output?
thanks

@KamilSJaron
Copy link
Owner

Hi, check this wiki page in the documentation: https://github.com/KamilSJaron/smudgeplot/wiki/chosing-L-and-U

If it was unclear, i will try to clarify it.

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

2 participants