---russianmorbidity.nlogo

NETLOGO

; breeds [race_one race_two]

globals [
time
Total_Healthy Total_Normal Total_Unhealthy Total_Kids Total_Adults
Percent_Healthy Percent_Normal Percent_Unhealthy Percent_Kids      
Total_Healthy_One Total_Normal_One Total_Unhealthy_One
Total_Healthy_Two Total_Normal_Two Total_Unhealthy_Two
Total_Healthy_Three Total_Normal_Three Total_Unhealthy_Three
Percent_Healthy_One Percent_Normal_One Percent_Unhealthy_One
Percent_Healthy_Two Percent_Normal_Two Percent_Unhealthy_Two
Percent_Healthy_Three Percent_Normal_Three Percent_Unhealthy_Three
number
default ; patch color default
youngshade oldshade
mean_age mean_expected_life mean_lifespan
shock_all? shock_1? shock_2? shock_3? shock_4?

region_1_healthy region_1_Normal region_1_unhealthy region_1_kids 
region_2_healthy region_2_Normal region_2_unhealthy region_2_kids 
region_3_healthy region_3_Normal region_3_unhealthy region_3_kids 
region_4_healthy region_4_Normal region_4_unhealthy region_4_kids 

region_1_total region_2_total region_3_total region_4_total

mean_expected_life1 mean_expected_life2 mean_expected_life3 mean_expected_life4
mean_lifespan1 mean_lifespan2 mean_lifespan3 mean_lifespan4
        ]      

turtles-own [ 
age
generation
status
income
random_death
Expected_Life
region_id

Exp_Utility_Healthy
Exp_Utility_Normal
Exp_Utility_Unhealthy 
lifestyle

          ]
            
patches-own [hood_lifespan death_total sum_death_age border? no_data? ]


to setup
  ca
  set-default-shape turtles "med_circle"
  set number ((density / 100) * count patches)
  regions
  ask random-n-of (number) patches with [border? = false]
    [ sprout 1 [ set status 65 set color status set age random 60 set generation 1 ] ]
  ask turtles [ ifelse age < 16 [set lifestyle 0][set lifestyle 2 set status 105 set color status] region-id]
  set time 0
  set shock_all? false set shock_1? false set shock_2? false set shock_3? false set shock_4? false   
  set mean_lifespan1 0 set mean_lifespan2 0 set mean_lifespan3 0 set mean_lifespan4 0
end

to regions

ask patches [
set border? false
if N-Regions = 2 [ if pycor = 0 [set pcolor 9.9999 set border? true]]
if N-Regions = 3 [ if pycor = 0 or (pycor < 1 and pxcor = 0) [set pcolor 9.9999 set border? true]]
if N-Regions = 4 [if pycor = 0 or pxcor = 0 [set pcolor 9.9999 set border? true]]
            ]               
end

to go
  set time time + 1
  ask turtles [without-interruption [

    wander
    region-id
    reproduce
                                    ]]
  ask patches [ update-patches ]

  ask turtles [without-interruption [

    kids-grow-old
    Healthy-grow-old
    Normal-grow-old
    Unhealthy-grow-old
    
                  ]]
  do-plots
  stats
  regional_stats
  regional-mean-life
  shock
  if count turtles = 0
    [ stop ]
end

to region-id

  if xcor > 0 and ycor > 0 [set region_id 1]
  if xcor < 0 and ycor > 0 [set region_id 2]
  if xcor < 0 and ycor < 0 [set region_id 3]
  if xcor > 0 and ycor < 0 [set region_id 4]
                     
end

