de_out1 <- differential_expression_wrapper(expr_median = expr_median_sample_cluster, 
  md = md, model = "lm", formula = formula_lm, K = K)

apply(de_out1$adjp < FDR_cutoff, 2, table)

adjp_BCRXLvsRef
FALSE              51
TRUE               42

de_out2 <- differential_expression_wrapper(expr_median = expr_median_sample_cluster, 
  md = md, model = "lmer", formula = formula_lmer, K = K)

apply(de_out2$adjp < FDR_cutoff, 2, table)

adjp_BCRXLvsRef
FALSE              23
TRUE               70