R/GSEA.write.gct.R

In GSEA-MSigDB/GSEA_R: Gene Set Enrichment Analysis (GSEA)

#' Write output .GCT files
#'
#' `GSEA.write.gct` produces properly formatted .GCT files of leading edge subsets
#'
#' Internal `GSEA.Analyze.Sets` function.
#'
#' @keywords internal
#'

GSEA.write.gct <-
function(gct, filename) {
 f <- file(filename, "w")
 cat("#1.2", "\n", file = f, append = TRUE, sep = "\t")
 cat(dim(gct)[1], "\t", dim(gct)[2], "\n", file = f, append = TRUE, sep = "\t")
 cat("Name", "\t", file = f, append = TRUE, sep = "\t")
 cat("Description", file = f, append = TRUE, sep = "\t")
 names <- names(gct)
 cat("\t", names[1], file = f, append = TRUE, sep = "\t")
 for (j in 2:length(names)) {
  cat("\t", names[j], file = f, append = TRUE, sep = "\t")
 }
 cat("\n", file = f, append = TRUE, sep = "\t")
 oldWarn <- options(warn = -1)
 m <- matrix(nrow = dim(gct)[1], ncol = dim(gct)[2] + 2)
 m[, 1] <- row.names(gct)
 m[, 2] <- row.names(gct)
 index <- 3
 for (i in 1:dim(gct)[2]) {
  m[, index] <- gct[, i]
  index <- index + 1
 }
 write.table(m, file = f, append = TRUE, quote = FALSE, sep = "\t", eol = "\n",
  col.names = FALSE, row.names = FALSE, na = "")
 close(f)
 options(warn = 0)
 return(gct)
}