Limma help (creating design matrix)
1
0
Entering edit mode
d808bc07 • 0
@d808bc07
Last seen 3.2 years ago

Hi, I'm trying to use limma to get a list of differentially expressed genes between subtypes (GCB and ABC) of B-cell lymphoma. I have a dataset from an Affymetrix array with gene expression reads for 271 samples. 71 of which were classified as ABC, 144 as GCB and 56 were unclassified.

I've managed to read my CEL files and perform normalisation on the gene expression data:

norm.batch = rma(batch)
dat = exprs(norm.batch)

But I'm just not sure where to start with creating a design matrix for the limma analysis

limma • 494 views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 1 hour ago
United States

There's no need to do the second line in your example. The limma package is quite happy to accept an ExpressionSet as input, so you don't need a matrix. In addition, you could hypothetically include gene information in your ExpressionSet which will then be included in your output.

Anyway, have you read the limma User's Guide? Section 9.3 seems to be exactly what you are after.

ADD COMMENT

Login before adding your answer.

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