[BioC] Beadarray import of iScan generated Illumina gene expression data
Aliaksei Holik
salvador at bio.bsu.by
Tue Dec 6 11:56:23 CET 2011
Dear members,
This is my first attempt to analyse a microarray and I would greatly
appreciate some help.
The data was generated using MOUSEREF-8_V2 BeadChip by a third party. I
provided RNA and got a CD back containing the data. Judging by the
presence of Swath 1 and 2 files, arrays were scanned using iScan system,
which prompts processing with processSwathData. As far as I understand,
I need three types of files for this processing, including *swath*.tiff
files, locs files and beadlevel.txt file. That's where the problem begins.
1) All images I received are in jpeg format. And even though Image
header files contain the line <Compressed>false</Compressed> I'm not
sure I'd be able to import the data into beadarray using these images.
2) I seem to miss .locs files. The files I've got in each chip's
directory (I had 2 chips hybridised in total) are:
- various *.jpeg files including files containing swath in the name
- ChipName_Sample_Grn.idat files
- Effective.cfg file
- Several Metrics.txt file (it appears there were several rounds of
scanning)
- Image header ChipName_Sample.xml files
- ChipName.sdf file
- ChipName_qc.txt file
- and finally 8 ChipName_Sample_perBeadFile.txt files of the following
format:
Code Grn GrnX GrnY
10008 102 621.3782 13898.54
10008 86 1875.903 23202.34
Etc...
My questions are:
- Do I need to wait for tiff files or can I use available jpegs?
- As far as I understand, locs files contain bead coordinates, but these
seem to be present in beadlevel txt files. I wonder, therefore, if
there's a workaround for the lack of locs files. Alternatively, can I
import data straight into a LumiBatch object providing the format of my
beadlevel txt file?
- Providing I go along with using beadarray for data import, how do I
find swathOverlap parameter required by processSwathData. I found
segment width and height in sdf file (only because they were
conveniently the same size as in the command manual), but can't find
anything remotely similar to the overlap parameter.
Your help is greatly appreciated.
Aliaksei.
More information about the Bioconductor
mailing list