Zhong model

Not finished, I will debug it tomorrow

Author: khvilaboa

epidemic_zhong_et_al/sirs_zhong.ma

@@ -0,0 +1,57 @@
+#include(sirs_zhong_macros.inc)
+[top]
+components : area
+
+[area]
+type : cell
+width : 50
+height : 50
+delay : transport
+defaultDelayTime : 1
+border : wrapped
+neighbors : 			area(-1,0) 
+neighbors : area(0,-1)  area(0,0)  area(0,1)
+neighbors : 			area(1,0) 
+initialvalue : -1
+
+
+localtransition : sirs-zhong-rule
+
+statevariables: population area p contact_rate cured_rate t_i t_p t_l t_r i_sus_0 i_inc_1 i_inc_2 i_inf_3 i_inf_4 i_lat_5 i_lat_6 i_rec_7 i_rec8 i_rec9
+statevalues: 100 1 1000 0.5 0.02 2 2 2 3 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+initialvariablesvalue: sirs_zhong.var
+
+neighborports : initial sus_0 inc_1 inc_2 inf_3 inf_4 lat_5 lat_6 rec_7 rec_8 rec_9 pop
+
+[sirs-zhong-rule]
+rule : {~initial := 0; ~pop := $population; ~sus_0 := $i_sus_0; ~inc_1 := $i_inc_1; ~inc_2 := $i_inc_2; ~inf_3 := $i_inf_3; ~inf_4 := $i_inf_4; ~lat_5 := $i_lat_5; ~lat_6 := $i_lat_6; ~rec_7 := $i_rec_7; ~rec_8 := $i_rec_8; ~rec_9 := $i_rec_9;}
+ 1 {(0,0)~initial = -1}
+
+rule : {
+~pop := $population; ~sus_0 := $i_sus_0; ~inc_1 := $i_inc_1; ~inc_2 := $i_inc_2; ~inf_3 := $i_inf_3; ~inf_4 := $i_inf_4; ~lat_5 := $i_lat_5; ~lat_6 := $i_lat_6; ~rec_7 := $i_rec_7; ~rec_8 := $i_rec_8; ~rec_9 := $i_rec_9;
+}
+{
+$i_rec_9 := $i_rec_8;
+$i_rec_8 := $i_rec_9;
+$i_rec_7 := $i_lat_6 + #macro(local_cured);
+
+$i_lat_6 := trunc(min((1 - $cured_rate) * $i_lat_5, 1)*100)/100;
+$i_lat_5 := trunc(min((1 - $cured_rate) * $i_inf_4, 1)*100)/100;
+$i_inf_4 := trunc(min((1 - $cured_rate) * $i_inf_3, 1)*100)/100;
+$i_inf_3 := trunc(min((1 - $cured_rate) * $i_inc_2, 1)*100)/100;
+$i_inc_2 := trunc(min((1 - $cured_rate) * $i_inc_1, 1)*100)/100;
+$i_inc_1 := #macro(internal_infected) + #macro(external_infected);
+
+$i_sus_0 := 1 - $i_inc_1 - $i_inc_2 - $i_inf_3 - $i_inf_4 - $i_lat_5 - $i_lat_6 - $i_rec_7 - $i_rec8 - $i_rec9;
+
+}
+1 { (0,0)~initial != -1 }
+
+
+
+
+
+
+
+
+

epidemic_zhong_et_al/sirs_zhong.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_zhong_et_al/sirs_zhong.var

@@ -0,0 +1 @@
+(25,25) = 100 1 1000 0.5 0.02 2 2 2 3 0.7 0.3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0

epidemic_zhong_et_al/sirs_zhong_macros.inc

@@ -0,0 +1,46 @@
+
+#BeginMacro(local_cured) 
+(( trunc(min($cured_rate * $i_lat_5, 1)*100) + 
+trunc(min($cured_rate * $i_inf_4, 1)*100) + 
+trunc(min($cured_rate * $i_inf_3, 1)*100) + 
+trunc(min($cured_rate * $i_inc_2, 1)*100) + 
+trunc(min($cured_rate * $i_inc_1, 1)*100) ) / 100 )
+#EndMacro
+
+#BeginMacro(internal_infected) 
+( truncate(min(( $contact_rate * 
+((%i_sus_0 * $population) / $area) *
+%i_sus_0 * 
+($i_lat_6 + $i_lat_5 + $i_inf_4 + $i_inf_3 + $i_incv_2 + $i_inc_1)
+) / $p, 1)*100) / 100 )
+#EndMacro
+
+#BeginMacro(external_infected) 
+( truncate(min(( $connection *
+$contact_rate * 
+((%i_sus_0 * $population) / $area) *
+%i_sus_0 * 
+((1,0)~lat_6 + (1,0)~lat_5 + (1,0)~inf_4 + (1,0)~inf_3 + (1,0)~inc_2 + (1,0)~inc_1)
+) / $p, 1)*100) / 100 )
++
+( truncate(min(( $connection *
+$contact_rate * 
+((%i_sus_0 * $population) / $area) *
+%i_sus_0 * 
+((0,-1)~lat_6 + (0,-1)~lat_5 + (0,-1)~inf_4 + (0,-1)~inf_3 + (0,-1)~inc_2 + (0,-1)~inc_1)
+) / $p, 1)*100) / 100 )
++
+( truncate(min(( $connection *
+$contact_rate * 
+((%i_sus_0 * $population) / $area) *
+%i_sus_0 * 
+((0,1)~lat_6 + (0,1)~lat_5 + (0,1)~inf_4 + (0,1)~inf_3 + (0,1)~inc_2 + (0,1)~inc_1)
+) / $p, 1)*100) / 100 )
++
+( truncate(min(( $connection *
+$contact_rate * 
+((%i_sus_0 * $population) / $area) *
+%i_sus_0 * 
+((-1,0)~lat_6 + (-1,0)~lat_5 + (-1,0)~inf_4 + (-1,0)~inf_3 + (-1,0)~inc_2 + (-1,0)~inc_1)
+) / $p, 1)*100) / 100 )
+#EndMacro