[BioC] EdgeR: coefficients for GLM for multiple groups

James W. MacDonald jmacdon at uw.edu
Thu Mar 21 21:27:26 CET 2013


Hi Capricy,

I'm not convinced you have read the user's guide carefully. There is, 
for example, a section called 'Quick start' that is only seven pages in, 
that shows the canonical glm analysis steps. Then there are two more 
worked examples starting on page 52 that show pretty much all the steps 
required to do a glm analysis.

However, in the code you show below you haven't followed the examples, 
nor have you used the accessor function that is intended to produce 
useful output.

Best,

Jim



On 3/21/2013 3:57 PM, capricy gao wrote:
> I think I have carefully read both the userguide and reference manual, 
> but couldn't figure it out. Could anybody help me out what arguments 
> should be included? Thanks a lot.
>
> In limma fit results: there is something like this:
> ===============
> $p.value
>             Grp1      Grp2vs1
> Gene 1 0.8604469 6.019156e-05
> Gene 2 0.2174605 1.673262e-05
> Gene 3 0.3571957 5.758422e-02
> Gene 4 0.6789641 1.758690e-01
> Gene 5 0.4589329 8.223679e-01
> 95 more rows ...
> ==============
>
> But in edgeR I couldn't find it as I posted previously....
>
>
> ------------------------------------------------------------------------
> *From:* James W. MacDonald <jmacdon at uw.edu>
> *To:* capricy gao <capricyg at yahoo.com>
> *Cc:* "bioconductor at r-project.org" <bioconductor at r-project.org>
> *Sent:* Thursday, March 21, 2013 1:45 PM
> *Subject:* Re: [BioC] EdgeR: coefficients for GLM for multiple groups
>
> Hi Capricy,
>
> Just like limma, edgeR has a very comprehensive user's guide that you
> can consult.
>
> library(edgeR)
> edgeRUsersGuide()
>
> Best,
>
> Jim
>
>
>
> On 3/21/2013 2:41 PM, capricy gao wrote:
> > I am looking at my fit results from GLM of edgeR for my data:
> >
> > and assume that $coefficients are related to the expression levels. 
> Just wondering how I can get the p values for those coefficients since 
> I remember in limma, each coefficient will be accompanied by the 
> corresponding t value and p value.
> >
> >
> > Thanks a lot for the help.
> >
> > capricy
> >
> > ===================
> >
> >> design=model.matrix(~0+regroup)
> >> y<- estimateGLMCommonDisp(y,design)
> >> y<- estimateGLMTagwiseDisp(y,design)
> >> fit<- glmFit(y,design)
> >> fit
> > An object of class "DGEGLM"
> > $coefficients
> >            regroupFR  regroupFA  regroupFM  regroupFP regroupFW  
> regroupMA
> > GS_14929  -9.543565  -9.267612  -9.310823  -9.167299 -10.62094  
> -9.260763
> > GS_09776 -10.304430 -10.985146 -10.644154 -10.640469 -10.59615 
> -10.798889
> > GS_18434 -11.327664 -11.786421 -11.643292 -11.902728 -12.50470 
> -11.654008
> > GS_08334 -10.789181 -10.271089 -10.511480 -10.375642 -10.39865 
> -10.836647
> > GS_09550 -10.564167 -10.152571 -10.410428 -10.098825 -10.36302  
> -9.722892
> >            regroupMM  regroupMP  regroupMW
> > GS_14929  -8.516256  -8.222225  -9.495394
> > GS_09776 -10.337434 -10.284703 -10.588544
> > GS_18434 -11.755905 -11.662276 -11.583950
> > GS_08334 -10.797529 -10.643409 -10.783104
> > GS_09550 -10.650495 -10.766462  -9.880359
> > 15457 more rows ...
> >
> > $fitted.values
> >                  MP        FA      FR        FW        MA      FW.1  
>       MM
> > GS_14929 596.53153 205.75696 84.61834 31.844025 169.27178  82.62258 
> 513.34590
> > GS_09776  75.81638  36.91137 39.53083 32.643715  36.33823  84.69746  
> 83.04882
> > GS_18434  19.09812  16.54813 14.19827  4.825192  15.43876  12.51946  
> 20.07922
> > GS_08334  52.95485  75.41272 24.33895 39.775089  34.99088 103.20054  
> 52.41066
> > GS_09550  46.82022  84.90529 30.48442 41.218519 106.62251 106.94567  
> 60.71721
> >                  FP      MA.1        MW        FM      MM.1      
> MW.1      MP.1
> > GS_14929 78.226172 123.52247 193.95899 122.61523 173.571103 
> 194.05252 568.23926
> > GS_09776 17.921685  26.51705  64.98434  32.30797  28.080239  
> 65.01567  72.22056
> > GS_18434  5.064548  11.26610  23.99487  11.88420  6.789129  
> 24.00644  18.19234
> > GS_08334 23.358730  25.53385  53.48880  36.89427  17.720948  
> 53.51460  50.44331
> > GS_09550 30.811447  77.80550 131.97315  40.81919  20.529537 
> 132.03679  44.59963
> >                FM.1      FP.1      FA.1      FR.1
> > GS_14929 259.44483 324.86875 256.36216 198.47462
> > GS_09776  68.36129  74.42772  45.98959  92.72063
> > GS_18434  25.14609  21.03277  20.61808  33.30242
> > GS_08334  78.06556  97.00745  93.96021  57.08766
> > GS_09550  86.37041 127.95815 105.78744  71.50204
> > 15457 more rows ...
> >
> > $counts
> >            MP  FA  FR FW  MA FW.1  MM  FP MA.1  MW  FM MM.1 MW.1 
> MP.1 FM.1 FP.1
> > GS_14929 221 284 170 23 267  105 209 106  52 218 158  277  170  926  
> 185  211
> > GS_09776  75  32  17 28  41  96  86  15  23  65  23  27  65  73  87  85
> > GS_18434  36  22  19  8  21    6  22  7    7  27  13    6  21    2  
> 23  15
> > GS_08334  44  77  23 50  41  78  78  14  21  61  29    8  46  59  
> 94  132
> > GS_09550  82  92  45 54 105  75  95  18  79 153  41    8  111  11  
> 86  178
> >          FA.1 FR.1
> > GS_14929  159    0
> > GS_09776  52  143
> > GS_18434  14  23
> > GS_08334  92  60
> > GS_09550  97  39
> > 15457 more rows ...
> >
> > $deviance
> >  GS_14929  GS_09776  GS_18434  GS_08334  GS_09550
> > 19.059229  3.629342  9.509356  4.268291  8.809585
> > 15457 more elements ...
> >
> > $df.residual
> > [1] 9 9 9 9 9
> > 15457 more elements ...
> >
> > $abundance
> > [1]  -9.081235 -10.552328 -11.717230 -10.579804 -10.233947
> > 15457 more elements ...
> >
> > $design
> >    regroupFR regroupFA regroupFM regroupFP regroupFW regroupMA regroupMM
> > 1        0        0        0        0        0        0        0
> > 2        0        1        0        0        0        0        0
> > 3        1        0        0        0        0        0        0
> > 4        0        0        0        0        1        0        0
> > 5        0        0        0        0        0        1        0
> >    regroupMP regroupMW
> > 1        1        0
> > 2        0        0
> > 3        0        0
> > 4        0        0
> > 5        0        0
> > 13 more rows ...
> >
> > $offset
> >          [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]
> > [1,] 14.61341 14.59445 13.9819 14.08233 14.39241 15.03576 14.75727 
> 13.52703
> > [2,] 14.61341 14.59445 13.9819 14.08233 14.39241 15.03576 14.75727 
> 13.52703
> > [3,] 14.61341 14.59445 13.9819 14.08233 14.39241 15.03576 14.75727 
> 13.52703
> > [4,] 14.61341 14.59445 13.9819 14.08233 14.39241 15.03576 14.75727 
> 13.52703
> > [5,] 14.61341 14.59445 13.9819 14.08233 14.39241 15.03576 14.75727 
> 13.52703
> >          [,9]    [,10]    [,11]    [,12]  [,13]    [,14]    [,15]    
> [,16]
> > [1,] 14.07733 14.76322 14.12002 13.67291 14.7637 14.56482 14.86951 
> 14.95085
> > [2,] 14.07733 14.76322 14.12002 13.67291 14.7637 14.56482 14.86951 
> 14.95085
> > [3,] 14.07733 14.76322 14.12002 13.67291 14.7637 14.56482 14.86951 
> 14.95085
> > [4,] 14.07733 14.76322 14.12002 13.67291 14.7637 14.56482 14.86951 
> 14.95085
> > [5,] 14.07733 14.76322 14.12002 13.67291 14.7637 14.56482 14.86951 
> 14.95085
> >          [,17]    [,18]
> > [1,] 14.81434 14.83441
> > [2,] 14.81434 14.83441
> > [3,] 14.81434 14.83441
> > [4,] 14.81434 14.83441
> > [5,] 14.81434 14.83441
> > 15457 more rows ...
> >
> > $dispersion
> > [1] 0.7464955 0.2734189 0.4034924 0.2832934 0.3819788
> > 15457 more elements ...
> >
> > $method
> > [1] "oneway"
> >
> > $samples
> >    group lib.size norm.factors
> > MP    MP  2220863            1
> > FA    FA  2179157            1
> > FR    FR  1181036            1
> > FW    FW  1305802            1
> > MA    MA  1780507            1
> > 13 more rows ...
> >     [[alternative HTML version deleted]]
> >
> >
> >
> > _______________________________________________
> > Bioconductor mailing list
> > Bioconductor at r-project.org <mailto:Bioconductor at r-project.org>
> > https://stat.ethz.ch/mailman/listinfo/bioconductor
> > Search the archives: 
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
> -- 
> James W. MacDonald, M.S.
> Biostatistician
> University of Washington
> Environmental and Occupational Health Sciences
> 4225 Roosevelt Way NE, # 100
> Seattle WA 98105-6099
>
>
>

-- 
James W. MacDonald, M.S.
Biostatistician
University of Washington
Environmental and Occupational Health Sciences
4225 Roosevelt Way NE, # 100
Seattle WA 98105-6099



More information about the Bioconductor mailing list