[BioC] hthgu133a + gcrma = pms error

Richard Finney rfinney5 at yahoo.com
Mon Jun 30 19:47:57 CEST 2008


Using a recent version of R, I get an error running gcrma on hthgu133a (high throughput Affymetrix U133A) chips. The error is 
    Error in pms[, i] <- GSB.adj(Yin = pms[, i], subset = index.affinities,  :
    number of items to replace is not a multiple of replacement length
    Calls: gcrma -> bg.adjust.gcrma -> gcrma.engine
    Execution halted
___
Here are the details ...

[root at cbiodev74 GCRMATEST]# uname -a
Linux cbiodev74.xxx.xxx.xxx 2.6.9-67.ELsmp #1 SMP Wed Nov 7 13:56:44 EST 2007 x86_64 x86_64 x86_64 GNU/Linux
[root at cbiodev74 GCRMATEST]# cat j2.R
library(affy)
library(hthgu133acdf)
library(hthgu133aprobe)
library(gcrma)
dat <- ReadAffy();
eset <- gcrma(dat)
write.exprs(eset,file="output.gcrma.txt")
bye("no")
[root at cbiodev74 GCRMATEST]# cat j2.R | R --vanilla

R version 2.7.1 (2008-06-23)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

<blah blah ... boilerplate startup messages removed to keep it short>
Type 'q()' to quit R.

> library(affy)
Loading required package: Biobase
Loading required package: tools

Welcome to Bioconductor

  Vignettes contain introductory material. To view, type
  'openVignette()'. To cite Bioconductor, see
  'citation("Biobase")' and for packages 'citation(pkgname)'.

Loading required package: affyio
Loading required package: preprocessCore
> library(hthgu133acdf)
> library(hthgu133aprobe)
Loading required package: matchprobes
> library(gcrma)
Loading required package: splines
> dat <- ReadAffy();
> eset <- gcrma(dat)
Adjusting for optical effect....Done.
Computing affinities.Done.
Adjusting for non-specific binding.Error in pms[, i] <- GSB.adj(Yin = pms[, i], subset = index.affinities,  :
  number of items to replace is not a multiple of replacement length
Calls: gcrma -> bg.adjust.gcrma -> gcrma.engine
Execution halted
[root at cbiodev74 GCRMATEST]#



More information about the Bioconductor mailing list