@@ -11,13 +11,15 @@ class CreateOrganAgeObject:
1111
1212 # init method or constructor
1313 def __init__ (self ,
14- path_version_scale_factors = 'v4_to_v4.1_scale_dict.json' ,
14+ path_version_scale_factors_v4 = 'v4_to_v4.1_scale_dict.json' ,
15+ path_version_scale_factors_v5 = 'v5_to_v4.1_scale_dict.json' ,
1516 path_organ_plist_dict1 = 'tissue_pproteinlist_5k_dict_gtex_tissue_enriched_fc4_stable_assay_proteins_seqid.json' ,
1617 path_organ_plist_dict2 = 'tissue_pproteinlist_5k_dict_dementia_optimized_WUADRC_WUADRC_trained_stableassayps_seqid.json' ,
1718 path_bootstrap_seeds = 'Bootstrap_and_permutation_500_seed_dict.json' ,
1819 ):
1920
20- self .data_and_model_paths = {"path_version_scale_factors" : path_version_scale_factors ,
21+ self .data_and_model_paths = {"path_version_scale_factors_v4" : path_version_scale_factors_v4 ,
22+ "path_version_scale_factors_v5" : path_version_scale_factors_v5 ,
2123 "path_organ_plist_dict1" : path_organ_plist_dict1 ,
2224 "path_organ_plist_dict2" : path_organ_plist_dict2 ,
2325 "path_bootstrap_seeds" : path_bootstrap_seeds ,
@@ -29,8 +31,10 @@ def __init__(self,
2931 def load_data_and_models (self ):
3032
3133 # Seqid:scale_factor dictionary
32- version_scale_factors = json .load (resources .open_text ("organage.data" , self .data_and_model_paths ["path_version_scale_factors" ]))
33- self .version_scale_factors = version_scale_factors
34+ version_scale_factors_v4 = json .load (resources .open_text ("organage.data" , self .data_and_model_paths ["path_version_scale_factors_v4" ]))
35+ self .version_scale_factors_v4 = version_scale_factors_v4
36+ version_scale_factors_v5 = json .load (resources .open_text ("organage.data" , self .data_and_model_paths ["path_version_scale_factors_v5" ]))
37+ self .version_scale_factors_v5 = version_scale_factors_v5
3438
3539 # organ:proteinlist dictionary
3640 organ_plist_dict1 = json .load (resources .open_text ("organage.data" ,
@@ -112,7 +116,10 @@ def normalize(self, assay_version):
112116 df_prot_norm = self .df_prot .copy ()
113117 if assay_version == "v4" :
114118 for prot in df_prot_norm .columns :
115- df_prot_norm [prot ] = df_prot_norm [prot ] * self .version_scale_factors [prot ]
119+ df_prot_norm [prot ] = df_prot_norm [prot ] * self .version_scale_factors_v4 [prot ]
120+ if assay_version == "v5" :
121+ for prot in df_prot_norm .columns :
122+ df_prot_norm [prot ] = df_prot_norm [prot ] * self .version_scale_factors_v5 [prot ]
116123 if assay_version == "v4.1" :
117124 pass
118125
0 commit comments