[BioC] edgeR

Gordon K Smyth smyth at wehi.EDU.AU
Thu Apr 10 01:35:09 CEST 2014


Aaron,

One reason for your surprise may be that you have read a lot of things 
into my reply that weren't there.

I simply showed how to use to software to collate results from multiple 
contrasts.  The contrasts could be anything including interactions.

Gordon

On Tue, 8 Apr 2014, Aaron Mackey wrote:

> I'm surprised that the given answer here was to check whether two fold
> changes are themselves significant or not, rather than what I tend to think
> of as a "contrast of contrasts" scenario to test whether the two contrasts'
> fold changes are significantly different from another between groups (i.e.
> testing for an interaction term)?  A transcript that fails to pass the
> threshold for dt1 may be insignificant at that particular threshold, yet
> still indistinguishable from the dt2 response that does happen to pass the
> threshold.  The absence of significance is not evidence for non-response.
> The OP seemed to be asking for a statistical test to look for differences
> in response between the two groups.  Another problem with !dt1 & dt2 is
> that dt1 and dt2 may both be "true" but of opposite signs (though I see
> this addressed in a later email).
>
> So in a nutshell: is there statistical reason to prefer independent FDR
> threshold tests vs. direct testing for an interaction term?
>
> Thanks,
> -Aaron
>
>
> On Mon, Apr 7, 2014 at 8:22 PM, Gordon K Smyth <smyth at wehi.edu.au> wrote:
>
>> If you want to select transcript that are DE for one contrast but not
>> another, first test each contrast:
>>
>>   lrt1 <- glmLRT(fit, contrast=mycontrast1)
>>   lrt2 <- glmLRT(fit, contrast=mycontrast2)
>>
>> Then apply significance thresholds:
>>
>>   dt1 <- decideTestsDGE(lrt1)
>>   dt2 <- decideTestsDGE(lrt2)
>>
>> Then select the transcripts you want:
>>
>>   selected <- !dt1 & dt2
>>
>> Best wishes
>> Gordon
>>
>>
>>  Date: Sun,  6 Apr 2014 03:32:49 -0700 (PDT)
>>> From: "Jahn Davik [guest]" <guest at bioconductor.org>
>>> To: bioconductor at r-project.org, jahn.davik at bioforsk.no
>>> Subject: [BioC] edgeR
>>>
>>>
>>> Hi there,
>>>
>>
>>  I have a question regarding edgeR - or it might actually be a more
>>> general statistical question. In any case, I am using edgeR to analyse my
>>> read counts and really would appreciate help.
>>>
>>
>>  My experimental setup is:
>>>
>>
>>  Two genotypes (B and S)
>>> Two treatments ('trt' vs 'ntrt')
>>> Two time points (0hs 8hs).
>>> (Three bio reps)
>>>
>>> Now, I would like to identify reads that are specific to either of the
>>> genotypes as their response to the treatment over the time points.
>>
>>  I expect that I can do pairwise comparisons like:
>>
>>  'B_tr_0hs' vs 'B_trt_8hs'), and 'B_ntr_0hs' vs 'B_ntrt_8hs'), and
>>> continuing doing the same with the S-genotype. Subsequently, using a
>>> suitable tool, I could filter out the transcripts for, say, B's response to
>>> treatment over these two time points that are not found in B. It is,
>>> however, a little tedious so my question here is whether this can be
>>> modeled and extracted in edgeR's GLM ?
>>>
>>> regards
>>> JD

______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}



More information about the Bioconductor mailing list