Question: Is there a way to plot both tumor and normal samples' copy number segments on the same plot using CNVkit scatter
gravatar for asjannu
2.9 years ago by
asjannu0 wrote:

Hi, I have been using CNVkit for 3 months now and discovering few interesting things.

These are my questions: 1) How can plot copy number segments of paired tumor and normal samples on the same scatter plot? Something like this enter image description here

2) I have done Segmentation and Segment Calling. Applied the default threshold values (-t=-1.1,-0.4,0.3,0.7) at the 'Call' step. But the number of Copy number positions did not change with the applied threshold. For example, My .segments.cns file had 100 CNV positions with the applied threshold, there should be a decrease in the number. But that didn't happen!! Am I doing something wrong??

3) How can use -t (threshold) at segmentation step? Does it need 4 threshold values like the call step?

I'd really appreciate if someone could answer my questions.

Thank you!

cnvkit • 1.1k views
ADD COMMENTlink modified 2.8 years ago by Eric T.2.6k • written 2.9 years ago by asjannu0
gravatar for Eric T.
2.8 years ago by
Eric T.2.6k
San Francisco, CA
Eric T.2.6k wrote:
  1. This feature isn't available on the command line, currently. It would be a nice feature, and you're not the first to request it. I've opened a ticket for it here:

However, within the Python interpreter or in a script you can use cnvlib.do_scatter to generate a plot with one segment set, and then use matplotlib functions and/or functions within cnvlib.scatter and cnvlib.plots to add another set of horizontal spans to the same matplotlib axis object.

  1. To merge segments with the same inferred integer copy number, you can use call --filter cn. Otherwise, the segment breakpoints are left alone and each segment is assigned an integer copy number (cn value, in the .cns file) independently. See here:

  2. It counts up by ploidy: the first threshold value specifies the cutoff between cn=0 and cn=1, the second specifies cn=1 vs. cn=2, etc. And beyond the specified cutoffs, for higher copy numbers, it defaults to the cutoffs that would be implied for log2 values in a fully clonal/germline/single-cell diploid sample.

ADD COMMENTlink written 2.8 years ago by Eric T.2.6k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1752 users visited in the last hour