EXOMECNV segmentation algorithm
0
0
Entering edit mode
8.0 years ago
ksivapra • 0

I am calling the CNVs from canine data using ExomeCNV. The organism has 38 pairs of chromosome and so, the CNV plot generated at the end could not acomodate all. Therefore, I am running the Rscript in batches of 2 where the each set has 20 chromosome. while running the second batch, I came across the following error: Error in 1:max(ecnv$copy.number, na.rm = TRUE) : result would be too long a vector Calls: do.plot.eCNV -> do.plot.one.eCNV In addition: Warning message: In max(ecnv$copy.number, na.rm = TRUE) : no non-missing arguments to max; returning -Inf Execution halted

This was encountered for just one sample out of 8. I narrowed down the problematic chromosome and found that it has 4347 calls which is NOT different from other samples. But I could not find the exact reason. Can anyone help me interpret and resolve it?

R • 1.5k views
ADD COMMENT

Login before adding your answer.

Traffic: 2810 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6