aaa---smart-diverse.nlogo

NETLOGO

globals
[
  generations          ; counts how many generations have passed
  red-count            ; population of red turtles
  blue-count           ; population of blue turtles
;  red-lag              ; lagged red population
;  blue-lag             ; lagged blue population
  parent-type          ; type of my parent
  environ              ; integer environ indicator
  time-in-environ      ; how long we stay in this environ
]

turtles-own
[
  my-type              ; turtle-type (to match with environ)
  time-til-repro       ; number of time steps remaining until reproduction (-1 means not counting down)
]

to setup
  clear-output
  setup-experiment
end

to setup-experiment
  cp ct
  clear-all-plots
  set generations 0
  set environ 1 + random environ-max
  set time-in-environ 1 + random time-in-environ-max
  cct carrying-capacity / 10
  [
    setxy (random-float screen-size-x)
          (random-float screen-size-y)    ; randomize turtle locations
    ifelse who < (carrying-capacity / 20)  ; start out with equal numbers of reds and blues
      [ set color blue ]
      [ set color red  ]
      set my-type 1 + random environ-max
      set time-til-repro -1
  ]
;  set red-lag [0 0 0 0 0 0 0 0 0 0]
;  set blue-lag [0 0 0 0 0 0 0 0 0 0]
  setup-plot
  plot-counts
end

to go
  wander
  live
  grim-reaper
  set generations generations + 1
  set time-in-environ time-in-environ - 1
  if time-in-environ = 0
    [ set environ 1 + random environ-max
      set time-in-environ 1 + random time-in-environ-max
    ]
;  set red-lag fput count turtles with [ color = red ] but-last red-lag
;  set blue-lag fput count turtles with [ color = blue ] but-last blue-lag
  plot-counts
end

;; to enable many repetitions with same settings
to go-experiment
  go
  if red-count < 10
  [
    output-print "red extinct after " + generations + " generations"
    setup-experiment
  ]
  if blue-count < 10
  [
    output-print "blue extinct after " + generations + " generations"
    setup-experiment
  ]
end

to wander
  ask turtles
  [
    rt random-float 30 - random-float 30
    fd 5
  ]
end

to live
  ask turtles
  [
    ifelse color = red
    [
      if my-type < environ
        [ set my-type my-type + learning-rate ]
      if my-type > environ
        [ set my-type my-type - learning-rate ]
      ifelse time-til-repro < 0
        [ if (abs (environ - my-type)) < 2
          [ set time-til-repro red-gestation ]
        ]
        [ if time-til-repro = 0
            [ set parent-type my-type
              hatch 1
                [ set color red
                  ifelse Memes
                    [ set my-type round random-normal parent-type red-spread ]
                    [ set my-type 1 + random environ-max ]
                  set time-til-repro -1
                  setxy (random-float screen-size-x)
                        (random-float screen-size-y)    ; randomize turtle locations
                ]
             ]
           set time-til-repro time-til-repro - 1
         ]
      ]
      ; blue:
      [ ifelse time-til-repro < 0
        [ set time-til-repro 1 + round ( 0.75 * min list (abs environ - my-type) (abs my-type - environ) ) ]
        [ if time-til-repro = 0
          [ set parent-type my-type
              hatch 1
                [ set color blue
                  ; set my-type parent-type + (3 - random 7)
                  set my-type 1 + random environ-max
                  set time-til-repro -1
                  setxy (random-float screen-size-x)
                        (random-float screen-size-y)    ; randomize turtle locations
                ]
           ]
           set time-til-repro time-til-repro - 1
         ]
       ]
   ]
end

;; kill turtles in excess of carrying capacity
;; note that reds and blues have equal probability of dying
to grim-reaper
  let num-turtles count turtles
  if num-turtles <= carrying-capacity
    [ stop ]
  let chance-to-die (num-turtles - carrying-capacity) / num-turtles
  ask turtles
  [
    if random-float 1.0 < chance-to-die
      [ die ]
  ]
end

