geney 1.2.38__py2.py3-none-any.whl → 1.2.39__py2.py3-none-any.whl

geney/oncosplice.py

@@ -382,7 +382,7 @@ def oncosplice(mut_id, splicing_threshold=0.5, protein_coding=True, primary_tran
 
         assert len(ref_protein) == len(cons_vector), f"Protein ({len(ref_protein)}) and conservation vector ({len(cons_vector)}) must be same length. {ref_protein}, \n>{cons_vector}\n>{transcript.cons_seq}"
 
-        missplicing = Missplicing(find_transcript_missplicing(transcript, mutations, engine=engine), threshold=splicing_threshold)
+        missplicing = Missplicing(find_transcript_missplicing(transcript, mutations, engine=engine, threshold=splicing_threshold), threshold=splicing_threshold)
         for mutation in mutations:
             transcript.pre_mrna += mutation
 

geney/pangolin_utils.py

@@ -49,6 +49,7 @@ def pang_one_hot_encode(seq):
 def pangolin_predict_probs(true_seq, models):
     # print(f"Running pangolin on: {true_seq}")
     model_nums = [0, 2, 4, 6]
+    model_nums = [0, 1, 2, 3, 4, 5, 6]
     INDEX_MAP = {0: 1, 1: 2, 2: 4, 3: 5, 4: 7, 5: 8, 6: 10, 7: 11}
 
     # seq = 'N'*5000 + true_seq + 'N'*5000
@@ -73,10 +74,7 @@ def pangolin_predict_probs(true_seq, models):
 
         scores.append(np.mean(score, axis=0))
 
-    # splicing_pred = np.array(scores).max(axis=0)
-    splicing_pred = np.array(scores).mean(axis=0)
+    splicing_pred = np.array(scores).max(axis=0)
     donor_probs = [splicing_pred[i] * donor_dinucleotide[i] for i in range(len(true_seq))]
     acceptor_probs = [splicing_pred[i] * acceptor_dinucleotide[i] for i in range(len(true_seq))]
-    # print(acceptor_probs)
-    # return donor_probs[context:-context], acceptor_probs[context:-context]
     return donor_probs, acceptor_probs

{geney-1.2.38.dist-info → geney-1.2.39.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: geney
-Version: 1.2.38
+Version: 1.2.39
 Summary: A Python package for gene expression modeling.
 Home-page: https://github.com/nicolaslynn/geney
 Author: Nicolas Lynn

{geney-1.2.38.dist-info → geney-1.2.39.dist-info}/RECORD

@@ -6,8 +6,8 @@ geney/graphic_utils.py,sha256=tjm6IDQ1BdfSeuPYzjlqAUHFQoDYH9jXTzJjKFS4Hh4,11078
 geney/gtex_utils.py,sha256=asL2lHyU5KsbWpV096vkf1Ka7hSo_RRfZqw7p5nERmE,1919
 geney/immune_utils.py,sha256=ZRni5ttrhpYBnmNr0d0ZatIbNPYs4nmQuoUO00SpsS4,5271
 geney/mutation_utils.py,sha256=C_kv2MB_L8LlhX3W2ooXjJ3uDoJ8zX1WeDtZKoBZJkI,1547
-geney/oncosplice.py,sha256=Ph8AMR2L9wjh3Gl3X58Dd8F1ea4on9jXIxg7EO0zVEU,22565
-geney/pangolin_utils.py,sha256=ETTGpuaQgdZ1v8H0NP8sbTEfGWu0VXUFUS7wsURsTc4,2991
+geney/oncosplice.py,sha256=J_nFs_xBSJtgMqeHv628QodRL0B2d-Zi1Ke7Pk7S4R4,22595
+geney/pangolin_utils.py,sha256=lLmnjJdJjqwWS85-1jlPLIjD2z14sWjzU87hS-8xxpQ,2873
 geney/power_utils.py,sha256=MehZFUdkJ2EFUot709yPEDxSkXmH5XevMebX2HD768A,7330
 geney/seqmat_utils.py,sha256=YV5DFLbfjXLIswPGvqK1-eEfwn9TUby0b2kewdGAKws,18372
 geney/spliceai_utils.py,sha256=gIGPC8u3J15A7EQrk2Elho5PbF9MmUUNopGGH-eEV8s,1873
@@ -20,7 +20,7 @@ geney/translation_initiation/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NM
 geney/translation_initiation/tis_utils.py,sha256=AF3siFjuQH-Rs44EV-80zHdbxRMvN4woLFSHroWIETc,4448
 geney/translation_initiation/resources/kozak_pssm.json,sha256=pcd0Olziutq-6H3mFWDCD9cujQ_AlZO-iiOvBl82hqE,1165
 geney/translation_initiation/resources/tis_regressor_model.joblib,sha256=IXb4DUDhJ5rBDKcqMk9zE3ECTZZcdj7Jixz3KpoZ7OA,2592025
-geney-1.2.38.dist-info/METADATA,sha256=a5g8X4PPpand3B4nvkbemKP4Sj066AVa7AqIjLBEGSU,948
-geney-1.2.38.dist-info/WHEEL,sha256=fS9sRbCBHs7VFcwJLnLXN1MZRR0_TVTxvXKzOnaSFs8,110
-geney-1.2.38.dist-info/top_level.txt,sha256=O-FuNUMb5fn9dhZ-dYCgF0aZtfi1EslMstnzhc5IIVo,6
-geney-1.2.38.dist-info/RECORD,,
+geney-1.2.39.dist-info/METADATA,sha256=i-1cERTVvL1H3UBw9RJbGNT6KaySdTAAHhXxL9mEo4M,948
+geney-1.2.39.dist-info/WHEEL,sha256=fS9sRbCBHs7VFcwJLnLXN1MZRR0_TVTxvXKzOnaSFs8,110
+geney-1.2.39.dist-info/top_level.txt,sha256=O-FuNUMb5fn9dhZ-dYCgF0aZtfi1EslMstnzhc5IIVo,6
+geney-1.2.39.dist-info/RECORD,,