Dear Donald,
Yes, the contrast to test the interaction effect should be DH-DC-CH+CC
((combo-control)-(x-control)-(y-control)).
If you don't get significant probes it means that the effect of the
combination of the two drugs is the sum of the marginal effects; that
could make sense depending on the nature of your experiment.
Hope it helps,
Oscar
On 29/3/11 16:27, "Donald Dunbar" <donald.dunbar@ed.ac.uk> wrote:
Hello,
I have a gene expression experiment with four replicates each for four
groups: control, two drugs and their combination.
Control (C.C)
Drug1 (D.C)
Drug2 (C.H)
Drug1+Drug2 (D.H)
I have imported the Affymetrix data, normalised and set up a design to
analyse in Limma.
C.C D.C C.H D.H
1 1 0 0 0
2 1 0 0 0
3 1 0 0 0
4 1 0 0 0
5 0 1 0 0
6 0 1 0 0
7 0 1 0 0
8 0 1 0 0
9 0 0 1 0
10 0 0 1 0
11 0 0 1 0
12 0 0 1 0
13 0 0 0 1
14 0 0 0 1
15 0 0 0 1
16 0 0 0 1
I initially set up contrasts and following lmFit, contrasts.fit,
eBayes,
and decideTests find differentially expressed probe-sets:
drug 1 v control contrast 'D.C-C.C' gives 46 probesets
drug 2 v control contrast 'C.H-C.C' gives 233 probesets
drug 1 + drug 2 combination contrast 'D.H-C.C' gives 3897 probesets
Clearly there are several thousand probe-sets differentially expressed
by the drug combination and not by the individual drugs. I can find
these crudely by filtering the p-values of the contrasts, but I would
rather do it properly.
I would like to define an interaction parameter and select
differentially expressed genes. I read Gordon Smyth's Limma user guide
and several helpful posts on this list but have so far failed to find
an
appropriate parametrization. The user guide's factorial example has
sepatare controls for each treatment group so isn't applicable for me
I
think. James MacDonald gave a helpful reply to Suzanne Szak in May
2009
and I thought I'd try that:
(combo-control)-(x-control)-(y-control) but I think I'm
misunderstanding
his answer: if I try this (D.H-C.C)-(D.C-C.C)-(C.H-C.C) I get no
significant probe-sets.
My own idea was (D.H-D.C)-(C.H-C.C) but that also gave no significant
probe-sets and I think was a bit silly.
Any help on this would be greatly appreciated.
With best wishes,
Donald
--
Donald Dunbar
University of Edinburgh
donald.dunbar@ed.ac.uk
Oscar M. Rueda, PhD.
Postdoctoral Research Fellow, Breast Cancer Functional Genomics.
Cancer Research UK Cambridge Research Institute.
Li Ka Shing Centre, Robinson Way.
Cambridge CB2 0RE
England
This communication is from Cancer Research UK. Our website is at
www.cancerresearchuk.org. We are a registered charity in England and
Wales (1089464) and in Scotland (SC041666) and a company limited by
guarantee registered in England and Wales under number 4325234. Our
registered address is Angel Building, 407 St John Street, London, EC1V
4AD. Our central telephone number is 020 7242 0200.
This communication and any attachments contain information which is
confidential and may also be privileged. It is for the exclusive use
of the intended recipient(s). If you are not the intended
recipient(s) please note that any form of disclosure, distribution,
copying or use of this communication or the information in it or in
any attachments is strictly prohibited and may be unlawful. If you
have received this communication in error, please notify the sender
and delete the email and destroy any copies of it.
E-mail communications cannot be guaranteed to be secure or error free,
as information could be intercepted, corrupted, amended, lost,
destroyed, arrive late or incomplete, or contain viruses. We do not
accept liability for any such matters or their consequences. Anyone
who communicates with us by e-mail is taken to accept the risks in
doing so.
[[alternative HTML version deleted]]
