[BioC] edgeR for entirely individual gene wise dispersion
Gordon K Smyth
smyth at wehi.EDU.AU
Mon Apr 29 01:45:03 CEST 2013
Dear Gu,
Yes, you can set prior.df=0 to get genewise dispersion estimates without
squeezing, although this not recommended except perhaps as a diagnostic.
In the formula G_0=20/df in the paper, df is the residual df and
prior.df=20. So
prior.df = G_0 * df.residual
Best wishes
Gordon
---------------------------------------------
Professor Gordon K Smyth,
Bioinformatics Division,
Walter and Eliza Hall Institute of Medical Research,
1G Royal Parade, Parkville, Vic 3052, Australia.
http://www.statsci.org/smyth
> Date: Sat, 27 Apr 2013 23:46:20 -0700
> From: Gu Mi <neo.migu at gmail.com>
> To: bioconductor at r-project.org
> Subject: [BioC] edgeR for entirely individual gene wise dispersion
>
> Dear All:
>
> I am using edgeR for RNA-Seq data analysis. According to the paper
> "Differential expression analysis of multifactor RNA-Seq experiments
> with respect to biological variation" by McCarthy, DJ et al. (2012), the
> tagwise dispersion model is recommended. It is a compromise between
> entirely individual genewise dispersion and the common dispersion
> models.
>
> My question is, is it possible to obtain a "pure genewise" model, i.e.
> there is NO shrinkage applied towards the common/trended dispersion? I
> think what I mean is to estimate \phi by maximizing APL_g (\phi) + G_0 *
> APL_s (\phi) where G_0 is set to zero (no shared part). Which argument
> in the estimateGLMTagwiseDisp function shall I change? Can I set
> prior.df = 0 to get a purely genewise model without any shrinkage? In
> the paper, G_0 = 20/df, but I am not sure if this "df" in the
> denominator is the prior.df argument in R, or something else.
>
> Thank you very much!
>
> Best,
> Gu
>
______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}
More information about the Bioconductor
mailing list