to kids-grow-old
  if lifestyle = 0 [
  set age age + 1
  ;;choose - shading or no shading
  ;; set color scale-color status age oldshade youngshade
  set color status
  ;set expected_life mean values-from neighbors with [border? = false] [hood_lifespan]
  ifelse (count neighbors with [border? = false and no_data? = false] > 0)
          [set expected_life mean values-from neighbors with [border? = false and no_data? = false ] [hood_lifespan]]
             [set expected_life def_life_exp]

  set random_death random 1000
  if ( Random_Death < (10 * Probability_Death_Kid + 1) )[ set death_total death_total + 1 set sum_death_age sum_death_age + age die]

  if age = 16 [ 
     choose-lifestyle
     ifelse sum values-from neighbors [count turtles-here with [lifestyle > 0]] < count neighbors with [border? = false] [go-if-empty-nearby][go-if-empty]
              ]
                    ]
end

to Healthy-grow-old
 if lifestyle = 1 [
  set age age + 1
  ;;choose - shading or no shading
  ;; set color scale-color status age oldshade youngshade
  set color status
  set random_death random 1000
  if (age < 50) and ( Random_Death < 10 * Probability_Death_Healthy )[ set death_total death_total + 1 set sum_death_age sum_death_age + age die ]
  if (age > 50) and (age < 90) and ( Random_Death < (20 + 10 * Probability_Death_Healthy) )[ set death_total death_total + 1 set sum_death_age sum_death_age + age die]
  if (age = 90 or age > 90)
          [ set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]
  set expected_life mean values-from neighbors [hood_lifespan]
  
                         ]
end

to Normal-grow-old
 if lifestyle = 2 [
  set age age + 1
  ;;choose - shading or no shading
  ;; set color scale-color status age oldshade youngshade
  set color status
  set random_death random 1000
  if (age < 50) and ( Random_Death < 10 * Probability_Death_Normal )[ set death_total death_total + 1 set sum_death_age sum_death_age + age die ]
  if (age > 50) and (age < 90) and ( Random_Death < (20 + 10 * Probability_Death_Normal) )[ set death_total death_total + 1 set sum_death_age sum_death_age + age die ]
  if (age = 90 or age > 90)
          [ set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]
  
  
                  ]
end

to Unhealthy-grow-old
 
 if lifestyle = 3 [
  set age age + 1
  ;;choose - shading or no shading
  ;; set color scale-color status age oldshade youngshade
  set color status
  set random_death random 1000
  if (age < 50) and ( Random_Death < 10 * Probability_Death_Unhealthy )[ set death_total death_total + 1 set sum_death_age sum_death_age + age die ]
  if (age > 50) and (age < 90) and ( Random_Death < (20 + 10 * Probability_Death_Unhealthy) )[ set death_total death_total + 1 set sum_death_age sum_death_age + age die ]
  if (age = 90 or age > 90)
          [ set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]
  
  
                     ]
end


to wander
locals [random_movement]

set random_movement random 100
if random_movement < movement [
  go-if-empty
                              ]
end

to reproduce
 
  if (age > 16 and age < 45 and random 1000 < repro_rate) [
    hatch 1 [
          set age 0 set generation generation + 1 set lifestyle 0 set status 65
          set color status 
       ;; ifelse sum values-from neighbors [count turtles-here] < 8 [go-if-empty-nearby][go-if-empty]
            ]
                                              ]
end

to go-if-empty  ;; turtle procedure
  locals [empty-patches target]
  without-interruption
    [ set empty-patches patches with [not any? turtles-here and border? = false]
      if any? empty-patches
        [ set target random-one-of empty-patches
          set heading towards target
          fd distance target ] 
    region-id
    ]
end

to go-if-empty-nearby  ;; turtle procedure
  locals [empty-patches target]
  without-interruption
    [ set empty-patches neighbors with [not any? turtles-here and border? = false]
      if any? empty-patches
        [ set target random-one-of empty-patches
          set heading towards target
          fd distance target ] 
    region-id
    ]
end

to choose-lifestyle

;set Exp_Utility_Healthy (((Expected_Life ) * b_healthy ) + a_healthy)
;set Exp_Utility_Normal (((Expected_Life ) * b_normal ) + a_normal)
;set Exp_Utility_Unhealthy (((Expected_Life ) * b_unhealthy ) + a_unhealthy)

