How to connect htseq-count to DESeq2?

Question

[takoyaki]

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.

[Nai]

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

Answer 1

[h.mon]

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.

[ATpoint]

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.