From e50f9e734260463c990a1fbdc088a22e5a49d00d Mon Sep 17 00:00:00 2001
From: Ben Green <ben@SEASITs-MacBook-Pro.local>
Date: Tue, 9 Jun 2015 17:42:03 -0400
Subject: corrected age to account for each node’s spawn date
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 R Scripts/min_inf_time.R | 15 +++++++++++++++
 1 file changed, 15 insertions(+)
 create mode 100644 R Scripts/min_inf_time.R

(limited to 'R Scripts/min_inf_time.R')

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')
-- 
cgit v1.2.3-70-g09d2