@@ -235,8 +235,8 @@ teEnrichment <- function(inputGenes = NULL, rnaSeqDataset = 1,
235
235
236
236
# ### Calculate the Hypergeometric P-Value#########
237
237
tissueNames <- as.character(tissueDetails $ RName )
238
- print(paste0(" Total background genes:" ,nrow(geneMappingForCurrentDataset )))
239
- print(paste0(" Total input genes:" ,length(inputEnsemblGenes )))
238
+ # print(paste0("Total background genes:",nrow(geneMappingForCurrentDataset)))
239
+ # print(paste0("Total input genes:",length(inputEnsemblGenes)))
240
240
x <- lapply(tissueNames , FUN = function (tissue ) {
241
241
242
242
tissueGenes <- finalTissueSpecificGenes %> %
@@ -309,7 +309,7 @@ teEnrichment <- function(inputGenes = NULL, rnaSeqDataset = 1,
309
309
nTeGenesInTissue <- nrow(tissueGenes )
310
310
nTotalGenes <- nrow(geneMappingForCurrentDataset )
311
311
nTotalInputGenes <- length(inputEnsemblGenes )
312
- print(paste0(" Total TS genes in " ,tissue ," " ,nTeGenesInTissue ))
312
+ # print(paste0("Total TS genes in ",tissue," ",nTeGenesInTissue))
313
313
pValue <- stats :: phyper(overlapGenes - 1 , nTeGenesInTissue ,
314
314
nTotalGenes - nTeGenesInTissue ,
315
315
nTotalInputGenes , lower.tail = FALSE )
0 commit comments