# NOT RUN {
data(design)
beta0 <- NULL
regression.output <- res$ori.res$newlm[1:50,]
ncores <- 1
C.matrix <- list()
C.matrix[[1]] <- limma::makeContrasts(line2, levels = design)
C.matrix[[2]] <- limma::makeContrasts(time2, time6, time24, levels = design)
names(C.matrix) <- c("line2","time")
teststatout <- rmRNAseq:::teststat(C.matrix, beta0, regression.output, ncores)
head(teststatout)
# }
Run the code above in your browser using DataCamp Workspace