[BioC] question about affymetrix array:ht-hg-u133acdf

Richard [guest] guest at bioconductor.org
Fri Oct 14 13:03:47 CEST 2011


hello,everyone!
I found a question when processing the affymetrix microarray data that I downloaded from NCBI. I had suffered from this question for two dayes, and I couldn't find useful information from the internet.

the affymetrix array i processed was ht-hg-u133acdf. And the R code was:
-----------------------------------------
#!/usr/bin/env Rscript
library("methods")
library("Biobase")
library("affy")
library("genefilter")
library("gcrma")
library("simpleaffy")
library("xtable")
library("preprocessCore")
library("affyPLM")
library("RColorBrewer")
library("lattice")
library("affyQCReport")
library("hthgu133acdf")
work.data<-ReadAffy("GSM637758.cel","GSM637759.cel","GSM637800.cel")
QCReport(work.data,file="t.pdf")
------------------------------------------

Then R give me thus information(Sorry some sentences were translated):
------------------------------------------
error at plot(qc(object)) : 
  something was wrong of function 'plot'when evaluating parameter 'x',at  setCEnvironment(cdfn) : 
  Could not find array definition file ' hthgu133acdf.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
------------------------------------------

I got the information of the variable "work.data" in R:
----------------------------------------
 AffyBatch object
size of arrays=744x744 features (17 kb)
cdf=HT_HG-U133A (22277 affyids)
number of samples=3
number of genes=22277
annotation=hthgu133a
------------------------------------------

I tried to solve this problem by all means:
changing the qc environment:setQCEnvironment("hthgu133acdf");
remove the "hthgu133acdf" package, then reinstall it;
updating the version of R from 2.10 to 2.13.2, then reinstall all the packages needed;
searching the google, a few questions that were similar were found, but no answers suitable.
when I changed the array data to "hgu133acdf" or "hgu95av2", everything was OK! So i don't know if there is some problem with that package,or anyone could help me solve this question!

thanks

 -- output of sessionInfo(): 

R version 2.13.2 (2011-09-30)
Platform: x86_64-pc-linux-gnu (64-bit)

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

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

--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list