-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathrunAgainstTCGA.py
25 lines (23 loc) · 1.33 KB
/
runAgainstTCGA.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
from subprocess import *
import os
# Run for each tumor type in TCGA
for tumor in ["ACC", 'BLCA', 'CESC', 'CHOL', 'COAD', 'DLBC', "PCPG","READ","SARC","SKCM","TGCT","THCA","THYM","UCEC","UCS","UVM","ESCA","GBM","HNSC","KICH","KIRC","KIRP","LGG","LIHC","LUAD","LUSC","MESO","OV","PAAD", "PRAD", "STAD"]:
print(tumor)
# Make the output directory if it doesn't exists already
# Run with both filters
print('\nRunning pq_mff...')
cmd4 = ['python oncoMergePathways.py',
'-gp C:/Users/nihan/Dropbox/OncoMerge_pathways/TCGA/GISTIC/'+tumor,
'-aaf C:/Users/nihan/Dropbox/OncoMerge_pathways/TCGA/OncoMerge_input_g2e_converter.csv',
'-ln "Locus ID"',
'-pam C:/Users/nihan/Dropbox/OncoMerge_pathways/TCGA/PAM/'+tumor+'_somMutMC3.csv',
'-mscv C:/Users/nihan/Dropbox/OncoMerge_pathways/TCGA/MutSig2cv/'+tumor+'_sig2cv.csv',
'-fus C:/Users/nihan/Dropbox/OncoMerge_pathways/TCGA/FUSIONS/'+tumor+'_fusions.csv',
'-op C:/Users/nihan/Documents/Oncomerge/output/'+tumor+'/pq_mff',
'-pq 0.1',
'-mlg 10',
'-mpf 0',
'-tcga True',
'-bl C:/Users/nihan/Dropbox/OncoMerge_pathways/TCGA/blacklist/blacklist_29850653_29625053.csv']
proc4 = Popen(' '.join(cmd4), shell=True, stdin=PIPE)
out = proc4.communicate()