Edger Dealing With Paired Experiment
1
0
Entering edit mode
8.2 years ago
liran0921 ▴ 110

Hi all,

I am using EdgeR for DE analysis. I would appreciate if somebody could help me check my experiment design.

I have 6 cows (2,63,66,67,73,74) for a treatment, each with 3 time points(m14,p7,p28). I want to compare between p28 vs m14, p7 vs m14 and p28 vs p7.

cow <- factor(c(2,2,2,63,63,63,66,66,66,67,67,67,73,73,73,74,74,74))
timep <- factor(c("m14","p7","p28","m14","p7","p28","m14","p7","p28","m14","p7","p28","m14","p7","p28","m14","p7","p28"))
data.frame(Sample=colnames(y),cow,timep)
design<-model.matrix(~cow+timep)

This follows the EdgerR user guide (31 March 2013 revised), case study 4.4 on page 53.

And here I have 18 libraries. So I would like to ask what if one library is missing? In this way, how to make the design? If I still use the way mentioned above, Edger cannot handle it.

Thanks.

edger paired • 2.1k views
ADD COMMENT
0
Entering edit mode
8.2 years ago

what if you try:

model.mat <- data.frame(Sample=colnames(y),cow,timep)
design <- model.matrix(~cow+timep, data=model.mat)

Addiitonally, you can try only using the paired samples (and ignoring the unpaired samples).

ADD COMMENT
0
Entering edit mode

Hi, thanks for your reply. The method you proposed get the same design as mine. So I suppose the matrix is right? Here it is:

 (Intercept) cow63 cow66 cow67 cow73 cow74 timepp28 timepp7
           1     0     0     0     0     0        0       0
           1     0     0     0     0     0        0       1
           1     0     0     0     0     0        1       0
           1     1     0     0     0     0        0       0
           1     1     0     0     0     0        0       1
           1     1     0     0     0     0        1       0
           1     0     1     0     0     0        0       0
           1     0     1     0     0     0        0       1
           1     0     1     0     0     0        1       0
           1     0     0     1     0     0        0       0
           1     0     0     1     0     0        0       1
           1     0     0     1     0     0        1       0
           1     0     0     0     1     0        0       0
           1     0     0     0     1     0        0       1
           1     0     0     0     1     0        1       0
           1     0     0     0     0     1        0       0
           1     0     0     0     0     1        0       1
           1     0     0     0     0     1        1       0
ADD REPLY

Login before adding your answer.

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