lumi - limma contrast design
1
0
Entering edit mode
@paolo-kunderfranco-5158
Last seen 7.2 years ago
Dear All, Could you please help me to design a contrast for Differential Analysis? .... constrast.matrix <- makeContrasts ( --------, levels=design) ... I have 4 groups. A - normal B - cancer C - normal + drug D - cancer + drug I would like to find genes modulated in D versus C. These genes should be not affected or change significantly in B versus A Any suggestion? Many thanks, Paolo [[alternative HTML version deleted]]
Cancer Cancer • 1.2k views
ADD COMMENT
0
Entering edit mode
Axel Klenk ★ 1.0k
@axel-klenk-3224
Last seen 34 minutes ago
UPF, Barcelona, Spain
Dear Paolo, sounds like you want: (D - C) - (B - A) Cheers, - axel On Tue, Apr 15, 2014 at 9:16 AM, Paolo Kunderfranco < paolo.kunderfranco@gmail.com> wrote: > Dear All, > > Could you please help me to design a contrast for Differential Analysis? > > .... > constrast.matrix <- makeContrasts ( > --------, > levels=design) > ... > > > I have 4 groups. > > A - normal > B - cancer > C - normal + drug > D - cancer + drug > > I would like to find genes modulated in D versus C. These genes should be > not affected or change significantly in B versus A > > > Any suggestion? > > Many thanks, > > Paolo > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor@r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > -- Axel Klenk Research Informatician Information Management Drug Discovery Actelion Pharmaceuticals Ltd. * Gewerbestrasse 16 * CH-4123 Allschwil * Switzerland G12.O1.R10 VOIP 92 63 67 * phone +41 61 565 63 67 axel.klenk@actelion.com * www.actelion.com Address for visitors: Hegenheimermattweg 91 -- The information of this email and in any file transmitted with it is strictly confidential and may be legally privileged. It is intended solely for the addressee. If you are not the intended recipient, any copying, distribution or any other use of this email is prohibited and may be unlawful. In such case, you should please notify the sender immediately and destroy this email. The content of this email is not legally binding unless confirmed by letter. Any views expressed in this message are those of the individual sender, except where the message states otherwise and the sender is authorised to state them to be the views of the sender's company. For further information about Actelion please see our website at http://www.actelion.com [[alternative HTML version deleted]]
ADD COMMENT
0
Entering edit mode
I don't think that contrast answers the question as asked. The interaction contrast finds genes whose normal vs cancer ratio is affected by the drug. But Paolo is asking for genes that are significant in the C vs D contrast but not in the A vs B contrast. The solution, as has been discussed several times recently on this list, is to conduct two separate tests and do a set difference operation on them. -Ryan On Tue 15 Apr 2014 02:26:03 AM PDT, Axel Klenk wrote: > Dear Paolo, > > sounds like you want: > (D - C) - (B - A) > > Cheers, > > - axel > > > > On Tue, Apr 15, 2014 at 9:16 AM, Paolo Kunderfranco < > paolo.kunderfranco at gmail.com> wrote: > >> Dear All, >> >> Could you please help me to design a contrast for Differential Analysis? >> >> .... >> constrast.matrix <- makeContrasts ( >> --------, >> levels=design) >> ... >> >> >> I have 4 groups. >> >> A - normal >> B - cancer >> C - normal + drug >> D - cancer + drug >> >> I would like to find genes modulated in D versus C. These genes should be >> not affected or change significantly in B versus A >> >> >> Any suggestion? >> >> Many thanks, >> >> Paolo >> >> [[alternative HTML version deleted]] >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor at r-project.org >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: >> http://news.gmane.org/gmane.science.biology.informatics.conductor >> > > >
ADD REPLY

Login before adding your answer.

Traffic: 642 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6