From 110069d77815a3d62e3526f18b2a34fb79beff1e Mon Sep 17 00:00:00 2001
From: Ben Green <ben@SEASITs-MacBook-Pro.local>
Date: Thu, 2 Jul 2015 23:10:54 -0400
Subject: just some minor tests with ml

---
 R Scripts/generate-network.R | 16 +++++++++-------
 1 file changed, 9 insertions(+), 7 deletions(-)

(limited to 'R Scripts/generate-network.R')

diff --git a/R Scripts/generate-network.R b/R Scripts/generate-network.R
index dab81a4..db7012d 100644
--- a/R Scripts/generate-network.R	
+++ b/R Scripts/generate-network.R	
@@ -3,9 +3,9 @@ setwd("~/Documents/Cascade Project/")
 source('criminal_cascades/R Scripts/temporal.R')
 source('criminal_cascades/R Scripts/structural.R')
 
-alpha = 1/10
-beta = 0.01
-delta = 0.25
+alpha = 0.1
+beta = 0.25
+delta = 0.1
 # lmbda = 1/10
 t_max = 1000
 
@@ -16,7 +16,7 @@ plot(g, vertex.size=3, vertex.label=NA)
 V(g)$seed = runif(vcount(g))<beta
 seeds = which(V(g)$seed)
 V(g)$vic = V(g)$seed
-V(g)$vic.day[V(g)$seed] = 1#sample(1:t_max, sum(V(g)$seed))
+V(g)$vic.day[V(g)$seed] = sample(1:t_max, sum(V(g)$seed), replace=T) #1 if testing ml2
 V(g)$spawn.date = 0
 V(g)$infector = NA
 
@@ -31,9 +31,11 @@ for (day in 1:t_max){
     inf.days = day + ceiling(rexp(length(infected),alpha))
     realized = ((inf.days <= V(g)$vic.day[infected]) %in% c(NA,T)) & (inf.days<=t_max)
     infected = infected[realized]
-    V(g)$vic[infected] = TRUE
-    V(g)$vic.day[infected] = inf.days[realized]
-    V(g)$infector[infected] = vic
+    if(sum(realized)){
+      V(g)$vic[infected] = TRUE
+      V(g)$vic.day[infected] = inf.days[realized]
+      V(g)$infector[infected] = vic
+    }
   }
 }
 
-- 
cgit v1.2.3-70-g09d2