microarray agilent data
1
0
Entering edit mode
evoke • 0
@evoke-8479
Last seen 8.7 years ago
European Union

Hi,

I have a question about my agilent microarray data.

It includes  before & after infection data of the same individuals and also replicates of non-infected other individuals as a kind of  control to these infected ones . There are 15 individuals and thus 30 microarray files. I am using "limma" package to read and normalize them. However, i have doubts if i should use the codes below to identifiy genes differently expressed, because when i cluster the datasets, the replicates do not group together:

 

RG_2= backgroundCorrect(z, method="minimum")
MA.exp= normalizeBetweenArrays(RG_2)
MA.exp_log= log(MA.exp$E)

SibShip = factor(targets_paired$SibShip)
Treat = factor(targets_paired$Treatment, levels=c("Not Inf.","Infected"))
design = model.matrix(~SibShip+Treat)
fit = lmFit(MA.exp_log, design)
fit = eBayes(fit)
topTable(fit, coef="TreatInfected")

In this case, should I analyze the data by just grouping as Infected & Non-infected?

Thank you in advance.

 

microarray agilent limma • 1.2k views
ADD COMMENT
0
Entering edit mode
arfranco ▴ 130
@arfranco-8341
Last seen 9 months ago
European Union

According limma, you need to go through one within and between normalization if using two color arrays, or a single normalization if using a single color. You need to provide with this information, and also with the targets definition

ADD COMMENT
0
Entering edit mode

One more thing, Why are you doing the MA.exp_log step ?

Have you read that in the limma reference document?. I don't think you need to do it, and if a log is calculated, it should be in base 2

ADD REPLY

Login before adding your answer.

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