[BioC] error in calculate integer copy numbers. bug?

Johnbosco [guest] guest at bioconductor.org
Fri Feb 7 17:11:46 CET 2014


I used getSegmentReadCountsFromBAM() to get a GRanges object from a list of tumor-normal bam files
Then I used referencecn.mops(). Whenever I try to calculate integer copy numbers using calcIntegerCopyNumbers() I get an error. This error can be reproduced using the data provided by cn.mops so this might be a bug? 

>data(cn.mops)
>resRef2 <- referencecn.mops (XRanges[,1], XRanges[,2])
Normalizing...
Starting local modeling, please be patient...
Reference sequence:  chrA
Starting segmentation algorithm...
Using "fastseg" for segmentation.
> resRef5 <- calcIntegerCopyNumbers(resRef2)
Error in `colnames<-`(`*tmp*`, value = "S_1") :
  length of 'dimnames' [2] not equal to array extent 

 -- output of sessionInfo(): 

R version 3.0.2 (2013-09-25)
Platform: x86_64-redhat-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets  methods  
[8] base     

other attached packages:
[1] snow_0.3-13          cn.mops_1.8.3        GenomicRanges_1.14.4
[4] XVector_0.2.0        IRanges_1.20.6       Biobase_2.22.0      
[7] BiocGenerics_0.8.0   BiocInstaller_1.12.0

loaded via a namespace (and not attached):
[1] Biostrings_2.30.1 bitops_1.0-6      Rsamtools_1.14.2  stats4_3.0.2     
[5] tools_3.0.2       zlibbioc_1.8.0 

--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list