Complex Limma design: technical replication, biological replication and repeated experiment
0
0
Entering edit mode
Guest User ★ 13k
@guest-user-4897
Last seen 9.6 years ago
Dear limma-users, I'm currently trying to create the correct design matrix for a limma analysis of my NimbleGen two-color microarray data. Basically, I have a control treatment (Ctr) and 5 Cy treatments (Cy). For each comparison Cy-Ctr, I have 2 biological replicates per treatment that are dye swapped on the array. These biological reps are also in duplicate (so each sample is on two arrays). For each comparison, I have also repeated the exposure experiment twice (i.e. the actual experiment was replicated, not the array, so new biological samples). So I actually have three "levels of replication": 1. technical replicates 2. biological replicates within a single experiment 3. replication of the entire experiment. I want to compare all my controls to each individual Cy treatment without the need to know differences between replicates. (i.e.I want to know the genes that are differentially expressed in all my samples of Cy treatment 1 for example versus all controls). I am now struggling how to make the correct design matrix with these levels of replication as I have several blocking factors. This is how my targets file looks like right now: SlideNumber SampleCy3 SampleCy5 Cy3 Cy5 TechnicalReps BiolocialReps 1 1.1 1.3 Ctr Cy1 1 1 2 1.1 1.3 Ctr Cy1 1 1 3 1.4 1.2 Cy1 Ctr 2 -1 4 1.4 1.2 Cy1 Ctr 2 -1 5 2.1 2.3 Ctr Cy1 3 2 6 2.1 2.3 Ctr Cy1 3 2 7 2.4 2.2 Cy1 Ctr 4 -2 10 2.4 2.2 Cy1 Ctr 4 -2 11 3.1 3.3 Ctr Cy2 5 3 12 3.1 3.3 Ctr Cy2 5 3 13 3.4 3.2 Cy2 Ctr 6 -3 14 3.4 3.2 Cy2 Ctr 6 -3 15 4.1 4.3 Ctr Cy2 7 4 16 4.1 4.3 Ctr Cy2 7 4 17 4.4 4.2 Cy2 Ctr 8 -4 18 4.4 4.2 Cy2 Ctr 8 -4 19 5.1 5.3 Ctr Cy3 9 5 20 5.1 5.3 Ctr Cy3 9 5 21 5.4 5.2 Cy3 Ctr 10 -5 22 5.4 5.2 Cy3 Ctr 10 -5 23 6.1 6.3 Ctr Cy3 11 6 24 6.1 6.3 Ctr Cy3 11 6 25 6.4 6.2 Cy3 Ctr 12 -6 26 6.4 6.2 Cy3 Ctr 12 -6 27 7.1 7.3 Ctr Cy4 13 7 28 7.1 7.3 Ctr Cy4 13 7 29 7.4 7.2 Cy4 Ctr 14 -7 30 7.4 7.2 Cy4 Ctr 14 -7 31 8.1 8.3 Ctr Cy4 15 8 32 8.1 8.3 Ctr Cy4 15 8 33 8.4 8.2 Cy4 Ctr 16 -8 34 8.4 8.2 Cy4 Ctr 16 -8 35 9.1 9.3 Ctr Cy5 17 9 36 9.1 9.3 Ctr Cy5 17 9 37 9.4 9.2 Cy5 Ctr 18 -9 38 9.4 9.2 Cy5 Ctr 18 -9 39 10.1 10.3 Ctr Cy5 19 10 40 10.1 10.3 Ctr Cy5 19 10 41 10.4 10.2 Cy5 Ctr 20 -10 42 10.4 10.2 Cy5 Ctr 20 -10 What would be the correct formula in limma to analysis this experiment? right now i'm just doing design <- modelMatrix(targets,ref="Ctr") fit <- lmFit(MAprobes.n,design) cont.matrix<-makeContrasts(Cy_1=Cy1, Cy_2=Cy2, Cy_3=Cy3, Cy_4=Cy4, Cy_5=Cy5, Cy_6=Cy6,levels=design) c.fit<-contrasts.fit(fit, contrasts=cont.matrix) but this does not take the technical replication into account nor the fact that the biological replicates are not on the "same level" (same experiment vs replicated experiment). I know the use of blocking factor and how to estimate correlation, but I really don't know how to integrate these three levels all together. Thanks for the help. -- output of sessionInfo(): no output -- Sent via the guest posting facility at bioconductor.org.
Microarray limma Microarray limma • 864 views
ADD COMMENT

Login before adding your answer.

Traffic: 582 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