# HG changeset patch # User iuc # Date 1548667660 18000 # Node ID bbcf5f7f2af20ca60aae6f2c712d76a37bbd16f7 # Parent ae39895af5fe825367132273d5abf197441a15a0 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/goseq commit 2dffcdff171ebf695ef0fc7eb8f8264f0d132115 diff -r ae39895af5fe -r bbcf5f7f2af2 goseq.r --- a/goseq.r Sun Sep 30 09:27:05 2018 -0400 +++ b/goseq.r Mon Jan 28 04:27:40 2019 -0500 @@ -134,7 +134,7 @@ sink(zz) GO.samp=goseq(pwf, genome = genome, id = gene_id, method="Sampling", repcnt=repcnt, use_genes_without_cat = use_genes_without_cat, gene2cat=go_map) sink() - + GO.samp$p.adjust.over_represented = p.adjust(GO.samp$over_represented_pvalue, method=p_adj_method) GO.samp$p.adjust.under_represented = p.adjust(GO.samp$under_represented_pvalue, method=p_adj_method) write.table(GO.samp, sampling_tab, sep="\t", row.names = FALSE, quote = FALSE) @@ -156,15 +156,15 @@ pdf("top10.pdf") for (m in names(results)) { p <- results[[m]] %>% - top_n(10, wt=-p.adjust.over_represented) %>% + top_n(10, wt=-over_represented_pvalue) %>% mutate(hitsPerc=numDEInCat*100/numInCat) %>% ggplot(aes(x=hitsPerc, y=substr(term, 1, 40), # only use 1st 40 chars of terms otherwise squashes plot - colour=p.adjust.over_represented, + colour=over_represented_pvalue, size=numDEInCat)) + geom_point() + expand_limits(x=0) + - labs(x="% DE in category", y="Category", colour="adj. P value", size="Count", title=paste("Top over-represented categories in", cats_title), subtitle=paste(m, " method")) + + labs(x="% DE in category", y="Category", colour="P value", size="Count", title=paste("Top over-represented categories in", cats_title), subtitle=paste(m, " method")) + theme(plot.title = element_text(hjust = 0.5), plot.subtitle = element_text(hjust = 0.5)) print(p) } diff -r ae39895af5fe -r bbcf5f7f2af2 goseq.xml --- a/goseq.xml Sun Sep 30 09:27:05 2018 -0400 +++ b/goseq.xml Mon Jan 28 04:27:40 2019 -0500 @@ -1,14 +1,14 @@ - + tests for overrepresented gene categories + bioconductor-goseq + bioconductor-org.hs.eg.db + bioconductor-org.dm.eg.db + bioconductor-org.dr.eg.db + bioconductor-org.mm.eg.db + r-dplyr + r-ggplot2 r-optparse - bioconductor-goseq - bioconductor-org.hs.eg.db - bioconductor-org.dm.eg.db - bioconductor-org.dr.eg.db - bioconductor-org.mm.eg.db - r-dplyr - r-ggplot2 - +