Question: How to do within-sample gene comparisons with limma-voom?
gravatar for samey
11 months ago by
samey10 wrote:

I have an experiment where I want to compare expression levels between homologous genes within the same sample. This has been done before using DESeq2 ( The idea is to make a count matrix where each homologous gene pair is placed side by side on the same row, and then use an offsets matrix incorporating both gene lengths and library sizes.

E.g. if I had 4 samples and 1000 gene pairs then the count matrix would have dimensions 1000 x 8, the design is then ~ 1 + homolog + sampleId (where levels(homolog) is c(0, 1)). This approach seems to work well with DESeq2.

However, since I have many samples DESeq2 is becoming too slow, so I wanted to try the (in my experience faster) limma-voom workflow instead. But I do not find any way to incorporate gene lengths into the limma-voom analysis. Gene lengths of homologous gene pairs often differ slightly, and not taking this into account leads to many false positives.

Is there any way to input offsets into limma-voom? Or to somehow account for both sequencing depth and gene length differences?

If my question is not clear, feel free to ask for clarification.



voom rna-seq limma • 329 views
ADD COMMENTlink written 11 months ago by samey10
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: 954 users visited in the last hour