beadarray: summarizing iScan data is not working
0
0
Entering edit mode
S. K. • 0
@s-k-12879
Last seen 7.4 years ago
Germany, Heidelberg

Hello,

I am trying to read Illumina HT12-v4 beadarrays with beadarray that were scanned with an iScan system. All necessary files are residing in the current working directory and I am not using a sample file.

library(beadarray)
processSwathData()
ExpBData <- readIllumina(illuminaAnnotation = "Humanv4")
ExpSetData <- summarize(ExpBData,useSampleFac = T)

Summarizing seems to work until the annotation step... then I get the following error message:

Making  summary object
Annotating control probes using package illuminaHumanv4.db Version:1.26.0

Error in value[[3L]](cond) : duplicate row.names: 200647550006
  AnnotatedDataFrame 'initialize' could not update varMetadata:
  perhaps pData and varMetadata are inconsistent?

Removing the useSampleFacs argument to summarize makes the error go away, but then I get an object that has separate entries for each Swath - so 24 "chips" instead of 12.

Can somebody point me in the right direction on how to get correct results?

Thanks,

Simon

R version 3.3.3 (2017-03-06)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Debian GNU/Linux 8 (jessie)

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

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

other attached packages:
 [1] illuminaHumanv4.db_1.26.0 org.Hs.eg.db_3.4.0        AnnotationDbi_1.36.2     
 [4] IRanges_2.8.2             S4Vectors_0.12.2          beadarray_2.24.0         
 [7] ggplot2_2.2.1             Biobase_2.34.0            BiocGenerics_0.20.0      
[10] illuminaio_0.16.0        

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.10         XVector_0.14.1       magrittr_1.5         GenomicRanges_1.26.4
 [5] zlibbioc_1.20.0      munsell_0.4.3        colorspace_1.3-2     stringr_1.2.0       
 [9] GenomeInfoDb_1.10.3  plyr_1.8.4           tools_3.3.3          grid_3.3.3          
[13] base64_2.0           gtable_0.2.0         DBI_0.6-1            lazyeval_0.2.0      
[17] openssl_0.9.6        digest_0.6.12        tibble_1.3.0         reshape2_1.4.2      
[21] bitops_1.0-6         RCurl_1.95-4.8       memoise_1.0.0        RSQLite_1.1-2       
[25] limma_3.30.13        stringi_1.1.5        BeadDataPackR_1.26.0 scales_0.4.1
beadarray iscan • 1.2k views
ADD COMMENT
0
Entering edit mode

Just a wild guess but might be related to another problem of mine due to longer chip identifiers since last October:

beadarray: readIllumina with sample sheet broken due to Illumina switching to 12 digit chip identifiers

ADD REPLY

Login before adding your answer.

Traffic: 1262 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