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

saving work for the day

Browse source
This commit is contained in:
Tanushree Tunstall 2020-04-11 19:00:39 +01:00
parent 7d2241ad81
commit 45889990e7
2 changed files with 1 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

3
mcsm/format_results.py
View file

@ -86,7 +86,7 @@ mcsm_data.rename(columns = my_colnames_dict, inplace = True)
mcsm_data.columns mcsm_data.columns
#%%=========================================================================== #%%===========================================================================
# populate mutationinformation column:mcsm style muts {WT}<POS>{MUT} # populate mutationinformation column:mcsm style muts {WT}<POS>{MUT}
print('populating column : Mutationinformation which is currently empty\n', mcsm_data['Mutationinformation']) print('Populating column : Mutationinformation which is currently empty\n', mcsm_data['Mutationinformation'])
mcsm_data['Mutationinformation'] = mcsm_data['Wild_type'] + mcsm_data['Position'].astype(str) + mcsm_data['Mutant_type'] mcsm_data['Mutationinformation'] = mcsm_data['Wild_type'] + mcsm_data['Position'].astype(str) + mcsm_data['Mutant_type']
print('checking after populating:\n', mcsm_data['Mutationinformation'] print('checking after populating:\n', mcsm_data['Mutationinformation']
, '\n===================================================================') , '\n===================================================================')
@ -114,7 +114,6 @@ else:
#%%=========================================================================== #%%===========================================================================
# create DUET_outcome column: classification based on DUET stability values # create DUET_outcome column: classification based on DUET stability values
print('Assigning col: DUET_outcome based on DUET stability values') print('Assigning col: DUET_outcome based on DUET stability values')
print('Sanity check:') print('Sanity check:')
# count positive values in the DUET column # count positive values in the DUET column
c = mcsm_data[mcsm_data['DUET_change_kcalpermol']>=0].count() c = mcsm_data[mcsm_data['DUET_change_kcalpermol']>=0].count()

3
mcsm/run_mcsm.py
View file

@ -168,9 +168,6 @@ def write_result_url(holding_page, out_result_url):
myfile.close() myfile.close()
print(myfile) print(myfile)
# return url # return url
#======================================================================= #=======================================================================
#%% call functions #%% call functions
mcsm_muts = format_data(infile_snps) mcsm_muts = format_data(infile_snps)