[BioC] different gcrma output on Linux and windows
Vladimir.Morozov at ari.altanapharma.com
Vladimir.Morozov at ari.altanapharma.com
Fri Nov 12 20:46:40 CET 2004
Hi
I have got different GCRMA results on Windows and Linux. Basicly I run
the standart GCRMA on three replicate U133A chips and get correlation 0.86
between the windows and linux results.
I have the same package versions on linux and windows:
> package.version("gcrma")
[1] "1.1.3"
> package.version("hgu133acdf")
[1] "1.4.3"
> package.version("hgu133aprobe")
[1] "1.0"
>
On windows
> R.Version()
$platform
[1] "i386-pc-mingw32"
.......
>data <-
ReadAffy(filenames=c(".............../GP00067.CEL","............/GP00025.CEL","....../GP00057.CEL"),phenoData="PC3.pdata")
> eset.gcrma <- gcrma(data)
Computing affinities.Done.
Adjusting for optical effect...Done.
Adjusting for non-specific binding...Done.
Normalizing
Calculating Expression
>write.table(exprs(eset.gcrma),file="PC3.eset.gcrma",quote=FALSE,sep="\t",
row.names = TRUE, col.names = TRUE)
copy PC3.eset.gcrma to PC3.eset.w.gcrma on Linux
on linux
> package.version("gcrma")
[1] "1.1.3"
>
> R.Version()
$platform
[1] "i686-pc-linux-gnu"
........
> data <-
ReadAffy(filenames=c("......./GP00067.CEL","..../GP00025.CEL","...../GP00057.CEL"),phenoData="PC3.pdata")
> eset.gcrma <- gcrma(data)
Computing affinities.Done.
Adjusting for optical effect...Done.
Adjusting for non-specific binding...Done.
Normalizing
Calculating Expression
#read the windows expression values
> eset.gcrm.w <- read.exprSet("PC3.eset.w.gcrma")
> cor(exprs(eset.gcrma),exprs(eset.gcrm.w))
GP00067.CEL GP00025.CEL GP00057.CEL
GP00067.CEL 0.8633871 0.8613008 0.8612103
GP00025.CEL 0.8629481 0.8615257 0.8608711
GP00057.CEL 0.8639332 0.8619907 0.8625088
#just correlation function check
> cor(exprs(eset.gcrma),exprs(eset.gcrma))
GP00067.CEL GP00025.CEL GP00057.CEL
GP00067.CEL 1.000000 0.999486 0.999442
GP00025.CEL 0.999486 1.000000 0.999133
GP00057.CEL 0.999442 0.999133 1.000000
Alspn u133B chip result are different.Thanks in advance for any
suggestions
Vladimir
[[alternative HTML version deleted]]
More information about the Bioconductor
mailing list