From 5985e217801370b967e572a8627ac2b8f1509afc Mon Sep 17 00:00:00 2001 From: Jannis Teunissen Date: Fri, 13 Oct 2023 20:17:21 +0200 Subject: [PATCH] Allow variable log size in sensitivity analysis --- tools/sensitivity_analyze_results.py | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/tools/sensitivity_analyze_results.py b/tools/sensitivity_analyze_results.py index 08f2d9fc..1607c21d 100755 --- a/tools/sensitivity_analyze_results.py +++ b/tools/sensitivity_analyze_results.py @@ -22,6 +22,13 @@ logs = sorted(args.logs) logs_df = [pd.read_csv(f, delim_whitespace=True) for f in args.logs] +log_sizes = np.array([len(df) for df in logs_df]) +max_size, min_size = log_sizes.max(), log_sizes.min() + +if max_size > min_size: + print(f'Warning: logs have different size, truncating to {min_size} rows') + logs_df = [df.head(min_size) for df in logs_df] + all_cases = {} for log, df in zip(logs, logs_df): @@ -44,7 +51,7 @@ effect_magnitudes = np.zeros(len(reaction_ix)) -print(f'Using data at time t = {times[args.time_index]}') +print(f'Using data at time t = {times[args.time_index]}\n') # Here mu indicates the mean derivative, mustar the mean absolute derivative, # and sigma the standard deviation in the derivative. All derivatives w.r.t. a