Question: Limma Design matrix for time series+treatment+batch effect
0
gravatar for senthil.duraikannukailasam
11 weeks ago by

dear all, I'm trying to get a DEG analysis of a microarray data.

I have 9 placebo and 9 treated mouse samples. The

3 replicates for each condition factors: 1. TimeOfTreatment <- "2days", " 7days", "12days" 2. TreatmentType <- "placebo" , "drug" 3. BatchID <- "001" , "002"

  TimeOfTreatment <- rep(c("2","4","6","2","4","6","2","4","6"),2)
  TreatmentType <- rep(c("placebo","drug"),9)
  BatchID  <- rep(c("1","1","2","1","1","1","2","1","1"),2)
  groupS <- paste(TimeOfTreatment,TreatmentType,sep= "_")

  model.matrix(~ 0 + groupS + Batch)

Now, I'm not sure how to get the contrast matrix!
For say , 1. I wanted to compare the DGE between all Treated and untreated samples. 2. Between Drug Treated 6th day with Drug Treated 4th day. 3. Between Drug Treated 4th day with Drug Treated 2nd day.

ADD COMMENTlink modified 11 weeks ago by Aaron Lun25k • written 11 weeks ago by senthil.duraikannukailasam0
Answer: Limma Design matrix for time series+treatment+batch effect
2
gravatar for Aaron Lun
11 weeks ago by
Aaron Lun25k
Cambridge, United Kingdom
Aaron Lun25k wrote:

Just use makeContrasts. You will have to fix your names so that they are syntactically valid, though.

I wanted to compare the DGE between all Treated and untreated samples.

Compare averages of the groupS levels, e.g., (A1 + A2 + A3)/3 - (B1 + B2 + B3)/3.

Between Drug Treated 6th day with Drug Treated 4th day. 3. Between Drug Treated 4th day with Drug Treated 2nd day.

Compare the relevant levels of groupS, e.g., Drug4 - Drug2, and so on.

ADD COMMENTlink modified 11 weeks ago • written 11 weeks ago by Aaron Lun25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 345 users visited in the last hour