lumiR {lumi} | R Documentation |
Read in Illumina expression data. We assume the data was saved in a comma or tab separated text file.
lumiR(fileName, sep = NULL, detectionTh = 0.01, na.rm = TRUE, lib = NULL)
fileName |
fileName of the data file |
sep |
the separation character used in the text file. |
detectionTh |
the p-value threshold of determining detectability of the expression. The detection value less than the detectionTh is claimed as detectable. |
na.rm |
determine whether to remove NA |
lib |
the annotation library used by addNuId2lumi |
The function can automatically determine the separation character if it is Tab or comma. Otherwise, the user should specify the separator manually. If the annotation library is provided, the Illumina Id will be replaced with nuID, which is used as the index Id for the lumi annotation packages.
Current version of lumiR can adaptively read the output of BeadStudio Verson 2 and 3. The format Version 3 made quite a few changes comparing with previous versions. One change is the detection value. It was called detectable when the detection value is close to one for Version 2 format. However, the detection value became a p-value in the Version 3. As a result, the detectionTh is automatically changed based on the version. The detectionTh 0.01 for the Version 3 will be changed as the detectionTh 0.99 for Version 2. Another big change is that Version 3 separately output the control probe (gene) information and a "Samples Table". As a result, the controlData slot in LumiBatch class was added to keep the control probe (gene) information, and a QC slot to keep the quality control information, including the "Sample Table" output by BeadStudio version 3.
return a LumiBatch object
Simon Lin, Pan Du
## specify the file name # fileName <- 'Barnes_gene_profile.txt' # Not Run ## load the data # x.lumi <- lumiR(fileName)