corrected age to account for each node’s spawn date

1 files changed, 15 insertions, 0 deletions

diff --git a/R Scripts/min_inf_time.R b/R Scripts/min_inf_time.R
new file mode 100644
index 0000000..914af61
--- /dev/null
+++ b/R Scripts/min_inf_time.R
@@ -0,0 +1,15 @@
+# analyze min possible infection time
+i = 1
+min_time = 0#rep(Inf,length(unique(dag_dat_vics$to)))
+min_time_dist = 0#rep(Inf,length(unique(dag_dat_vics$to)))
+for(to in unique(dag_dat_vics$to)){
+  rows = which(dag_dat_vics$to==to & dag_dat_vics$dist<2)
+  if(length(rows)>0){
+    min_time[i] = min(dag_dat_vics$t2[rows]-dag_dat_vics$t1[rows])
+    min_time_dist[i] = dag_dat_vics$dist[rows[which.min(dag_dat_vics$t2[rows]-dag_dat_vics$t1[rows])]]
+    i = i + 1
+  }
+}
+median(min_time)
+mean(min_time<100)
+save(min_time_1,min_time_2,min_time_3,file='Results/min_inf_time.RData')