Voom without design matrix
2
0
Entering edit mode
Ninni Nahm ▴ 50
@ninni-nahm-6347
Last seen 9.9 years ago
Germany

Hi all!

I have some RNA-seq data with only one phenotype (disease)  and no controls. I want to use voom to 'convert' the data to 'microarray' like data. Is this possible? 

I do not want to do a differential expression analysis. 

Best 

Nina

limma voom • 3.4k views
ADD COMMENT
0
Entering edit mode

Hi! 

Thank you! I want to use elastic net with the output. My plan was to use the fpkm table with voom and then apply elastic net.

Nina

ADD REPLY
2
Entering edit mode
@james-w-macdonald-5106
Last seen 19 minutes ago
United States

Certainly. You don't ever need to specify a design matrix; voom() will by default just use a design matrix that only has an intercept term. From ?voom:

 design: 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.

 

ADD COMMENT
1
Entering edit mode
@gordon-smyth
Last seen 28 minutes ago
WEHI, Melbourne, Australia

Are you planning to analyze the output values using software tools that understand precision weights? If not, then you won't be fully using the voom output, so you may be better off using the cpm() or rpkm() functions of the edgeR package.

Edit: given your plan to use elastic net, I would use rpkm() with a large prior.count, say prior.count=5.

ADD COMMENT

Login before adding your answer.

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