diff --git a/run_analysis.R b/run_analysis.R
index affe7bd89a5d8ed37897f21b5418ad1ab389e857..c83c6fa1b18e2b90db789fd1865ccec0bc1e4d1d 100644
--- a/run_analysis.R
+++ b/run_analysis.R
@@ -111,3 +111,6 @@ names(run0)
 N <- run0$N ; dim(N)
 plot_traj(N, xlab = "Annee", ylab = "Taille de population (totale)")
 abline(h = K)
+
+out = run0
+get_metrics(N = out$N)$scenario$impact[time_horzion, "avg",-1]