[BioC] backgroundCorrect offset value
Gordon K Smyth
smyth at wehi.EDU.AU
Tue Aug 30 01:02:30 CEST 2011
Dear Prasad,
> Date: Mon, 29 Aug 2011 02:17:01 +0000
> From: Prasad Siddavatam <siddavatam at gmail.com>
> To: <bioconductor at stat.math.ethz.ch>
> Subject: Re: [BioC] backgroundCorrect offset value
>
> Dear Dr. Smyth,
>
> Thank you very much for your detailed response. I am going to read the
> references. I understand why I am getting more differential genes with bigger
> offset value. Basically we are bringing the variance close to zero (with more
> uniformity across probes).
>
> Are you suggesting to delete the probes before doing the background correct?
No, I'm not. All probes should be retained for background correction.
Non-expressed probes should be filtered before using eBayes().
> If yes, is there any limit on the maximum number of probes to be deleted?
>
> In the variance stabilization, which maximum value (approximately) of
> fit$df.prior is considered high enough to call a good variance stabilization is
> achieved.
There is no maximum value. Higher is better.
I think you're worrying about this more than is necessary. A decent value
like offset=50 will give good results in a wide variety of situations.
You can even use
fit <- eBayes(fit, trend=TRUE)
which will makes uniformity of the variance less important. Again, use
plotSA(fit)
to see what this does.
Best wishes
Gordon
> Appreciate your time and help.
>
> Prasad
______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}
More information about the Bioconductor
mailing list