ifelse discounting? = false 
    [ set Exp_Utility_Healthy (((Expected_Life - 16) *  b_healthy) + a_healthy)
      set Exp_Utility_Normal (((Expected_Life - 16) * b_normal) + a_normal) 
      set Exp_Utility_Unhealthy (((Expected_Life - 16) * b_unhealthy) + a_unhealthy)
    ]
    [ set Exp_Utility_Healthy (a_healthy + ((b_healthy * (1 - exp(- rate * (Expected_Life - 16 )))) / (exp(rate) - 1)))
      set Exp_Utility_Normal (a_normal + ((b_normal * (1 - exp(- rate * (Expected_Life - 16)))) / (exp(rate) - 1)))   
      set Exp_Utility_Unhealthy (a_unhealthy + ((b_unhealthy * (1 - exp(- rate * (Expected_Life - 16)))) / (exp(rate) - 1)))
    ]

ifelse (Exp_Utility_Healthy >= Exp_Utility_Normal) and (Exp_Utility_Healthy >= Exp_Utility_Unhealthy and Healthy_Option = true)
       [set lifestyle 1 set status 89]
       [ ifelse (Exp_Utility_Normal >= Exp_Utility_Unhealthy ) [set lifestyle 2 set status 105][ set lifestyle 3 set status 26]]    

end

to update-patches
set no_data? false 
ifelse death_total > 0 [set hood_lifespan (sum_death_age / death_total)] [set no_data? true set hood_lifespan mean_expected_life]


end

to shock

if timedshock1 = true and time > 39 and time < 50 [set shock_1? true]

if shock_all? [ ask random-n-of shock_size turtles 
              [ set shock_all? false set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]  ]

if shock_1? [ ifelse N-regions = 2 or N-regions = 3 [ ask random-n-of shock_size turtles with [ycor > 0]
              [ set shock_1? false set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]  ]
                                                   [ ask random-n-of shock_size turtles with [ycor > 0 and xcor > 0]
              [ set shock_1? false set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]  ]
            ]
if shock_2?  [ ifelse N-regions = 2 or N-regions = 3 [ ask random-n-of shock_size turtles with [ycor < 0]
              [ set shock_2? false set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]  ]
                                                    [ ask random-n-of shock_size turtles with [ycor > 0 and xcor < 0]
              [ set shock_2? false set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]  ]
            ]
if shock_3?  [ ifelse N-regions = 3 [ ask random-n-of shock_size turtles with [ycor < 0 and xcor > 0]
              [ set shock_3? false set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]  ]
                                   [ ask random-n-of shock_size turtles with [ycor < 0 and xcor < 0]
              [ set shock_3? false set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]  ]
            ]
if shock_4?  [ if N-regions = 4 [ ask random-n-of shock_size turtles with [ycor < 0 and xcor > 0]
              [ set shock_4? false set death_total death_total + 1 set sum_death_age sum_death_age + age set pcolor default die ]  ]
            ]                                                                       
end


to regional_stats

set region_1_kids count turtles with [lifestyle = 0 and region_id = 1] 
set region_1_healthy count turtles with [lifestyle = 1 and region_id = 1]
set region_1_Normal count turtles with [lifestyle = 2 and region_id = 1]
set region_1_unhealthy count turtles with [lifestyle = 3 and region_id = 1]
set mean_expected_life1 mean values-from turtles with [region_id = 1] [expected_life]

set region_2_kids count turtles with [lifestyle = 0 and region_id = 2] 
set region_2_healthy count turtles with [lifestyle = 1 and region_id = 2]
set region_2_Normal count turtles with [lifestyle = 2 and region_id = 2]
set region_2_unhealthy count turtles with [lifestyle = 3 and region_id = 2]
set mean_expected_life2 mean values-from turtles with [region_id = 2] [expected_life]

