model replicates the graph in the paper

Author: cruizm

epidemic_hoya_et_al/SIR_Hoya.ma

@@ -4,30 +4,40 @@ components : area
 
 [area]
 type : cell
-width : 21
-height : 21
+width : 50
+height : 50
 delay : transport
-defaultDelayTime : 100
+defaultDelayTime : 1
 border : wrapped
 neighbors : 			area(-1,0) 
 neighbors : area(0,-1)  area(0,0)  area(0,1)
 neighbors : 			area(1,0) 
-initialvalue : 0
+initialvalue : -1
 
 
 localtransition : sir-hoya-rule
 
-statevariables: population v c m
-statevalues: 10 0.6 1 1
+statevariables: population v c m i_infec i_rec i_sus
+statevalues: 100 0.6 1.0 0.5 0.0 0.0 1.0
+initialvariablesvalue: SIR_Hoya.var
 
-
-neighborports : value infec rec pop sus
+neighborports : initial infec rec pop sus
 
 [sir-hoya-rule]
-rule : { ~pop := $population;
-~infec:= (1 - #macro(e)) * (0,0)~infec + $v*(0,0)~sus*(0,0)~infec + (0,0)~sus*#macro(effect_neighbours);
-~sus := (0,0)~sus - $v*(0,0)~sus*(0,0)~infec - (0,0)~sus*#macro(effect_neighbours);
-~rec := (0,0)~rec + #macro(e)*(0,0)~infec;
-} 100 { t }
+rule : {~initial := 0; ~pop := $population; ~infec := $i_infec; ~sus := $i_sus; ~rec := $i_rec; }
+ 1 {(0,0)~initial = -1}
+
+rule : {~pop := $population; ~infec := $i_infec; ~sus := $i_sus; ~rec := $i_rec; }
+{$i_infec:= trunc(min((1 - #macro(e)) * (0,0)~infec + $v*(0,0)~sus*(0,0)~infec + (0,0)~sus*#macro(effect_neighbours),1)*100)/100;
+$i_rec := trunc(min((0,0)~rec + #macro(e)*(0,0)~infec,1)*100)/100;
+$i_sus := 1 - $i_rec - $i_infec;} 
+1 { (0,0)~initial != -1 }
+
+
+
+
+
+
+
 
 

epidemic_hoya_et_al/SIR_Hoya.pal

@@ -0,0 +1,7 @@
+[-100; -0.001] 255 0 0
+[0; 0.200] 0 0 0
+[0.201; 0.400] 50 50 50
+[0.401; 0.600] 100 100 100
+[0.601; 0.800] 150 150 150
+[0.801; 1.00] 200 200 200
+[1.001;100] 255 0 0

epidemic_hoya_et_al/SIR_Hoya.var

@@ -0,0 +1 @@
+(25,25) = 10 0.6 1 1 0.7 0.0 0.3

epidemic_hoya_et_al/SIR_Hoya_macros.inc

@@ -1,5 +1,5 @@
 #BeginMacro(e) 
-0.5
+0.4
 #EndMacro