13.9 Load libraries and datasets

library(tidyverse)
library(rafalib)

# BiocManager::install(c("genefilter", "qvalue"))
library(genefilter) # rowttests function
library(qvalue) # to calculate a qvalue and compare with p.adjust(...method="fdr")

geneExpression <- readRDS("data/13_geneExpression.rds")
sampleInfo<- readRDS("data/13_sampleInfo.rds")