-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathMUTPARAMS.sh~
executable file
·83 lines (56 loc) · 2.72 KB
/
MUTPARAMS.sh~
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
# Centralized parameter file
# ************************************************************************************
# INPUT OUTPUT
# ************************************************************************************
export MUTATIONMATRIXPATH=/Users/jlu96/conte/jlu/mutex/data/ACC_broad-som.m2
# This will be prefix all of the output files.
export OUTPUTPREFIX=/Users/jlu96/conte/jlu/mutex/output/ACC_broad-som
# ************************************************************************************
# LOADING PARAMETERS
# ************************************************************************************
# Limit to genes mutated in at least this # of patients
export MINFREQ=
# Limit to genes mutated in at least this % of patients
export MINPERCENTILE=2
# Limit to this # of top most frequently mutated genes
export TOP_NUMBER=
# List of patients to limit loading to
export PATIENTFILE=$FOLDER/shared_patients.plst
# List of Genes to limit loading to
export GENEFILE=COSMICGenes_OnlyLoss.txt
# Blacklists
export PATIENTBLACKLIST=
export GENEBLACKLIST=
# ************************************************************************************
# LIMIT WHICH PAIRS TO TEST
# ************************************************************************************
# If GENELIST1 and GENELIST2 are both set, only pairs containing one gene in GENELIST1 and one in
# GENELIST2 will be tested. If only GENELIST1 is set, only pairs containing two genes in GENELIST1
# Will be tested.
export GENELIST1=
export GENELIST2=
# If set, only test pairs from PAIRLIST. The column header must be "Gene0 Gene1"
export PAIRLIST=
# If set, only test genes that are at least this distance away.
# This will only apply to genes ending in "gain" or "loss".
# E.G. Say TP53 and MYC are 30 bp away. If the threshold is set to 1000000, the pair (TP53gain, MYCloss) will not
# be tested. But the pair (TP53, MYCloss) would be tested. This allows only Copy Number pairs to be
# filtered out
export MINDISTANCETHRESHOLD=100000000
# ************************************************************************************
# CALCULATION PARAMETERS
# ************************************************************************************
# Mutex Pair PValue Threshold. Default is 0.05
export MPROB=
# Maximum Overlap between mutually exclusive pairs.
export MAXOVERLAP=
# Cooccur Pair PValue Threshold. Default is 0.05
export CPROB=
# Minimum overlap between cooccurring pairs
export MINCOOCCUR=
# minimum co-occurrence ratio for cooccuring pairs. This is the overlap/
export MINRATIO=
# If "Triplet", search for triplets among pairs. If "Network", calculate Network attributes
export GROUPTYPE=Triplet
# Number of processes to start and divide parameter testing across.
export NUMPROCESSES=3