[BioC] beadarray example data problem
Keith Satterley
keith at wehi.EDU.AU
Fri Oct 20 04:06:52 CEST 2006
Hi Mark or anyone else who can help with the beadarray package,
I'm having trouble getting the beadarray example data to work as described:
> sessionInfo()
R version 2.4.0 (2006-10-03)
i386-pc-mingw32
locale:
LC_COLLATE=English_Australia.1252;LC_CTYPE=English_Australia.1252;LC_MONETARY=English_Australia.1252;LC_NUMERIC=C;LC_TIME=English_Australia.1252
attached base packages:
[1] "tools" "methods" "stats" "graphics" "grDevices" "utils"
"datasets" "base"
other attached packages:
beadarray affy affyio Biobase limma
"1.2.1" "1.12.0" "1.2.0" "1.12.0" "2.9.1"
When I load the BSData item and look at BSData, I have the following:
=====================
> data(BSData)
> BSData
ExpressionSetIllumina (storageMode: list)
assayData: 47293 features, 18 samples
element names: exprs, BeadStDev, NoBeads, Detection
phenoData
rowNames: I.1, IC.1, ..., P42.2 (18 total)
varLabels and varMetadata:
Sample_Name: Sample_Name
Sample_Well: Sample_Well
...: ...
Sentrix_Position: Sentrix_Position
(7 total)
featureData
featureNames: GI_10047089-S, GI_10047091-S, ..., trpF (47293 total)
varLabels and varMetadata: none
experimentData: use 'experimentData(object)'
Annotation [1] "Illumina"
QC Information
Available Slots: Signal StDev Detection
featureNames: 1475542110_F, 1475542113_E, 1475542114_A, ..., 1475542113_D,
1475542113_F
sampleNames: Biotin, cy3_high, cy3_low, ..., pm, negative
=====================
When I extract the data files from Rdata.zip in the package's data directory and
run the following commands from the Example.R file, I have the following:
=====================
> targets = readBeadSummaryTargets("targets.txt")
> targets
DataFile SampleSheet QCInfo
1 raw_data.csv raw_data_sample_sheet.csv raw_data_qcinfo.csv
> BSData <- readBeadSummaryData(targets)
Reading file raw_data.csv
Reading file raw_data.csv
> BSData
ExpressionSetIllumina (storageMode: list)
assayData: 47293 features, 18 samples
element names: exprs, BeadStDev, NoBeads, Detection
phenoData
rowNames: IH.1, IC.1, ..., P42.2 (18 total)
varLabels and varMetadata:
Sample_Name: Sample_Name
Sample_Well: Sample_Well
...: ...
Sentrix_Position: Sentrix_Position
(7 total)
featureData
featureNames: 1, 2, ..., 47293 (47293 total)
varLabels and varMetadata: none
experimentData: use 'experimentData(object)'
Annotation character(0)
QC Information
Available Slots:
featureNames: 1475542110_F, 1475542113_E, 1475542114_A, ..., 1475542113_D,
1475542113_F
sampleNames: Biotin, cy3_high, cy3_low, ..., pm, negative
=====================
My problem is that I don't have the Target ID's in this second approach. In
particular note that featureNames in the second approach are not set to the
TargetIDs as they are in the BSData item.
Can you help me solve this problem,
cheers,
Keith
========================
Keith Satterley
Bioinformatics Division
The Walter and Eliza Hall Institute of Medical Research
Parkville, Melbourne,
Victoria, Australia
More information about the Bioconductor
mailing list