File:Combat DE mod.png: Revision history

20 May 2025

• curprev 12:4412:44, 20 May 2025 Brb talk contribs 1,670 bytes +1,670 Demonstrate the importance of the 'mod' parameter in sva::ComBat() <syntaxhighlight lang='r'> library(sva) library(limma) library(ggplot2) library(tidyr) set.seed(42) n_genes <- 100 n_samples <- 20 # Group: 10 Case and 10 Control — spread across batches group <- sample(rep(c("Control", "Case"), each=10)) # now randomized batch <- rep(c(1, 2), each=10) # Simulate expression data expr <- matrix(rnorm(n_genes * n_samples, mean=5, sd=1), nrow=n_genes) # Add true signal to 10 genes for Case...