[BioC] Different RMA expression values in different BioC versions?
Quanli Wang
quanli at stat.duke.edu
Wed Jan 10 19:51:04 CET 2007
Thank you all for the quick responses. Now seems that I figured out the
problem. It has nothing todo with BioC but seems to me that it has
something to do with the different implementations of "merge" function in
R2.3.1 and R2.4.0. When I tried to reorder the expression values using
merge function from R based on one prefered probeid list with "sort"
option set to FALSE,R2.3.1 gives an order based on the x, while R2.4.0
gives an order based on y.
Here is the little legacy script that does the ordering:
eset <- read.table("eset.RData")
geneID <- geneNames(eset)
genedescriptorfile = "C:\\CAGP\\9999\\cel\\h133+.genedescriptors"
genedes <- read.table(genedescriptorfile,sep="\t",header=FALSE, row.names=NULL)
rmamatrix2 <- cbind(geneID,rmamatrix)
cprmamatrix <-
merge(x=rmamatrix2,y=genedes,by.x=1,by.y=1,all.x=FALSE,all.y=FALSE,sort=FALSE)
rmamatrix <- cprmamatrix[,2:ncol(rmamatrix2)]
Can someone verify that it is indeed the case even it is more like an R
problem, not a BioC problem.
Best, Quanli
On
Wed, 10 Jan 2007, James W. MacDonald wrote:
> Quanli Wang wrote:
>>> Different than what? An earlier version with the same chips?
>>
>> Yes. For the same chips using justRMA, the values from R2.3.1/BioC1.8 will
>> be very different from that from R2.4.0/BioC1.9.
>> Quanli
>
> I don't see that here. I get identical results for the same chips using these
> two versions of R/BioC
>
>> eset <- justRMA()
> Background correcting
> Normalizing
> Calculating Expression
>> save(eset, file="eset2.3.1.Rdata")
>> sessionInfo()
> Version 2.3.1 (2006-06-01)
> i386-pc-mingw32
>
> attached base packages:
> [1] "tools" "methods" "stats" "graphics" "grDevices" "utils"
> "datasets"
> [8] "base"
>
> other attached packages:
> hgu133plus2cdf affy affyio Biobase
> "1.12.0" "1.10.0" "1.0.0" "1.10.1"
>
> [close R-2.3.1, start R-2.4.0]
>
>> eset2 <- justRMA()
> Background correcting
> Normalizing
> Calculating Expression
>> load("eset2.3.1.Rdata")
>> all.equal(exprs(eset), exprs(eset2))
> [1] TRUE
>> sessionInfo()
> R version 2.4.1 (2006-12-18)
> 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] "tools" "stats" "graphics" "grDevices" "utils" "datasets"
> "methods"
> [8] "base"
>
> other attached packages:
> hgu133plus2cdf affy affyio Biobase
> "1.14.0" "1.12.2" "1.2.0" "1.12.2"
>
> Can you show the code you used to get different results, and the output of
> sessionInfo(), as requested in my last email?
>
> Best,
>
> Jim
>
>
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>
>
More information about the Bioconductor
mailing list