# HG changeset patch # User proteore # Date 1542806533 18000 # Node ID c8bff18c102b372cf4fd2f73cbdefa847080c7de # Parent 3a66e1e65c0276686b53bb21345c82528374c393 planemo upload commit d884b75cfc292105807c42c58f864d536dc9e431-dirty diff -r 3a66e1e65c02 -r c8bff18c102b GO-enrich.R --- a/GO-enrich.R Tue Nov 06 08:56:05 2018 -0500 +++ b/GO-enrich.R Wed Nov 21 08:22:13 2018 -0500 @@ -264,12 +264,14 @@ for (onto in ontology) { if (args$go_represent == "true") { ggo<-repartition.GO(gene, orgdb, onto, level, readable=TRUE) + ggo <- as.data.frame(apply(ggo, c(1,2), function(x) gsub("^$|^ $", NA, x))) #convert "" and " " to NA output_path = paste("cluster_profiler_GGO_",onto,".tsv",sep="") write.table(ggo, output_path, sep="\t", row.names = FALSE, quote = FALSE ) } if (args$go_enrich == "true") { ego<-enrich.GO(gene, universe_gene, orgdb, onto, pval_cutoff, qval_cutoff,plot) + ego <- as.data.frame(apply(ego, c(1,2), function(x) gsub("^$|^ $", NA, x))) #convert "" and " " to NA output_path = paste("cluster_profiler_EGO_",onto,".tsv",sep="") write.table(ego, output_path, sep="\t", row.names = FALSE, quote = FALSE ) } diff -r 3a66e1e65c02 -r c8bff18c102b cluster_profiler.xml --- a/cluster_profiler.xml Tue Nov 06 08:56:05 2018 -0500 +++ b/cluster_profiler.xml Wed Nov 21 08:22:13 2018 -0500 @@ -1,4 +1,4 @@ - + (Human, Mouse, Rat) (clusterProfiler) R