From 1739e9f5706bb8a73de5dbf0b467de49ea040898 Mon Sep 17 00:00:00 2001
From: Ben Green <bgreen@g.harvard.edu>
Date: Mon, 8 Jun 2015 15:21:51 -0400
Subject: added my R scripts

---
 R Scripts/-compare-sim.R | 34 ++++++++++++++++++++++++++++++++++
 1 file changed, 34 insertions(+)
 create mode 100755 R Scripts/-compare-sim.R

(limited to 'R Scripts/-compare-sim.R')

diff --git a/R Scripts/-compare-sim.R b/R Scripts/-compare-sim.R
new file mode 100755
index 0000000..d96bd81
--- /dev/null
+++ b/R Scripts/-compare-sim.R	
@@ -0,0 +1,34 @@
+
+vics = vic_ids
+graph = lcc
+
+n.infections = length(vics)
+vic.neighbors = rep(0,n.infections)
+for (i in 1:n.infections){
+  u = vics[i]
+  nbhd = unlist(neighborhood(graph, nodes=u, order=3))
+  nbhd = intersect(vic_ids,nbhd)
+  vic.neighbors[i] = sum(nbhd %in% vics)
+}
+
+
+n.infections = length(vics)
+vic.time = rep(0,n.infections)
+for (i in 1:n.infections){
+  if (i%%1000==0) print(i)
+  u = vics[i]
+  nbhd = unlist(neighborhood(graph, nodes=u, order=1))
+  nbhd = intersect(vic_ids,nbhd)
+  nbhd = setdiff(nbhd,u)
+  tu = as.numeric(V(graph)$vic_date[u])
+  tvs = as.numeric(V(graph)$vic_date[nbhd])
+  tvs = tvs[tvs<tu]
+  if (length(tvs)>0){
+    vic.time[i] = min(tu-tvs)
+  }
+}
+vic.time = vic.time[vic.time>0]
+mean(vic.time)
+median(vic.time)
+mean(vic.time<100)
+sum(vic.time>0)
-- 
cgit v1.2.3-70-g09d2