Update prioritize_variants.py

Bugfix

aidiva/variant_prioritization/prioritize_variants.py

@@ -225,11 +225,11 @@ def compute_hpo_relatedness_and_final_score(variant):
 
                 if gene_distances or gene_distances_interacting:
                     # take only the maximum HPO relatedness to prevent downvoting of genes if no HPO relation in interacting genes is observed
-                    hpo_relatedness = max(max(gene_distances, default=0.0), max(gene_distances_interacting, default=0.0))
-                    #hpo_relatedness_interacting = max(gene_distances_interacting, default=0.0)
+                    hpo_relatedness = max(gene_distances, default=0.0)
+                    hpo_relatedness_interacting = max(gene_distances_interacting, default=0.0)
                     ## TODO: try different weighting of AIDIVA_SCORE and HPO_RELATEDNESS and HPO_RELATEDNESS_INTERACTING (eg 0.6 and 0.3 and 0.1)
                     # predicted pathogenicity has a higher weight than the HPO relatedness
-                    final_score = (float(variant["AIDIVA_SCORE"]) * 0.67 + float(hpo_relatedness) * 0.33) #+ float(hpo_relatedness_interacting) * 0.1)# / 3
+                    final_score = (float(variant["AIDIVA_SCORE"]) * 0.67 + float(max(hpo_relatedness,hpo_relatedness_interacting)) * 0.33) #+ float(hpo_relatedness_interacting) * 0.1)# / 3
                     #final_score = (float(variant["AIDIVA_SCORE"]) + float(hpo_relatedness) + float(hpo_relatedness_interacting)) / 3
             else:
                 final_score = np.nan