to plot-counts
  set-current-plot "Populations"

  set red-count count turtles with [ color = red ]
  set-current-plot-pen "Reds"
  plot red-count

  set blue-count count turtles with [ color = blue ]
  set-current-plot-pen "Blues"
  plot blue-count

  set-current-plot-pen "Total"
  plot count turtles
  
  set-current-plot "Type Distributions"
  
  set-current-plot-pen "my-type-red"
  histogram-from turtles with [ color = red ] [my-type]
  
  set-current-plot-pen "my-type-blue"
  histogram-from turtles with [ color = blue ] [my-type]
  
;  if generations > 60
;    [ set-current-plot "Lagged"
;  
;      set-current-plot-pen "Red-lag"
;      plotxy item Lag red-lag item 0 red-lag
;  
;      set-current-plot-pen "Blue-lag"
;      plotxy item Lag blue-lag item 0 blue-lag
;    ]
  
end


to setup-plot
  set-current-plot "Populations"
  set-plot-y-range 0 floor ((carrying-capacity / environ-max ) * 1.2)

  set-current-plot "Type Distributions"
  set-plot-x-range 0 environ-max
  set-plot-y-range 0 round (10 * (count turtles) / environ-max)
  set-histogram-num-bars environ-max
  
;  set-current-plot "Lagged"
;  set-plot-x-range 0 1.025 * carrying-capacity
;  set-plot-y-range 0 1.025 * carrying-capacity
end



; *** NetLogo Model Copyright Notice ***
;
; Portions of the code for this model were derived from a Simple
; Birth Rates Model.
;
; That model was created as part of the project: CONNECTED MATHEMATICS:
; MAKING SENSE OF COMPLEX PHENOMENA THROUGH BUILDING OBJECT-BASED PARALLEL
; MODELS (OBPML).  The project gratefully acknowledges the support of the
; National Science Foundation (Applications of Advanced Technologies
; Program) -- grant numbers RED #9552950 and REC #9632612.
;
; Copyright 1998 by Uri Wilensky. All rights reserved.
;
; Permission to use, modify or redistribute this model is hereby granted,
; provided that both of the following requirements are followed:
; a) this copyright notice is included.
; b) this model will not be redistributed for profit without permission
;    from Uri Wilensky.
; Contact Uri Wilensky for appropriate licenses for redistribution for
; profit.
;
;
; To refer to this model in academic publications, please use:
; Carter, T. (2005).  NetLogo Smart-Diverse model.
;
;
; *** End of NetLogo Model Copyright Notice ***
@#$#@#$#@
GRAPHICS-WINDOW
290
10
603
344
50
50
3.0
0
10
1
1
1
0
1
1
1

CC-WINDOW
5
478
960
573
Command Center
0

BUTTON
133
10
254
43
run-experiment
go-experiment
T
1
T
OBSERVER
T
NIL

SLIDER
9
46
246
79
learning-rate
learning-rate
0
1.0
0.33
0.01
1
NIL

BUTTON
5
10
65
43
setup
setup
NIL
1
T
OBSERVER
T
NIL

BUTTON
69
10
129
43
go
go
T
1
T
OBSERVER
T
NIL

PLOT
4
245
284
464
Populations
Generations
Population
0.0
50.0
0.0
1200.0
true
true
PENS
"Reds" 1.0 0 -2674135 true
"Blues" 1.0 0 -13345367 true
"Total" 1.0 0 -10899396 true

MONITOR
8
193
84
242
Num Reds
red-count
3
1

MONITOR
93
193
170
242
Num Blues
blue-count
3
1

MONITOR
177
193
267
242
Generations
generations
3
1

OUTPUT
291
351
603
445

SLIDER
9
118
155
151
environ-max
environ-max
0
100
70
5
1
NIL

SLIDER
159
118
286
151
time-in-environ-max
time-in-environ-max
1
100
50
1
1
NIL

SLIDER
9
82
154
115
carrying-capacity
carrying-capacity
0
20000
4000
1000
1
NIL

PLOT
632
68
951
300
Type Distributions
Types
Population
1.0
70.0
0.0
100.0
true