Brian Tsai
2011-Oct-16 14:51 UTC
[R] background normalization in rma() in the affy package
Hi, i was looking into the documentation for the rma() function in affy() package, and was trying to figure out how exactly the background normalization is done. I read all three papers cited in the rma() documentation, but the most detailed explanation i could find was in Irizary et al., 2003, where they state that they compute B(PM_{ijn}) = E[s_{ijn} | PM_{ijn}] where s_{ijn} is assumed to be exponential, and bg_{ijn} is normal. I still don't understand what value is being computed here, neither am i clear on what the correction looks like. i.e. if s_{ijn} is an exponentially-distributed random variable, how is bg_{ijn} fit into this? thanks! [[alternative HTML version deleted]]