Main public logs
Appearance
Combined display of all available logs of 太極. You can narrow down the view by selecting a log type, the username (case-sensitive), or the affected page (also case-sensitive).
- 12:44, 20 May 2025 Brb talk contribs created page File:Combat DE mod.png (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...)
- 12:44, 20 May 2025 Brb talk contribs uploaded File:Combat DE mod.png (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...)