-
Notifications
You must be signed in to change notification settings - Fork 5
/
Copy pathQTL_config_template.R
34 lines (34 loc) · 1.09 KB
/
QTL_config_template.R
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
#ID for user's trait of interest. (Can be any string)
trait = "TRAITNAME"
#path to the input files
traitFilePath = "TRAITPATH"
#column IDs from trait file
trait_A1col = "A1COL"
trait_A2col = "A2COL"
trait_SNPcol = "SNPCOL"
trait_CHRcol = "CHRCOL"
trait_BPcol = "BPCOL"
trait_Pcol = "PCOL"
trait_Ncol = "NCOL"
trait_MAFcol = "MAFCOL"
#trait info not in the input file
#traitType is set either to "cc" or "quant"
traitType = "TRAITTYPE"
#This is the proportion of samples that are cases in a case control GWAS, if you are using a quantitative trait this should be set to "". traitProp = cases / case + controls
traitProp = TRAITPROP #look this up
#locus information for running coloc.
chrom = CHROMOSOME
colocStart = STARTBP
colocStop = STOPBP
#reference genome build: "hg19" or "hg38"
build = "BUILD"
lead_SNP = "SNPNUMBER"
#"eqtl" or "sqtl"
qtlType = "QTLTYPE"
#plink parameters
clump_P1 = CLUMPP1
clump_KB = CLUMPKB
clump_R2 = CLUMPR2
#config file with paths to the qtl data and plink files
#setup_config_R = "/project/voight_GWAS/bychen9/eQTL_colocalizer/setup_config.R"
setup_config_R = "SETUPCONFIGR"