[BioC] ExomeCopy package error

nat [guest] guest at bioconductor.org
Tue Dec 18 17:04:56 CET 2012


Hi,
I am testing ExomeCopy on a set of illumina sequencing data generated after agilent custum gene set  capture, only  35 gene exons overall (~ 300 Mb total).
1-My first question would be : is the method reliable for smaller capture?
2-Secondly.,
I sucessfully used  the method on the data provided with the vignette. I then tried with my own set of data
in total 10 samples, across  2282 rows , I have  calculated the GC and bg values nad went on to the next step where I ran exomeCopy on one sample and one chromosome

###this seems to work for certains samples, expect the warning for small subset
> fit <- exomeCopy(rdata["Chr10"], sample.name = exome.samples[10],
+ X.names = c("bg", "GC", "GC.sq", "width"), S = 0:6,
+ d = 2)
Warning message:
In exomeCopy(rdata["Chr10"], sample.name = exome.samples[10], X.names = c("bg",  :
  exomeCopy was tested for thousands of ranges covering a targeted region on a single chromosome.  The results might not be reliable for less than a hundred ranges.


###then I try with another sample and got this error message 
Error in if (goto.cnv < 1/(nstates - 1)) { : 
  missing value where TRUE/FALSE needed

I am not sure where this comes from and would like some help on this

>  fit <- exomeCopy(rdata["Chr10"], sample.name = exome.samples[7],
+ X.names = c("bg", "GC", "GC.sq", "width"), S = 0:6,
+ d = 2)
Error in if (goto.cnv < 1/(nstates - 1)) { : 
  missing value where TRUE/FALSE needed
In addition: Warning message:
In exomeCopy(rdata["Chr10"], sample.name = exome.samples[7], X.names = c("bg",  :
  exomeCopy was tested for thousands of ranges covering a targeted region on a single chromosome.  The results might not be reliable for less than a hundred ranges.
> exome.samples
 [1] "sample9"  "sample8"  "sample7"  "sample6"  "sample5"  "sample1" 
 [7] "sample16" "sample15" "sample14" "sample13"

 
example if my data

#####
RangedData with 2282 rows and 15 value columns across 18 spaces
        space                 ranges   |   sample9   sample8   sample7
     <factor>              <IRanges>   | <integer> <integer> <integer>
1        Chr1 [115247069, 115247168]   |      4369      4446      3675
2        Chr1 [115247169, 115247268]   |      8222      8081      6978
3        Chr1 [115247269, 115247368]   |      9548     10026      8690

 -- output of sessionInfo(): 

R version 2.15.0 (2012-03-30)
Platform: x86_64-unknown-linux-gnu (64-bit)

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

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

other attached packages:
[1] exomeCopy_1.2.0      Rsamtools_1.8.6      Biostrings_2.24.1   
[4] GenomicRanges_1.8.13 IRanges_1.14.4       BiocGenerics_0.2.0  

loaded via a namespace (and not attached):
[1] bitops_1.0-4.1 stats4_2.15.0  tools_2.15.0   zlibbioc_1.2.0


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



More information about the Bioconductor mailing list