Question: How to create "design" for RNA-seq experiment using limma?
gravatar for Parham
3.8 years ago by
Parham1.4k wrote:


I am using limma to analyse RNA-seq data. I have a problem I cannot sort it out by myself. This might be a very beginner question. I cannot figure out how to create design for voom . Its in chapter 15 of limma manual. What I created is: 

> dge
An object of class "DGEList"
            Rep1.NoB1 Rep1.B1 Rep2.NoB1 Rep2.B1 Rep3.NoB1 Rep3.B1
SPAC212.11          1       0         1       0         3       0
SPAC212.10          1       0         0       2         0       0
SPAC212.09c         0       0         0       0         1       0
SPNCRNA.70          2       1         4       1         3       1
SPAC212.08c       216     274       229     312       267     435
7010 more rows ...

          group lib.size norm.factors
Rep1.NoB1     1 12610138    0.9719185
Rep1.B1       2  9903882    1.0414141
Rep2.NoB1     1  8344326    0.9787473
Rep2.B1       2 11509548    1.0291455
Rep3.NoB1     1 13078338    0.9643794
Rep3.B1       2 18742285    1.0170714

The experiment design is NoB1 is control and B1 is the treatment and it is done in three replicates. So Rep1.NoB1, Rep2.NoB1 and Rep3.NoB1 are replicates of control and the other three are treatments.

Next step is:

v <- voom(dge,design,plot=TRUE)

That I can't make it work! My question is what do I need to make the design matrix of? According to voom documentation it says "design matrix with rows corresponding to samples and columns to coefficients to be estimated. Defaults to the unit vector meaning that samples are treated as replicates."  

What is this "coefficients"? And how can I group replicates together? 


Thanks a lot in advance, I appreciate your help!


rna-seq limma design • 1.7k views
ADD COMMENTlink modified 18 months ago by sneha108ss10 • written 3.8 years ago by Parham1.4k
gravatar for sneha108ss
18 months ago by
sneha108ss10 wrote:


I am also facing a similar problem. How did you create the design matrix??

ADD COMMENTlink written 18 months ago by sneha108ss10
Please log in to add an answer.


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