Question: two-color microarray - technical replicates
0
21 months ago by
Sanches0
Sanches0 wrote:

Hi all,

I am using Limma for the first time to analyze my two-color microarray data. So, I would like to know the best way to proceed to create my Targets file. My experiment has the following technical replicates:

CY3                                  CY5

1. gene of interest                       control

2. gene of interest                       control

3. control                                     gene of interest

4. control                                     gene of interest

5. gene of interest (treated)         control (treated)

6. gene of interest (treated)          control (treated)

7. control (treated)                        gene of interest(treated)

.  1=2 / 3=4 / 5=6

My question is: how can I proceed in this part below?

corfit<- duplicateCorrelation(MA, design, ndups = 1, block=biolrep)
fit <- lmFit(MA, design, block = biolrep, cor = corfit\$consensus)
fit<-eBayes(fit)

cheers,

Sanches

microarray limma target file • 570 views
modified 21 months ago • written 21 months ago by Sanches0

Dear Gordon Smyth,

So, in this case, do I have to write my script the following way below?

targets <- cbind(Cy3=c("Gene of interest","Gene of interest", "Control", "Control", "Gene of interest_Treatment", "Gene of interest_Treatment", "Control_Treatment"),  Cy5=c("Control", "Control", "Gene of interest", "Gene of interest", "Control_Treatment",

"Control_Treatment", "Gene of interest_Treatment"))

rownames(targets) <- paste("Array",1:7)

parameters <- cbind(C=c(-1,-1, 0, 0, 1, 1, 2), T=c(0, 0, -1, -1, 2 , 2, 1))
rownames(parameters) <- c("Reference","Control","Treatment")  #  What is the right way to write this part?

modelMatrix(targets, parameters)

modelMatrix(targets, ref="Reference")

design<- model.matrix(targets, parameters)

# OBS.: I ran the script but it appears an error.

No, that isn't right, for a few different reasons. If you use the 'parameters' argument to modelMatrix() then it has to refer to samples that appear in your experiment. You can't just make up names like "Reference" that don't actually appear in your targets file.

Answer: two-color microarray - technical replicates
1
21 months ago by
Gordon Smyth37k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth37k wrote:

You have 7 arrays, so the design matrix has to have 7 rows and biolrep has to have 7 entries.

Why don't you use modelMatrix(), which makes the design matrix for you?

Edit:

You ask how to make the targets file. The targets file just contains the two columns Cy3 and Cy5 that you already posted as part of your question. So you seem to have made the targets file ok already.

I see now that your experiment is a little more complicated in that you have two controls rather than one. modelMatrix() can make the design matrix for you, but it's probably easier to do it manually:

design <- cbind(Dye=1, GeneVsCtrl.Untreated=c(-1,-1,1,1,0,0,0), GeneVs.Ctrl.Treated=c(0,0,0,0,-1,-1,1))

With this design matriix, the second coefficient compares your gene of interest to the control without treatment, while the second coefficient commpares your gene of interest to the control with treatment. The first column is a good idea to account for gene-specific dye effects.

Hi Gordon,

Actually, I would like to compare:  Gene Vs Control, Gene Vs Control(treated), Gene Vs Gene(treated) and, Gene(treated) Vs Control (Treated). So I have been having trouble designing the matrix.

Problem 1: my target file is composed by two block corresponding to biological replicates (Array1=Array2 / Array 3=Array4 / Array5=Array6).  Problem 2: There are also technical replicates in dye-swap pairs: Array12 vs Array34 / Array 56 vs Array7.

Question: Is it necessary use "duplicateCorrelation function? Is there another fast way to do that?

Dear Gordon,

I think I can use separate channel analysis to analyze my data.  What is your opinion about that?