readQC {beadarray} | R Documentation |
Reads the standard format of Illumina quality control information and produces diagnostic plots
readQC(file, columns = list(exprs = "AVG_Signal", NoBeads = "Avg_NBEADS", Detection="Detection", BeadStDev="BEAD_STDERR", Narrays="NARRAYS", arrayStDev = "ARRAY_STDEV", controlID = "ProbeID", controlType="TargetID"),sep="\t",skip=7,header=T)
file |
name of file containing qc information |
columns |
a vector of column names to read from the file |
skip |
number of lines of header information to ignore in the file |
header |
if TRUE the column names in the file are read |
sep |
a character string for the file separator |
The format of the quality control files differs slightly between BeadStudio versions 1 and 2. Version 1 of the software gives one averaged value for each control type, whereas version 2 gives suumarised values for each control of a particular type. The user does not have to know the version of BeadStudio used to generate the file.
Once read in, the quality control information can be very useful for diagnostic purposes.
readQC produces an assayData object with list items defined by the columns parameter.
Mark Dunning
##Code to read the example quality control file included with the #package. #QC = readQC(sep=",", skip=7, columns=list(exprs="AVG.Signal", BeadStDev="SeqVAR")) #the average expression of each control can then be accessed by the $ operator #QC$exprs