Question

two-color microarray - technical replicates

0

Entering edit mode

Sanches • 0

@sanches-13853

Last seen 6.1 years ago

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?

biolrep<- c(1, 1, 2,2, 3) ## I am not sure about this part
design <- c(-1,1,-1,1) ## I am not sure about this part too
corfit<- duplicateCorrelation(MA, design, ndups = 1, block=biolrep)
fit <- lmFit(MA, design, block = biolrep, cor = corfit$consensus)
fit<-eBayes(fit)
topTable(fit, adjust = "BH")

Thanks in advance

cheers,

Sanches

microarray limma TARGET FILE • 1.3k views

ADD COMMENT • link 6.7 years ago Sanches • 0

0

Entering edit mode

Dear Gordon Smyth,

Thanks for your reply.

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.

Thanks in advance.

ADD REPLY • link updated 6.7 years ago by Gordon Smyth 50k • written 6.7 years ago by Sanches • 0

0

Entering edit mode

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.

Anyway, see my answer below.

ADD REPLY • link 6.7 years ago Gordon Smyth 50k

score 1 · Answer 1 · 2017-08-28

1

Entering edit mode

Gordon Smyth 50k

@gordon-smyth

Last seen 7 hours ago

WEHI, Melbourne, Australia

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.

ADD COMMENT • link 6.7 years ago Gordon Smyth 50k

0

Entering edit mode

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?

Thanks in advance

ADD REPLY • link 6.7 years ago Sanches • 0

0

Entering edit mode

Dear Gordon,

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

Thanks in advance

ADD REPLY • link 6.7 years ago Sanches • 0