[BioC] results of justPlier vs. RMA
He, Yiwen (NIH/CIT) [C]
heyiwen at mail.nih.gov
Thu Nov 9 22:37:57 CET 2006
When comparing the results of RMA from the affy package and justPlier
from the plier package, I found that the summarized data are clustered
very differently. The cluster of the RMA results is very close to MAS5
values, showing 5 groups. However, PLIER results show a very different
and strange clustering pattern, almost like using single linkage. (See
I used Pearson correlation coefficient distance and average linkage for
I remember attending a web talk about PLIER and the speaker mentioned
post-processing scaling in PLIER. Is that done in justPlier? Or is this
an additional step that needs to follow justPlier? Or is it something
Thank you for your help!
Here is my code:
R version 2.4.0 (2006-10-03)
attached base packages:
 "tools" "methods" "stats" "graphics" "grDevices" "utils"
 "datasets" "base"
other attached packages:
plier affy affyio Biobase
"1.4.0" "1.12.0" "1.2.0" "1.12.2"
> eset1 <- rma(myData)
> exprs(eset1) <- 2^exprs(eset1)
> eset2 <- justPlier(myData, normalize=T)
> exprs(eset2) <- 2^exprs(eset2)
More information about the Bioconductor