Skip to main content

Table 2 Summary of model diagnostics and differential expression for each method when applied to the real RNA-Seq data from the case study

From: lmerSeq: an R package for analyzing transformed RNA-Seq data with linear mixed effects models

 

lmerSeq

DREAM

rmRNAseq

Singular fits

200 (1.5%)

96 (0.7%)*

Not available

Levene’s test for heteroskedsticity

   

 p vaule \(<0.05\)

206 (1.6%)

149 (1.1%)

Not available

 BH adjusted p value \(<0.05\)

0

0

Not available

Kolmogorov–Smirnov test for normality

   

 p value \(<0.05\)

24 (0.2%)

17 (0.1%)

Not available

 BH adjusted p value \(<0.05\)

0

0

Not available

Number of DEGs

   

 48 h versus baseline

1

0

0

 1 week versus baseline

1452

1828

142

 1 week versus 48 h

362

476

182

  1. BH Benjamini–Hochberg, DEGs differentially expressed genes
  2. *Models had to be refit with lmerSeq using weights since DREAM does not report or exclude singular fits