This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to connect htseq-count to DESeq2?

I want to know how to use htseq-count data for detecting DEGs by DESeq2.

We can get gene expression count matrix by htseq-count, and also we can use this output in DESeq2. However, the last 4 rows of htseq-count matrix are not count data but information about counting result.

Do everyone remove last 4 rows by themselves in R ?

Or are there any functions in DESeq2 to pass the raw result of htseq-count ?

So, in conclusion, I want to know how to connect htseq-count to DESeq2. Thank you.

r rna-seq next-gen deg mapping

HI,

I have two separate matrix for normal and control. DeSeq2fromHTSeqcount , may i have to run separately for each matrix. I want GE.txt for further analysis. Please help me to input the data in Deseq2

1 answer

Or are there any functions in DESeq2 to pass the raw result of htseq-count ?

Yes, there is a function called DESeqDataSetFromHTSeqCount(). You can read about it in a section of the manual titled htseq-count input.

scheme4193, DESeq2 exists for quite a while now, the manual is outstandingly comprehensive. Please read it carefully, you'll find answers there, starting from input preparations over design formulas etc and troubleshooting.

Log in to answer this question.