This article produces a data table for presentation purposes (see results/data-table-for-regression.csv). Note that these data were used for previous regression analyses.

source("../../src/import-data.R")
gene.ids <- unlist(read.csv("../../data/genes.regression.new", as.is = TRUE))
names(gene.ids) <- gene.ids
Y <- data.frame(lapply(get.readcounts(gene.ids)[gene.ids], getElement, "S"))
X <- get.predictors()[1:13]
write.csv(Z <- cbind(Y, X), "../../results/data-table-for-regression.csv")