Question: Different versions of r/affy combination may lead to differences in rma normalization results
gravatar for yshuang
3 months ago by
yshuang0 wrote:

Our group revisited some old codes for microarray data analysis in our group. We re-ran the code in a system with old versions of r and affy and a system with the up-to-date versions. We are aware of some differences in results of rma normalization when using different versions of r, affy package (and Bioconductor may be also involved in). Some of the differences are noticed in the first decimal place. We would like to share our findings for your reference even though we have no idea how and why it happened. The following is the details of the tests.

Test data: GSE11001

#rma normalization
data <- ReadAffy()
eset <- rma(data)
write.exprs(eset, file="output.txt")

Environment and results:

settingstest1 test2
OSMac OS X 10.6.8 Snow Leopard CentOS 7.6
r2.13.1 3.5.3
affy1.30.0 1.60.0
sample (cel.gz)GSM278159 GSM278159
probe set200784_s_at 200784_s_at
rma result4.903961609 4.896936087
ADD COMMENTlink modified 3 months ago • written 3 months ago by yshuang0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 284 users visited in the last hour