[BioC] Pulling significant gene ids from beadarray data
Johnstone, Alice
Alice.Johnstone at esr.cri.nz
Fri Aug 10 04:01:49 CEST 2007
I have been analysing my Illumina data from rat ref12 using beadarray
and the summary data. I have managed to get a list of significantly
different probes, but am at a loss as to how to get the gene info this
corresponds to. The ExpressionSetIllumina object doesn't contain any
annotation data, neither the ProbeID or the TargetID appear to be
incorporated following the code in the vignette (march 26 2007). Is
there a way to include this to ensure I am following the right probe?
Then is there a way to get the gene information directly to R, such as
when using the lumi package and nuID etc?
Any help gratefully appreciated!
Alice Johnstone
PhD Candidate
Institute of Environmental Science and Research Ltd
Kenepuru Science Centre
34 Kenepuru Drive
PO Box 50-348
Porirua
New Zealand
Tel: + 64 4 914-0717
Fax: + 64 4 914-0770
dataFile="ROD-AJ-BeadStudioV2_Sample_Probe_Profile.csv"
qcFile="ControlGeneProfile.csv"
sampleSheet="sampleSheet.csv"
BSData<-readBeadSummaryData(dataFile,qcFile=qcFile,sampleSheet=sampleShe
et,
header=TRUE,skip=7,columns=list(exprs="AVG_Signal",BeadStDev="BEAD_STDEV
",
NoBeads="Avg_NBEADS"),qc.columns=list(exprs="AVG.Signal",Detection="Dete
ction"),
qc.skip=0,qc.sep=",",sep=",")
> exprs(BSData)[1:10,1:2]
SAL SAL.1
1 151.8 151.2
2 142.6 133.4
3 143.6 151.0
4 135.6 138.2
5 149.5 153.6
6 133.1 155.1
7 141.2 157.4
8 144.2 135.4
9 168.0 132.7
10 126.5 144.7
> sessionInfo()
R version 2.5.1 (2007-06-27)
i386-pc-mingw32
locale:
LC_COLLATE=English_New Zealand.1252;LC_CTYPE=English_New
Zealand.1252;LC_MONETARY=English_New
Zealand.1252;LC_NUMERIC=C;LC_TIME=English_New Zealand.1252
attached base packages:
[1] "grid" "tools" "stats" "graphics" "grDevices" "utils"
"datasets" "methods" "base"
other attached packages:
beadarray beadarraySNP quantsmooth lodplot quantreg
SparseM affy affyio geneplotter
"1.4.0" "1.2.0" "1.2.0" "1.1" "4.08"
"0.73" "1.14.2" "1.4.1" "1.14.0"
lattice annotate Biobase limma
"0.15-11" "1.14.1" "1.14.1" "2.10.5"
More information about the Bioconductor
mailing list