Question

Linear regression in DESeq2

0

Entering edit mode

chimeric • 0

@chimeric-9840

Last seen 6.6 years ago

Is it possible to conduct a linear regression analysis in DESeq2, where one could get a slope, r-squared, and p-value? Based on this post (DESEQ2 linear model of dose) it looks like I can get a slope (though it confuses me how log2FoldChange could be interpreted as a slope...clarification would be helpful here).

From my previous analysis in DESeq, I've identified which genes change over time relative to control using a LRT analysis. Now I would like to identify the pattern of gene expression over time, irrespective of the control. For example, did expression start low, and then end high? I had planned to do a simple linear regression analysis on the normalized counts, but was wondering if it would be more appropriate to do this in DESeq2. From the post above, it sounds like I would be able to tell the direction of the slope, but it is unclear to me if it would also have the information to tell the goodness of fit of the data, and whether the slope is significantly different from zero.

Thank you!

Erin

deseq • 3.3k views

ADD COMMENT • link updated 6.8 years ago by Michael Love 41k • written 6.8 years ago by chimeric • 0

Wolfgang Huber · Answer 1 · 2017-07-27

We have a FAQ about continuous covariates in the DESeq2 vignette. Briefly, yes, you can include them, and this is modeling the log of counts over changes in the covariates. The LFC is the change for each unit in the covariate. You can see this from the formula log(q) = X beta in the vignette and DESeq2 paper. Like any other covariate, you can obtain Wald or LRT results tables for this coefficient.