# HG changeset patch # User johnheap # Date 1559765850 14400 # Node ID 92c5ecb025b6510fe53807e43f11f41777e8557b # Parent 7983a9305b5332555a249f4b1feaf1c1e02f5ff9 Uploaded diff -r 7983a9305b53 -r 92c5ecb025b6 Tryp_V_T.py --- a/Tryp_V_T.py Wed Jun 05 05:54:10 2019 -0400 +++ b/Tryp_V_T.py Wed Jun 05 16:17:30 2019 -0400 @@ -238,16 +238,16 @@ return int(sac[1:i]) def combineFPMK(tdict): - dir_path = os.path.dirname(os.path.realpath(__file__))+'/' - - fpkm_df = pd.read_csv(dir_path+tdict['name']+'.cuff/genes.fpkm_tracking', sep='\t') + # dir_path = os.path.dirname(os.path.realpath(__file__))+'/' + cur_path = os.getcwd()+'/' + fpkm_df = pd.read_csv(cur_path+tdict['name']+'.cuff/genes.fpkm_tracking', sep='\t') #fpkm_df = pd.read_csv('genes.fpkm_tracking',sep='\t') #print(fpkm_df.head()) fpkm_df['locus'] = fpkm_df['locus'].apply(lambda names: names[:names.find(':')]) #print(fpkm_df.head()) - reducedBlast_df = pd.read_csv(dir_path + tdict['name']+'_transcript.csv') + reducedBlast_df = pd.read_csv(cur_path + tdict['name']+'_transcript.csv') # reducedBlast_df = pd.read_csv('TrinityVT_transcript.csv') saccverSet = set(reducedBlast_df['saccver']) saccverList = list(saccverSet) @@ -359,7 +359,7 @@ plt.savefig(tdict['html_resource'] + tdict['name']+"_phylotypes.png") if tdict['pdf'] == 'PDF_Yes': plt.savefig(tdict['html_resource'] + tdict['name']+"phylotypes.pdf") - plt.show() +# plt.show() pass