Mercurial > repos > iuc > fgsea
comparison fgsea.R @ 2:17eb1e0d711f draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/fgsea commit 7f92023de5d3d922ad6cfd22f4de263b56127a1c
author | iuc |
---|---|
date | Mon, 06 May 2019 04:45:25 -0400 |
parents | 9bb7943b5263 |
children |
comparison
equal
deleted
inserted
replaced
1:101b208a3e1a | 2:17eb1e0d711f |
---|---|
54 } else { | 54 } else { |
55 pathways <- load(sets_file) | 55 pathways <- load(sets_file) |
56 pathways <- get(pathways) | 56 pathways <- get(pathways) |
57 } | 57 } |
58 | 58 |
59 # need to set seed for reproducibility https://github.com/ctlab/fgsea/issues/12 | |
60 set.seed(42) | |
59 fgseaRes <- fgsea(pathways, ranks, minSize=min_size, maxSize=max_size, nperm=n_perm) | 61 fgseaRes <- fgsea(pathways, ranks, minSize=min_size, maxSize=max_size, nperm=n_perm) |
60 fgseaRes <- fgseaRes[order(pval), ] | 62 fgseaRes <- fgseaRes[order(pval), ] |
61 # Convert leadingEdge column from list to character to output | 63 # Convert leadingEdge column from list to character to output |
62 fgseaRes$leadingEdge <- sapply(fgseaRes$leadingEdge, toString) | 64 fgseaRes$leadingEdge <- sapply(fgseaRes$leadingEdge, toString) |
63 | 65 |