[BioC] arrayQualityMetrics error and question
David Westergaard
david at harsk.dk
Sun Mar 4 11:02:52 CET 2012
Hello,
I am trying to use the BioC package arrayQualityMetrics to asses the
quality of a range of experiments. I am, however, running into errors,
when processing Affymetrix Hugene designs, particularly
http://www.ebi.ac.uk/arrayexpress/experiments/E-GEOD-33443
I'm using the following code:
library('arrayQualityMetrics'); library('affy')
mixture.batch <- ReadAffy()
arrayQualityMetrics(expressionset=mixture.batch,
outdir=paste("/home/projects/rgroups/cheminf/dwe/figures/self_processed/QA/",veggie[i],sep=""),
do.logtransform=TRUE,
force=TRUE
)
And get the error:
Error in setQCEnvironment(cdfn) :
Could not find array definition file ' hugene10stv1cdf.qcdef '.
Simpleaffy does not know the QC parameters for this array type.
See the package vignette for details about how to specify QC
parameters manually.
Error in density.default(as.matrix(log2(mm(expressionset))), ...) :
'x' contains missing values
In addition: Warning message:
In aqm.qcstats(expressionset) :
'plot(qcStats)' from the package 'simpleaffy' failed for this dataset.
Error in obj[[i]]$legend : $ operator is invalid for atomic vectors
In addition: Warning message:
In arrayQualityMetrics(expressionset = mixture.batch, outdir =
paste("/home/projects/rgroups/cheminf/dwe/figures/self_processed/QA/",
:
Cannot draw the Perfect Match versus MisMatch plot
Now, from what I've read in an older mail to this list, Hugene is a
PM-only type design, building my own .qcdef wouldn't really do no
good. Googling around didn't really give any useful information. So
how do I get by this problem?
Also, why do only some array types output a complete index.html,
providing the numbers for outlier detection?
Best Regards,
David Westergaard
More information about the Bioconductor
mailing list