Based on this post Question: Differential gene expression based on read counts using DESeq package.
Using htseq_count, I generated the raw read count for each gene for my cancer and normal sample. Then I used DESeq and DESeq2 package to identify differential gene expression.
1) First 4 rows are from DESeq output
2) Next 4 rows are from DESeq 2 output
3) Then followed by HTseq_count output
Just to understand, I manually computed the log2(fold_change) for both DESeq and DESeq2 R package output.
For DESeq output, I first divided (value2/value1) and then took log(value2/value1,2). I received the same value as computed by DESeq.
But for DESeq2 output, I divided (val2/val1) and then took log(val2/val1,2). The DESeq2 computed log2 fold change is different from mine. Am I missing something?