(config, indir, out_dir, label)
| 157 | return str("{0:0.2f}".format(stats.pearsonr(Y, pred[i].squeeze())[0])) |
| 158 | |
| 159 | def main(config, indir, out_dir, label): |
| 160 | data = pd.read_table(indir) |
| 161 | params = LoadConfig(config) |
| 162 | |
| 163 | X_train, Y_train = prepare_input(data[data['set'] == "Train"], params) |
| 164 | X_valid, Y_valid = prepare_input(data[data['set'] == "Val"], params) |
| 165 | X_test, Y_test = prepare_input(data[data['set'] == "Test"], params) |
| 166 | |
| 167 | DeepSTARR(params)[0].summary() |
| 168 | DeepSTARR(params)[1] |
| 169 | main_model, main_params = DeepSTARR(params) |
| 170 | main_model, my_history = train(main_model, X_train, Y_train, X_valid, Y_valid, main_params) |
| 171 | |
| 172 | endTime=time.time() |
| 173 | seconds=endTime-startTime |
| 174 | print("Total training time:",round(seconds/60,2),"minutes") |
| 175 | |
| 176 | dev_results = summary_statistics(X_test, Y_test[0], "test", "Dev", main_model, main_params, out_dir) |
| 177 | hk_results = summary_statistics(X_test, Y_test[1], "test", "Hk", main_model, main_params, out_dir) |
| 178 | |
| 179 | result = { |
| 180 | "AutoEAP": { |
| 181 | "means": { |
| 182 | "PCC(Dev)": dev_results, |
| 183 | "PCC(Hk)": hk_results |
| 184 | } |
| 185 | } |
| 186 | } |
| 187 | |
| 188 | with open(f"{out_dir}/final_info.json", "w") as file: |
| 189 | json.dump(result, file, indent=4) |
| 190 | |
| 191 | main_model.save(out_dir + '/' + label + '.h5') |
| 192 | |
| 193 | if __name__ == "__main__": |
| 194 | try: |
no test coverage detected