[BioC] affy: adding new probe set definitions and making rma expression summaries
Mette Langaas
mettela at math.ntnu.no
Fri Jul 8 01:29:24 CEST 2005
1) Problem: want to add five probe sets (from existing probes) in
hgfocus-array and make expression summaries for all probesets (existing
+ five new) using rma (medianpolish)
2) Have: AffyBatch object, dsnorm, of normalized probes.
ds <- ReadAffy(filenames=paste(celcatalog,celfilenames,sep=""))
dsbgcorr <- bg.correct(ds,method="rma")
dsnorm <- normalize(dsbgcorr,method="quantiles")
3) Have: definition of the new probesets:
hlag <- "211528_x_at"
locpm <- unlist(pmindex(ds,hlag))
locmm <- unlist(mmindex(ds,hlag))
# the five new (overlapping...) probsets by redefining the original
probesets
hlaglocation <- list(HLAG1to7=cbind(pm=locpm[1:7],mm=locmm[1:7]),
HLAG8=cbind(pm=locpm[8],mm=locmm[8]),
HLAG9=cbind(pm=locpm[9],mm=locmm[9]),
HLAG10to11=cbind(pm=locpm[10:11],mm=locmm[10:11]),
HLAG9to11=cbind(pm=locpm[9:11],mm=locmm[9:11]))
4) Howto: how is it possible to "computeExprSet" (or related) with
method "medianpolish" for the new probesets located as given in
hlaglocation together with all the existing probesets in the AffyBatch
object dsnorm?
I have read through Bioconductor documentation and mail-list-questions
about affy, altcdfenvs, and sadly I am not able to understand
how to do this :-(
I'm running
> version
_
platform i686-redhat-linux-gnu
arch i686
os linux-gnu
system i686, linux-gnu
status
major 2
minor 1.1
year 2005
month 06
day 20
language R
Your help is greatly appreciated!
Best regards
Mette
--
Mette Langaas
Department of Mathematical Sciences
Norwegian University of Science and Technology
NO-7491 Trondheim, Norway
Voice: +47-7359-1697 URL : http://www.math.ntnu.no/~mettela
Fax : +47-7359-3524 Email: mette.langaas at math.ntnu.no
--
Mette Langaas
Department of Mathematical Sciences
Norwegian University of Science and Technology
NO-7491 Trondheim, Norway
Voice: +47-7359-1697 URL : http://www.math.ntnu.no/~mettela
Fax : +47-7359-3524 Email: mette.langaas at math.ntnu.no
More information about the Bioconductor
mailing list