[BioC] combineAffyBatch & GCRMA?
Michael Seewald
mseewald at gmx.de
Thu Jun 17 09:43:04 CEST 2004
Dear Robert, dear Bioconductor people,
I would like to use combineAffyBatch to merge a dataset of Affymetrix
HG-U133Plus2 and HG-U133A arrays. The main point here is the comparison of the
two arrays - and possibly generation of a unified dataset that is valid for
further analysis. (Disposing half of the Plus2 information is not of concern
for now..)
So far I used the following to combine and normalize with RMA (works great!):
# DATA.133A, DATA.133Plus are generated by ReadAffy()
> DATA.combine <- combineAffyBatch( list(DATA.133A,DATA.133Plus),
+ c("hgu133aprobe","hgu133plus2probe"),
+ newcdf="newcdfenv", verbose=T
+ )
> newcdfenv <- DATA.combine$cdf
> DATA.A.RMA <- rma(DATA.combine$dat)
Then I run into trouble because (I believe) I have to work with the cdf
environment:
> DATA.A.GCRMA <- gcrma(DATA.A)
Computing affinitiesError in library(cdfpackagename, character.only = TRUE) :
There is no package called 'newcdfenvcdf'
I tried also:
> ai <- compute.affinities(cdfname="newcdfenv",verbose=T)
Error in library(cdfpackagename, character.only = TRUE) :
There is no package called 'newcdfenvcdf'
How can I create the correct "newcdfenvcdf"? And why doesn't
specifying "newcdfenv" not help?
Best wishes,
Michael
More information about the Bioconductor
mailing list