Method creation for ssGSEA

Jigyasa-G · Aug 8, 2023 · 6fd974a · 6fd974a
1 parent c2ea40a
commit 6fd974a
Showing 1 changed file with 14 additions and 0 deletions.
diff --git a/R/ssGSEA.R b/R/ssGSEA.R
@@ -1 +1,15 @@
 #Perform ssGSEA Analysis
+
+ssgsea <- function(gene_expression, genes, num_gene_sets, genes_per_set) {
+  gs <- list()
+  for (i in 1:num_gene_sets) {
+    sampled_genes <- sample(genes, genes_per_set, replace = FALSE)
+    gs[[i]] <- sampled_genes
+  }
+  names(gs) <- paste0("gs", 1:length(gs))
+
+  ssGSEA_scores <- gsva(gene_expression, gs, method = "ssgsea")
+  return(ssGSEA_scores)
+}
+
+# ssgsea(gene_expression, genes, num_gene_sets, genes_per_set)