# Produce a design matrix
design <- model.matrix(~ 0 + factor(rep(1:2, each = 3)))
colnames(design) <- paste0("ng", c(50, 100))
# Normalize and log transform the data
yeast_norm <- psrn(yeast, "identifier")
# Generate the plot
yeast_norm %>%
calculate_mean_sd_trends(design) %>%
plot_gamma()
Run the code above in your browser using DataLab