Skip to content

Instantly share code, notes, and snippets.

@pcantalupo
Created January 10, 2020 12:16
Show Gist options
  • Save pcantalupo/1cd4241a766abe5e783651642794f38c to your computer and use it in GitHub Desktop.
Save pcantalupo/1cd4241a766abe5e783651642794f38c to your computer and use it in GitHub Desktop.
██████╗ ███████╗███████╗ ██████╗ █████╗
██╔════╝ ██╔════╝██╔════╝██╔════╝██╔══██╗
██║ ███╗███████╗█████╗ ██║ ███████║
██║ ██║╚════██║██╔══╝ ██║ ██╔══██║
╚██████╔╝███████║███████╗╚██████╗██║ ██║
╚═════╝ ╚══════╝╚══════╝ ╚═════╝╚═╝ ╚═╝
[*] loading datasets ...
[*] GSECA running ...
=== FMM and DD ===
[*] NE threshold: 0.01
[*] Number of Expression Classes: 7
-- Parallelizing 2 cores ...
[*] FMM Getting thresholds ...
[*] DD Setting expression classes ...
|======================================================================| 100%
|======================================================================| 100%
=== Gene Set Enrichment Class Analysis ===
[*] GSECA core ...
[*] GSECA scoring ...
|======================================================================| 100%
[*] Empirical P-Value estimation ...
R Version: R version 3.6.1 (2019-07-05)
snowfall 1.84-6.1 initialized (using snow 0.4-3): parallel execution on 2 CPUs.
Library parallel loaded.
Library parallel loaded in cluster.
Killed
The command '/bin/sh -c R -e "source('Scripts/config.R');M = read.delim('Examples/PRAD.ptenloss.M.tsv');L = read.delim('Examples/PRAD.ptenloss.L.tsv')[,1];pl = read.gmt.file('gene_sets/cereda.158.KEGG.gmt');res = GSECA_executor(M, L, pl, outdir = 'Results', analysis = 'my_analysis', N.CORES = 2, EMPIRICAL = T, BOOTSTRP = T, nsim = 2, AS = 0.25, PEMP = 1, SR = 0.7, toprank = 10, iphen = c('CASE', 'CNTR'));message(typeof(res));"' returned a non-zero code: 137
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment