Entering edit mode
yshuang
•
0
@yshuang-20457
Last seen 5.3 years ago
Hello
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
library(affy)
data <- ReadAffy()
eset <- rma(data)
write.exprs(eset, file="output.txt")
Environment and results:
settings | test1 | test2 |
OS | Mac OS X 10.6.8 Snow Leopard | CentOS 7.6 |
r | 2.13.1 | 3.5.3 |
affy | 1.30.0 | 1.60.0 |
sample (cel.gz) | GSM278159 | GSM278159 |
probe set | 200784_s_at | 200784_s_at |
rma result | 4.903961609 | 4.896936087 |