comparison topGO_enrichment.R @ 12:8eaa43ba1bfc draft

planemo upload commit 4ba1ebe7b3f5e3fabf78b5fed7ed0b92e2cbf9e5-dirty

11:fa2e27165d5d

  if (length(intersect(sample,allGenes))==0){
    print("None of the input ids can be found in the org package data, enrichment analysis cannot be realized. \n The inputs ids probably have no associated GO terms.")
    return(c(NULL,NULL))
  }
  
  geneList = factor(as.integer(allGenes %in% sample)) 
  if (length(levels(geneList)) == 1 ){
    stop("All or none of the background genes are found in tested genes dataset, enrichment analysis can't be done")
  }
  names(geneList) <- allGenes
  

12:8eaa43ba1bfc

  if (length(intersect(sample,allGenes))==0){
    print("None of the input ids can be found in the org package data, enrichment analysis cannot be realized. \n The inputs ids probably have no associated GO terms.")
    return(c(NULL,NULL))
  }
  
  geneList = factor(as.integer(allGenes %in% sample)) #duplicated ids in sample count only for one
  if (length(levels(geneList)) == 1 ){
    stop("All or none of the background genes are found in tested genes dataset, enrichment analysis can't be done")
  }
  names(geneList) <- allGenes