8000 No circles in spider plot · Issue #2 · dariober/cnv_facets · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

No circles in spider plot #2

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.

Sign up for GitHub

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

Jump to bottom
Closed
ahwanpandey opened this issue Feb 5, 2019 · 6 comments
Closed

No circles in spider plot #2

ahwanpandey opened this issue Feb 5, 2019 · 6 comments
8000

Comments

@ahwanpandey
Copy link
ahwanpandey commented Feb 5, 2019 •
edited
Loading

Hello,

Thanks for this tool!

I ran cnv_facets on some whole genome samples with your recommended values. None of the spider plots have circles in them.

cnv_facets.R \
     -t tumor.bam \
     -n normal.bam \
     -vcf common_all_20180423.vcf.gz \
     -o output \
     --gbuild hg19 \
     --cval 25 400 \
     --nbhd-snp 500 \
     --depth 25 4000 \
     --snp-nprocs 8 \
     --snp-baq 10 \
     --snp-mapq 5

I have attached all the plots generated.

Note: I have also run the tool with various combinations of cval [25 400] [25 1000] and [50 1000] and they all have the same problem.

aocs-116-1-3 cval_25_400 nbhd-snp_500 cnv
AOCS-116-1-3.cval_25_400.nbhd-snp_500.cov.pdf
AOCS-116-1-3.cval_25_400.nbhd-snp_500.spider.pdf

Thanks!
@ahwanpandey
Copy link
Author

Also, are there any other parameters you would recommend for WGS data to get cleaner segmentation? Most of my data is 40x Germline and 80x Tumor. I don't really want to loose small focal changes but just looking at the plot, it seems to be pretty noisy.

Thanks!

@ahwanpandey ahwanpandey changed the title No circles in spider plot and "est_insert_size" is NA No circles in spider plot Feb 6, 2019
@ahwanpandey
Copy link
Author

I loaded the pileup created by cnv_facets into the original FACETS R library and was successful in generating the circles in the spider plot. It's just cnv_facets not showing them for some reason.

aocs-116-1-3 logrlogorspider

@dariober
Copy link
Owner
dariober commented Feb 6, 2019

Hi- thanks for reporting this and apologies for the trouble. Could you post the version of cnv_facets from:

cnv_facets.R -v

Also it would be great if you could send me the csv pileup by putting it on e.g. dropbox so I can have a look and see if I can reporduce the bug. (To anonimise the data you can drop the REF and ALT columns with gunzip -c my_pileup.csv.gz | cut -f 1,2,5,6,9,10 -d',' | gzip > my_pileup.x.csv.gz).

The version of facets that is built in cnv_facets may not be the latest available. If that is the case I can rebuild cnv_facets with the latest one.

@ahwanpandey
Copy link
Author

Hey Dario,

Thanks for the reply. Here is the version:

cnv_facets.R -v
0.12.1; facets=0.5.14

And the data is here:

https://www.dropbox.com/sh/xeou3ildyhiqcd6/AAADFy0tKwj_f1EkKAsj_mM8a?dl=0

Thanks!

@ahwanpandey
Copy link
Author

I just manually updated FACETS:

cnv_facets.R -v
0.12.1; facets=0.6.0

Still the same issue, but the plot has a legend now.

AOCS-116-1-3.cval_25_400.nbhd-snp_500.spider.pdf

@dariober
Copy link
Owner
dariober commented Feb 9, 2019

Indeed, this was an embarrassing bug now fixed (apologies). I'm going to upload version 0.13.0 to bioconda and it should be available there in a day or two.

@dariober dariober closed this as completed Feb 9, 2019
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
2 participants
@dariober @ahwanpandey
0