getting error in processSwathData command from the beadarray package
2
0
Entering edit mode
@harpreetsaini-7493
Last seen 9.7 years ago
United Kingdom

Hi,

I am using the "processSwathData" command from the beadarray package to
process Illumina iScan files and I am getting the following error:

processSwathData("./3999538058", outputDir="./3999538058", twoColour=NULL, textstring="_perBeadFile.txt", segmentHeight=326, segmentWidth=397, fullOutput=T, newTextString=".txt")

Error in unlist(scan(file, sep = sep, what = list(integer(), integer(),  :
  error in evaluating the argument 'x' in selecting a method for function
'unlist': Error in scan(file, what, nmax, sep, dec, quote, skip, nlines,
na.strings,  :
  scan() expected 'a real', got '8880.3552.084'

Any idea? Any help would be greatly appreciated.

 

Here is the sessionInfo()

R version 3.1.1 (2014-07-10)
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=en_GB.UTF-8  
 [7] LC_PAPER=en_GB.UTF-8       LC_NAME=C                
 [9] LC_ADDRESS=C               LC_TELEPHONE=C           
[11] LC_MEASUREMENT=en_GB.UTF-8 LC_IDENTIFICATION=C      

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

other attached packages:
[1] beadarray_2.14.1    ggplot2_1.0.1       Biobase_2.24.0    
[4] BiocGenerics_0.10.0

loaded via a namespace (and not attached):
 [1] AnnotationDbi_1.26.1 base64_1.1           BeadDataPackR_1.16.0
 [4] colorspace_1.2-6     DBI_0.3.1            digest_0.6.8       
 [7] GenomeInfoDb_1.0.2   GenomicRanges_1.16.4 grid_3.1.1         
[10] gtable_0.1.2         illuminaio_0.6.1     IRanges_1.22.10    
[13] limma_3.20.9         MASS_7.3-39          munsell_0.4.2      
[16] plyr_1.8.1           proto_0.3-10         Rcpp_0.11.5        
[19] reshape2_1.4.1       RSQLite_1.0.0        scales_0.2.4       
[22] stats4_3.1.1         stringr_0.6.2        XVector_0.4.0      


Thanks,
Harpreet

microarray beadarray illumina • 1.8k views
ADD COMMENT
0
Entering edit mode

Hi Harpreet,

Sorry you're having issues with this function.  I suspect this may be to do with the character used for decimal points.  Is there anywhere you can host the "_perBeadFile.txt" file so that I can have look at it?

Thanks,

Mike

ADD REPLY
0
Entering edit mode
@harpreetsaini-7493
Last seen 9.7 years ago
United Kingdom

Hi Mike,

Many thanks for your prompt reply.

the "_perBeadFile.Txt" is here,

https://www.dropbox.com/sh/he86j5wjdwrsrny/AAD_BVBNnWXo23foCHFzlW9Ba?dl=0

Let me know if there is any problem in accessing the file.

Thanks,

Harpreet

ADD COMMENT
0
Entering edit mode

Hi Harpreet,

Thanks for uploading that.  It looks like the file you've put there is the "_beadTypeFile.txt", which contains summarised values for each bead type and isn't intended to work with the processSwath() function.

Do you have one that ends with "_perBeadFile.Txt"?  These contain the bead-level data that beadarray works with.

Thanks,

Mike

ADD REPLY
0
Entering edit mode

Apologies about this.

I have now uploaded the "_perBeadFile.txt" file.

https://www.dropbox.com/sh/he86j5wjdwrsrny/AAD_BVBNnWXo23foCHFzlW9Ba?dl=0

Many thanks,

Harpreet

ADD REPLY
0
Entering edit mode

Hi Mike,

I forgot to mention that I have uploaded two "_perBeadFile.txt" files, one for the "G" array and one for the "H" array.

The "processSwathData" command worked fine with the "G" array , but gave the error message with the "H" array.

Many thanks for looking into this.

Harpreet

 

ADD REPLY
0
Entering edit mode
Mike Smith ★ 6.6k
@mike-smith
Last seen 4 hours ago
EMBL Heidelberg

Hi Harpreet,

The problem seems to be with line 969239 in "3999538058_H_perBeadFile.txt".  For some reason the fourth column contains the entry 8880.3552.084, which is clearly not a coordinate.  The simplest fix is to manually edit that entry to be 8880.355 (there's almost always 7 sig. fig.  in those columns).  The file can then be read normally.

However, you might want to do some more digging.  There's a blank entry at line 969342 and the very final line of the file ends abruptly.  It looks to me like several thousand lines have been cut out and pasted into the end of the file.  If you have another copy of the file somewhere you might want to check that there hasn't been some corruption introduced in copying it between devices.  

Mike

ADD COMMENT
0
Entering edit mode

Thanks Mike.

I have manually edit the file and rerunning the command. I may ask the service provider to re-scan the files, hope the new files will be okay.

 

Thanks again,

Harpreet

 

ADD REPLY

Login before adding your answer.

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