# Setup model matrix
design <- model.matrix(~ 0 + factor(rep(1:2, each = 3)))
colnames(design) <- paste0("ng", c(50, 100))
yeast %>%
# Normalize data
psrn("identifier") %>%
# Get mean-variance trends
calculate_mean_sd_trends(design)
Run the code above in your browser using DataLab