Entering edit mode
3.9 years ago
Argonaut
•
0
Hi,
I am using DESeq2 Likelihood Ratio Test for multi-group comparison (6 different genotypes).
###Example
##### dds <- DESeqDataSetFromMatrix(countData = data, colData = meta, design = ~ genotype)
### Likelihood ratio test
dds_lrt <- DESeq(dds, test="LRT", reduced = ~ 1)
As far as I know LRT doesn't use log2foldchange to calculate the p.value. Thus, I am interested to know how DESeq2 calculates these fold changes ( 6 genotypes in my case), and is there any way to interpret them.
Thanks! PS. I would appreciate any link to the tutorial/paper on how to perform multi-group RNA-seq analysis and interpret it.