added my R scripts

1 files changed, 34 insertions, 0 deletions

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)