tanu/LSHTM_analysis
1
0
Fork 0
Code Issues Pull requests Releases Wiki Activity

refactoring or_kin script minor changes only

Browse source
This commit is contained in:
Tanushree Tunstall 2020-07-10 12:37:41 +01:00
parent aaf3f5e084
commit c80faef0bf
4 changed files with 10 additions and 15 deletions
Download patch file Download diff file Expand all files Collapse all files

6
scripts/data_extraction.py
View file

@ -69,7 +69,9 @@ from tidy_split import tidy_split
#%% command line args
arg_parser = argparse.ArgumentParser()
arg_parser.add_argument('-d', '--drug', help='drug name', default = None)
arg_parser.add_argument('-g', '--gene', help='gene name (case sensitive)', default = None) # case sensitive
arg_parser.add_argument('-g', '--gene', help='gene name (case sensitive)', default = None)
args = arg_parser.parse_args()
#=======================================================================
#%% variable assignment: input and output paths & filenames
@ -82,7 +84,7 @@ gene = 'pncA'
gene_match = gene + '_p.'
print('mut pattern for gene', gene, ':', gene_match)
nssnp_match = gene_match+'[A-Z]{3}[0-9]+[A-Z]{3}'
nssnp_match = gene_match +'[A-Z]{3}[0-9]+[A-Z]{3}'
print('nsSNP for gene', gene, ':', nssnp_match)
wt_regex = gene_match.lower()+'(\w{3})'

6
scripts/dssp_df.py
View file

@ -52,7 +52,7 @@ drug = args.drug
gene = args.gene
gene_match = gene + '_p.'
data_dir = args.datadir
datadir = args.datadir
indir = args.input_dir
outdir = args.output_dir
@ -63,9 +63,7 @@ DEBUG = args.debug
#=============
# directories
#=============
if data_dir:
datadir = data_dir
else:
if not datadir:
datadir = homedir + '/' + 'git/Data'
if not indir:

6
scripts/kd_df.py
View file

@ -60,7 +60,7 @@ drug = args.drug
gene = args.gene
gene_match = gene + '_p.'
data_dir = args.datadir
datadir = args.datadir
indir = args.input_dir
outdir = args.output_dir
@ -72,9 +72,7 @@ DEBUG = args.debug
#============
# directories
#============
if data_dir:
datadir = data_dir
else:
if not datadir:
datadir = homedir + '/' + 'git/Data'
if not indir:

7
scripts/rd_df.py
View file

@ -51,7 +51,7 @@ drug = args.drug
gene = args.gene
gene_match = gene + '_p.'
data_dir = args.datadir
datadir = args.datadir
indir = args.input_dir
outdir = args.output_dir
@ -62,9 +62,7 @@ DEBUG = args.debug
#============
# directories
#============
if data_dir:
datadir = data_dir
else:
if not datadir:
datadir = homedir + '/' + 'git/Data'
if not indir:
@ -79,7 +77,6 @@ if not outdir:
if rd_filename:
in_filename_rd = rd_filename
else:
#in_filename_rd = '3pl1_rd.tsv'
in_filename_rd = gene.lower() + '_rd.tsv'
infile_rd = outdir + '/' + in_filename_rd