Hi, i am trying to use beadarray package for my illumina expression data. Thing is whenever I tried readBeadSummaryData function with codes below, just to find empty QC information likie below.
My qcFile is controlprobeprofile.txt acquired from genomestudio. Is there anybody that has similar problem and solution? Thanks in advance!!
Data=readBeadSummaryData(dataFile, qcFile=qcFile, sampleSheet=NULL,
sep="\t", skip=0, ProbeID="Probe_ID",
columns = list(exprs = "AVG_Signal", se.exprs="BEAD_STDERR",
nObservations = "Avg_NBEADS", Detection="Detection Pval"),
qc.sep="\t", qc.skip=0, controlID="TargetID",
qc.columns = list(exprs="AVG_Signal", se.exprs="BEAD_STDERR",
nObservations="Avg_NBEADS", Detection="Detection Pval"),
illuminaAnnotation="Humanv3", dec=".", quote="", annoCols = c("TargetID", "PROBE_ID","SYMBOL"))
ExpressionSetIllumina (storageMode: list)
assayData: 13447 features, 5 samples
element names: exprs, se.exprs, nObservations, Detection
protocolData: none
phenoData
rowNames: 9683758072_E 9683758071_J ...
9999593135_G (5 total)
varLabels: sampleID
varMetadata: labelDescription
featureData
featureNames: 60689 1580504 ... NEGATIVE (13447
total)
fvarLabels: ProbeID TargetID ... Status (5 total)
fvarMetadata: labelDescription
experimentData: use 'experimentData(object)'
Annotation: Humanv3
QC Information
Available Slots:
QC Items:
sampleNames: