[BioC] Removing genes from linear model fit advice sought
Yong Li
yong.li at zbsa.uni-freiburg.de
Thu Feb 11 10:00:43 CET 2010
Dear Stephen,
I think a previous post by Gordon Smyth talking about filtering should
answer most of your questions. You can find it at:
https://stat.ethz.ch/pipermail/bioconductor/2009-January/025827.html
Hope it helps.
Yong
stephen sefick wrote:
> I have done this with all of the genes in a microarray experiment
> including blanks, negative controls, empties, hk genes, and spike ins,
> genes of interest
> 1. RG <- backgroundCorrect(RG, method="normexp", offset=50)
>
> 2. MA <- maNorm(as(RG, "marrayRaw"), norm="twoD")
>
> 3. WA <- normalizeBetweenArrays(as(MA, "MAList"), method="scale")
>
> This seems sensible to me. Is it?
>
> I am now thinking that I should remove everything except for the know
> differentially expressed genes and the genes of interest before
> fitting the linear model, contrasts, bayesian smotthing. Is this a
> sensible coarse of action? Thanks for all of your help in advance.
> kind regards,
>
More information about the Bioconductor
mailing list