[BioC] problem with function rma

Juan Lin julin at aecom.yu.edu
Wed May 30 02:09:11 CEST 2007


Hello Bioconductorers,

I am working with 12 affy rat2302 chips recently. It was OK when I tried to
read my cel files into R:

> RatData<-ReadAffy()
> list.celfiles()
 [1] "Rat 16hrs-1-2113.CEL"  "Rat 16hrs-2-2114.CEL"  "Rat 16hrs-3-2115.CEL"
 [4] "Rat 30hrs-1-2116.CEL"  "Rat 30hrs-2-2117.CEL"  "Rat 30hrs-3-2118.CEL"
 [7] "Rat 4hrs-1-2110.CEL"   "Rat 4hrs-2-2111.CEL"   "Rat 4hrs-3-2112.CEL"
[10] "Rat Normal-1-2107.CEL" "Rat Normal-2-2108.CEL" "Rat Normal-3-2109.CEL"

> phenoData(RatData) <- read.phenoData("phenoData.txt",header=T)
> show(pData(phenoData(RatData)))
                        ID time replicate
Rat Normal-1-2107.CEL 2107    0         1
Rat Normal-2-2108.CEL 2108    0         2
Rat Normal-3-2109.CEL 2109    0         3
Rat 4hrs-1-2110.CEL   2110    4         1
Rat 4hrs-2-2111.CEL   2111    4         2
Rat 4hrs-3-2112.CEL   2112    4         3
Rat 16hrs-1-2113.CEL  2113   16         1
Rat 16hrs-2-2114.CEL  2114   16         2
Rat 16hrs-3-2115.CEL  2115   16         3
Rat 30hrs-1-2116.CEL  2116   30         1
Rat 30hrs-2-2117.CEL  2117   30         2
Rat 30hrs-3-2118.CEL  2118   30         3

However, when I was going to do normalization, I got an error message as
following:

> esetRat <- rma(RatData)
Background correcting
Normalizing
Calculating Expression
Error in validObject(.Object) : invalid class "exprSet" object: sampleNames
different from names of phenoData rows

Can anyboday give me a clue on how to solve this problem?

BTW. my sessionInfor is as following:
> sessionInfo()
R version 2.4.0 (2006-10-03)
i386-pc-mingw32

locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

attached base packages:
[1] "splines"   "tools"     "methods"   "stats"     "graphics"  "grDevices"
[7] "utils"     "datasets"  "base"

other attached packages:
  rat2302cdf affyQCReport   simpleaffy   genefilter     survival
affy
    "1.14.0"     "1.12.0"      "2.8.0"     "1.12.0"       "2.31"
"1.12.1"
      affyio      Biobase
     "1.2.0"     "1.12.2"

Thanks a lot!
Juan



More information about the Bioconductor mailing list