[BioC] Create affyBatch from mouse exon array data using ReadAffy or extractAffyBatch() from aroma.affymetrix
De Bondt, An-7114 [PRDBE]
ADBONDT at PRDBE.jnj.com
Thu Jul 17 09:33:06 CEST 2008
Thanks, Christian and Henrik, for your feedback!
With respect to setting up the AffymetrixCelSet in aroma.affymetrix, I used checkChipType=FALSE because of the use of the alternative chipType (MmEx10stv1_Mm_ENSE in stead of MoEx-1_0-st-v1). If I use checkChipType=TRUE in this setup, I get the following:
Error in list("AffymetrixCelSet$byName(projectName, chipType = chipType, checkChipType = T" = <environment>, :
[2008-07-17 08:25:21] Exception: Invalid name of directory containing CEL files. The name of the directory (MmEx10stv1_Mm_ENSE) must be the same as the chip type used for the CEL files (MoEx-1_0-st-v1) unless using argument 'checkChipType=FALSE': rawData/myDataSet/MmEx10stv1_Mm_ENSE
at throw.default("Invalid name of directory containing CEL files. The name of
at throw("Invalid name of directory containing CEL files. The name of the dire
at fromFiles.AffymetrixCelSet(static, path = path, cdf = cdf, ...)
at fromFiles(static, path = path, cdf = cdf, ...)
at withCallingHandlers(expr, warning = function(w) invokeRestart("muffleWarnin
at method(static, ...)
at AffymetrixCelSet$byName(projectName, chipType = chipType, checkChipType = T
With respect to xps, your reference to script4xps.R is really helpful. Do you have a URL from where the files "MoEx-1_0-st-v1.r2.clf", "MoEx-1_0-st-v1.r2.pgf", "MoEx-1_0-st-v1.na25.mm9.probeset.csv", "MoEx-1_0-st-v1.na25.mm9.transcript.csv" can be downloaded? I searched on the Affy site but did not find it, sorry.
Does xps work with AffyBatches? If not, is it possible to create an AffyBatch with the raw data?
From: henrik.bengtsson at gmail.com [mailto:henrik.bengtsson at gmail.com]On
Behalf Of Henrik Bengtsson
Sent: Wednesday, 16 July 2008 21:48
Cc: De Bondt, An-7114 [PRDBE]; bioconductor at stat.math.ethz.ch
Subject: Re: [BioC] Create affyBatch from mouse exon array data using
ReadAffy or extractAffyBatch() from aroma.affymetrix
I never received the original message for this one - was it posted to
BioC? Anyway, my comments below.
On Wed, Jul 16, 2008 at 11:59 AM, cstrato <cstrato at aon.at> wrote:
> Dear An
> I cannot answer your question regarding aroma.affymetrix,
> but since you also mention "xps":
> Please note that xps can handle mouse exon arrays,
> see the file "script4xps.R" in directory examples of how
> to import the necessary clf, pgf and annotation files.
> Please let me know if you experience any problems.
> Best regards
> C.h.i.s.t.i.a.n S.t.r.a.t.o.w.a
> V.i.e.n.n.a A.u.s.t.r.i.a
> e.m.a.i.l: cstrato at aon.at
> De Bondt, An-7114 [PRDBE] wrote:
>> Dear UseRs,
>> I am analysing a dataset from mouse exon arrays using aroma.affymetrix. I
>> can read the raw data using following code.
>> chipType <- "MmEx10stv1_Mm_ENSE"
>> cdf <- AffymetrixCdfFile$fromChipType(chipType = chipType)
>> # setup the CEL set; read the raw data
>> projectName <- "myDataSet"
>> cs <- AffymetrixCelSet$byName(projectName, chipType = chipType,
>> checkChipType=FALSE, cdf = cdf)
Actually, that is not doing anything but setting up the
AffymetrixCelSet. It does not read in the data, except validating
that the CEL files are consistent with each other (and the CDF). BTW,
you should only use 'checkChipType=FALSE', if you really know what you
are doing; if it gives an error otherwise, there is often a good
reason for it.
>> Next, I would like to make an AffyBatch from these raw data but I stumble
>> at a memory message (see below). This is the same message as directly
>> affyBatchRaw <- ReadAffy(filenames =
>> paste("./rawData/myDataSet/MmEx10stv1_Mm_ENSE/", celfiles, sep = ""))
I see that you previously/below tried:
affyBatchRaw <- extractAffyBatch(cs)
which is pretty much the same as the above. In aroma.affymetrix we
use prefix "extract..." on method names to make it explicit that you
load all data into memory and that any changes done on the obtained
object will *not* be reflected in the underlying data files.
Having said all this, what you do above is not really utilizing the
aroma.affymetrix package at all. All your problems are unrelated to
that package and has to do with the 'affy' package.
Your alternative is to do your exon analysis in 'aroma.affymetrix'
(see online Vignettes), or use 'xps' as Christian suggests.
>> I am working on a linux machine with 70GB of memory.
>> Did anyone experience this before? Is this typical for mouse exon arrays?
>> I tried using exonmap as well as xps but, as far as I experienced, they are
>> not yet adjusted for mouse exon arrays.
>> Thanks in advance for your help!
>> Kind regards,
>>> affyBatchRaw <- extractAffyBatch(cs)
>> Error in read.affybatch(filenames = l$filenames, phenoData = l$phenoData,
>> : Calloc could not allocate (-1889533886 of 48) memory
>> 5: .Call("read_abatch", filenames, rm.mask, rm.outliers, rm.extra,
>> ref.cdfName, dim.intensity, verbose, PACKAGE = "affyio")
>> 4: read.affybatch(filenames = l$filenames, phenoData = l$phenoData,
>> description = l$description, notes = notes, compress = compress, rm.mask
>> = rm.mask, rm.outliers = rm.outliers, rm.extra = rm.extra, verbose =
>> verbose, sd = sd, cdfname = cdfname)
>> 3: ReadAffy(filenames = filenames, sampleNames = sampleNames, ...,
>> verbose = as.logical(verbose))
>> 2: extractAffyBatch.AffymetrixCelSet(cs)
>> 1: extractAffyBatch(cs)
>> R version 2.6.2 (2008-02-08) x86_64-unknown-linux-gnu
>> attached base packages:
>>  tools stats graphics grDevices utils datasets methods
>>  base
>> other attached packages:
>>  mmex10stv1mmensecdf_10.0.0 farms_1.3.1 
>> MASS_7.2-42 preprocessCore_1.0.0  affyio_1.6.1
>> Biobase_1.16.3  aroma.affymetrix_0.9.3
>> aroma.apd_0.1.3  R.huge_0.1.5 affy_1.16.0
>>  affxparser_1.10.2 aroma.core_0.9.3 
>> sfit_0.1.5 aroma.light_1.8.1  digest_0.3.1
>> matrixStats_0.1.2  R.rsp_0.3.4
>> R.cache_0.1.7  R.utils_1.0.2 R.oo_1.4.3
>>  R.methodsS3_1.0.1
>> loaded via a namespace (and not attached):
>>  rcompgen_0.1-17
>> [[alternative HTML version deleted]]
>> Bioconductor mailing list
>> Bioconductor at stat.math.ethz.ch
>> Search the archives:
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> Search the archives:
More information about the Bioconductor