Mercurial > repos > bimib > cobraxy

comparison COBRAxy/flux_simulation_beta.py @ 480:ecaf3c9852b3 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author luca_milaz
date Mon, 22 Sep 2025 16:29:01 +0000
parents a036b8d3642f
children 1c6fcf97cabe
comparison
equal deleted inserted replaced
479:a036b8d3642f 480:ecaf3c9852b3
210 # Merge all batches into a single DataFrame 210 # Merge all batches into a single DataFrame
211 all_samples = [] 211 all_samples = []
212 212
213 for i in range(n_batches): 213 for i in range(n_batches):
214 batch_filename = f"{ARGS.output_path}/{model_name}_{i}_OPTGP.npy" 214 batch_filename = f"{ARGS.output_path}/{model_name}_{i}_OPTGP.npy"
215 batch_data = np.load(batch_filename) 215 batch_data = np.load(batch_filename, allow_pickle=True)
216 all_samples.append(batch_data) 216 all_samples.append(batch_data)
217 217
218 # Concatenate all batches 218 # Concatenate all batches
219 samplesTotal_array = np.vstack(all_samples) 219 samplesTotal_array = np.vstack(all_samples)
220 220
287 # Merge all batches into a single DataFrame 287 # Merge all batches into a single DataFrame
288 all_samples = [] 288 all_samples = []
289 289
290 for i in range(n_batches): 290 for i in range(n_batches):
291 batch_filename = f"{ARGS.output_path}/{model_name}_{i}_CBS.npy" 291 batch_filename = f"{ARGS.output_path}/{model_name}_{i}_CBS.npy"
292 batch_data = np.load(batch_filename) 292 batch_data = np.load(batch_filename, allow_pickle=True)
293 all_samples.append(batch_data) 293 all_samples.append(batch_data)
294 294
295 # Concatenate all batches 295 # Concatenate all batches
296 samplesTotal_array = np.vstack(all_samples) 296 samplesTotal_array = np.vstack(all_samples)
297 297
298 # Convert back to DataFrame with proper column names 298 # Convert back to DataFrame with proper column namesq
299 samplesTotal = pd.DataFrame(samplesTotal_array, columns=reaction_ids) 299 samplesTotal = pd.DataFrame(samplesTotal_array, columns=reaction_ids)
300 300
301 # Save the final merged result as CSV 301 # Save the final merged result as CSV
302 write_to_file(samplesTotal.T, model_name, True) 302 write_to_file(samplesTotal.T, model_name, True)
303 303