**0**wrote:

Hi,

I am wondering which test is suitable for my case two condition i.e. C (3 replicates) and T (3 replicates)) and what is basically makes difference between "Wald" and "LRT" test implemented in DESeq2. I tested both with my dataset but i was surprised when i end up with much differences in terns of DE tags. Then i tested with example dataset and found there is really differences in both tests.

here i m presenting both tests using example of deseq and got differences in terms of DE tags:

# count tables from RNA-Seq data cnts <- matrix(rnbinom(n=1000, mu=100, size=1/0.5), ncol=10) cond <- factor(rep(1:2, each=5)) cond # object construction dds <- DESeqDataSetFromMatrix(cnts, DataFrame(cond), ~ cond) # standard analysis (Default Wald test) dds <- DESeq(dds) res <- results(dds) sum(na.omit(res$pvalue <= 0.05)) ###7 sum(na.omit(res$pvalue <= 0.05 & abs(res$log2FoldChange) >= 1)) ###4 sum(na.omit(res$pvalue <= 0.05 & abs(res$log2FoldChange) >= 1.5)) ###0 sum(na.omit(res$padj <= 0.05)) ###0 sum(na.omit(res$padj <= 0.05 & abs(res$log2FoldChange) >= 1)) ###0 sum(na.omit(res$padj <= 0.05 & abs(res$log2FoldChange) >= 1.5)) ###0 # an alternate analysis: likelihood ratio test ddsLRT <- DESeq(dds, test="LRT", reduced= ~ 1) resLRT <- results(ddsLRT) sum(na.omit(resLRT$pvalue <= 0.05)) ###7 sum(na.omit(resLRT$pvalue <= 0.05 & abs(resLRT$log2FoldChange) >= 1)) ###7 sum(na.omit(resLRT$pvalue <= 0.05 & abs(resLRT$log2FoldChange) >= 1.5)) ###2 sum(na.omit(resLRT$padj <= 0.05)) ###0 sum(na.omit(resLRT$pvalue <= 0.05 & abs(resLRT$log2FoldChange) >= 1)) ###7 sum(na.omit(resLRT$padj <= 0.05 & abs(resLRT$log2FoldChange) >= 1.5)) ###0

**I read the manual and understand the concept behind the wald (applicable when we have two condition) and LRT test (when we have more than 2 conditions and intent to check interaction between conditions). But what if i only have two condition like Male (15 or more replicate) and Female (15 replicates or more replicate).**

Then in this case what test would be more appropriate using DESeq2 ?

Thanks in advance for any advice

**25k**• written 3.3 years ago by unique379 •

**0**