Question

Detect CNVs in target sample

0

Entering edit mode

ashkot09 • 0

@ashkot09-9452

Last seen 8.3 years ago

Hi all,

Using cn.MOPS for CNV detection and it requires optimally 6 samples (6 BAM files) to effectively run the analysis. I am generating CNV plots for each chromosome but I believe it outlines the CNVs for all samples.

What I want is generate plots for CNVs detected in any one sample i.e. generate CNV plots for every chromosome for any sample. How can that be done?

Thanks in advance.

CNV cn.mops • 1.6k views

ADD COMMENT • link updated 8.3 years ago by Günter Klambauer ▴ 540 • written 8.3 years ago by ashkot09 • 0

score 0 · Answer 1 · 2016-01-13

0

Entering edit mode

Günter Klambauer ▴ 540

@gunter-klambauer-5426

Last seen 3.2 years ago

Austria

Have you tried segplot(resCNMOPS,...) ? This function should do exactly what you request...

ADD COMMENT • link 8.3 years ago Günter Klambauer ▴ 540

0

Entering edit mode

Thanks, this worked, how to save the segplot as an image?

ADD REPLY • link 8.3 years ago ashkot09 • 0

1

Entering edit mode

pdf("mySegplot.pdf")

segplot(resCNMOPS) # or: segplot(resCNMOPS,toFile=TRUE)

dev.off()

ADD REPLY • link 8.3 years ago Günter Klambauer ▴ 540

0

Entering edit mode

Hi,

if I want to plot all the chromosomes together in one graph, what argument should I set? I can not find it in help. it only says

plot.type

the type of plot. (Default = "s"). what other options? Thanks!

ADD REPLY • link 7.4 years ago tangming2005 ▴ 190