[BioC] arrayQualityMetrics problem with Agilent 2-colour

Edwin Groot edwin.groot at biologie.uni-freiburg.de
Wed Jun 23 12:22:16 CEST 2010


Hello all,
Would anyone know of a simple way to plot the heatmap in
arrayQualityMetrics on a per-channel basis, rather than the per-array
basis?

I have 16 2-colour Agilent arrays read in by the LIMMA package, but
arrayQualityMetrics shows only 16 arrays in the heatmap, rather than
the 32 red and green samples that I would expect. When one array is an
outlier in the heatmap, I do not know which sample (red or green) is at
fault.

The code that I used is quite simple:
> library(limma)
> library(arrayQualityMetrics)
> targets <- readTargets("qc.exp")
> RG <- read.maimages(targets$FileName, source="agilent")
> arrayQualityMetrics(expressionset=RG, outdir="default", force = TRUE,
do.logtransform = TRUE)
The directory 'default' has been created.
KernSmooth 2.23 loaded
Copyright M. P. Wand 1997-2009
(loaded the KernSmooth namespace)
[[1]]
[[2]]

I also attempted a custom call to aqm.* functions. However, the
preparatory step fails:
> RGprep = aqm.prepdata(expressionset=RG, do.logtransform = TRUE)
Error in function (classes, fdef, mtable)  : 
  unable to find an inherited method for function "aqm.prepdata", for
signature "RGList"

> sessionInfo()
R version 2.11.1 (2010-05-31) 
i486-pc-linux-gnu 

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=C              LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base
    

other attached packages:
[1] arrayQualityMetrics_2.6.0 affyPLM_1.24.0           
[3] preprocessCore_1.10.0     gcrma_2.20.0             
[5] affy_1.26.1               Biobase_2.8.0            
[7] limma_3.4.2              

loaded via a namespace (and not attached):
 [1] affyio_1.16.0        annotate_1.26.0      AnnotationDbi_1.10.1
 [4] beadarray_1.16.0     Biostrings_2.16.3    DBI_0.2-5           
 [7] genefilter_1.30.0    grid_2.11.1          hwriter_1.2         
[10] IRanges_1.6.4        lattice_0.18-8       latticeExtra_0.6-11 
[13] marray_1.26.0        RColorBrewer_1.0-2   RSQLite_0.9-1       
[16] simpleaffy_2.24.0    splines_2.11.1       stats4_2.11.1       
[19] survival_2.35-8      tools_2.11.1         vsn_3.16.0          
[22] xtable_1.5-6        

Thanks in advance,
Edwin
-- 
Dr. Edwin Groot, postdoctoral associate
AG Laux
Institut fuer Biologie III
Schaenzlestr. 1
79104 Freiburg, Deutschland
+49 761-2032945



More information about the Bioconductor mailing list