Mercurial > repos > kpbioteam > clusterprofiler_go

annotate clusterProfiler_go.R @ 7:e63c8d858a65 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2-dirty
author kpbioteam
date Mon, 16 Jul 2018 06:20:31 -0400
parents 6d08a0432d3f
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
1 require("clusterProfiler", quietly = TRUE)
1
49ae9762eb96 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2-dirty
kpbioteam
parents: 0
diff changeset
2 require("DO.db", quietly = TRUE)
0
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
3 require("org.Hs.eg.db", quietly = TRUE)
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
4 require("ggplot2", quietly = TRUE)
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
5
1
49ae9762eb96 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2-dirty
kpbioteam
parents: 0
diff changeset
6
0
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
7 args <- commandArgs(trailingOnly = TRUE)
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
8
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
9 input1 = args[1]
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
10 output1 = args[2]
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
11 output2 = args[3]
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
12
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
13 gene <- read.table(input1)
1
49ae9762eb96 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2-dirty
kpbioteam
parents: 0
diff changeset
14 ego <- enrichGO(gene = c(gene$V1),
7
e63c8d858a65 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2-dirty
kpbioteam
parents: 6
diff changeset
15 keyType = "ENTREZID",
0
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
16 OrgDb = org.Hs.eg.db,
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
17 ont = "BP",
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
18 pAdjustMethod = "BH",
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
19 pvalueCutoff = 0.01,
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
20 qvalueCutoff = 0.05,
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
21 readable = TRUE)
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
22
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
23 write.table(ego, output1)
5aba765e27b4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2
kpbioteam
parents:
diff changeset
24 dotplot(ego)
5
a12bbeb898a4 planemo upload for repository https://github.com/kpbioteam/clusterProfiler_go commit 19aa66aed9f704dc98680d2573ff7b13c81c0ca2-dirty
kpbioteam
parents: 4
diff changeset
25 ggsave(file = output2, device = "pdf", width = 60, height = 30, units = "cm")