set region_3_kids count turtles with [lifestyle = 0 and region_id = 3] 
set region_3_healthy count turtles with [lifestyle = 1 and region_id = 3]
set region_3_Normal count turtles with [lifestyle = 2 and region_id = 3]
set region_3_unhealthy count turtles with [lifestyle = 3 and region_id = 3]
set mean_expected_life3 mean values-from turtles with [region_id = 3] [expected_life]

set region_4_kids count turtles with [lifestyle = 0 and region_id = 4] 
set region_4_healthy count turtles with [lifestyle = 1 and region_id = 4]
set region_4_Normal count turtles with [lifestyle = 2 and region_id = 4]
set region_4_unhealthy count turtles with [lifestyle = 3 and region_id = 4]
set mean_expected_life4 mean values-from turtles with [region_id = 4] [expected_life]

set region_1_total count turtles with [region_id = 1]
set region_2_total count turtles with [region_id = 2]
set region_3_total count turtles with [region_id = 3]
set region_4_total count turtles with [region_id = 4]

end

to regional-mean-life

if time > 20 [

set mean_lifespan1 (sum values-from patches with [pxcor > 0 and pycor > 0] [sum_death_age] / sum values-from patches with [pxcor > 0 and pycor > 0] [death_total])
set mean_lifespan2 (sum values-from patches with [pxcor < 0 and pycor > 0] [sum_death_age] / sum values-from patches with [pxcor < 0 and pycor > 0] [death_total])
set mean_lifespan3 (sum values-from patches with [pxcor < 0 and pycor < 0] [sum_death_age] / sum values-from patches with [pxcor < 0 and pycor < 0] [death_total])
set mean_lifespan4 (sum values-from patches with [pxcor > 0 and pycor < 0] [sum_death_age] / sum values-from patches with [pxcor > 0 and pycor < 0] [death_total])

]

end


to stats

ifelse time > 2 [
set mean_age mean values-from turtles [age]
set mean_expected_life mean values-from turtles [expected_life]
            ] [set mean_expected_life def_life_exp]
if time > 20 [set mean_lifespan (sum values-from patches [sum_death_age] / sum values-from patches [death_total])] 
             
set Total_Healthy count turtles with [lifestyle = 1]
set Total_Normal count turtles with [lifestyle = 2]
set Total_Unhealthy count turtles with [lifestyle = 3]
set Total_Kids count turtles with [lifestyle = 0]
set Total_Adults count turtles with [lifestyle > 0]
             
set percent_kids (Total_Kids / count turtles) set percent_healthy (Total_healthy / Total_Adults)
set percent_normal (Total_normal / Total_Adults) set percent_unhealthy (Total_unhealthy / Total_Adults)
             
end

to do-plots
 
  set-current-plot "Mean Lifespan and Expected Lifespan"
      set-current-plot-pen "Mean Lifespan"
      plot mean_lifespan
      set-current-plot-pen "Expected Lifespan"
      plot mean_expected_life
  set-current-plot "Unhealthy"
      set-current-plot-pen "Unhealthy 1"
      plot region_1_unhealthy 
      set-current-plot-pen "Unhealthy 2"
      plot region_2_unhealthy
      set-current-plot-pen "Unhealthy 3"
      plot region_3_unhealthy 
      set-current-plot-pen "Unhealthy 4"
      plot region_4_unhealthy  
      
 end
@#$#@#$#@
GRAPHICS-WINDOW
261
10
702
472
20
20
10.5122
1
20
0
0
0
0

CC-WINDOW
5
720
887
815
Command Center

BUTTON
8
10
93
43
Setup
setup
NIL
1
T
OBSERVER
T
NIL

MONITOR
98
11
182
60
total population
count turtles
3
1

BUTTON
7
44
93
77
Go/Pause
go
T
1
T
OBSERVER
NIL
NIL

SLIDER
704
227
877
260
Probability_Death_Healthy
Probability_Death_Healthy
0
25
1.0
0.1
1
NIL

SLIDER
7
80
150
113
density
density
0
100
70
1
1
NIL

MONITOR
185
11
242
60
NIL