poblacion

globals [cont blues blacks greens reds time]

to setup
  clear-all
    ask n-of (S)patches
  [
    sprout 1
    [
      set shape "person"
      set color blue
    ]
  ]
    ask n-of (I)patches
  [
    sprout 1
    [
      set shape "person"
      set color red
    ]
  ]
    ask n-of (R)patches
  [
    sprout 1
    [
      set shape "person"
      set color green
    ]
  ]
end 

to go
  ask turtles with [color = blue]
  [
  set cont 0
    ask neighbors
    [
      set cont cont + count turtles-here with [color = red]
    ]
      if cont != 0 and random-float 1.0 < alpha
      [
        set color red
      ]
  ]
  ask turtles with [color = red]
  [
    if random-float 1.0 < beta
    [
      set color green
    ]
    if random-float 1.0 < gamma
    [
      set color black
    ]
  ]
  set blacks count turtles with [color = black]
  set greens count turtles with [color = green]
  set reds count turtles with [color = red]
  set blues count turtles with [color = blue]
  set time time + 1
  set-current-plot "poblacion"
  set-current-plot-pen "blues"
  plotxy time blues
  set-current-plot-pen "reds"
  plotxy time reds
  set-current-plot-pen "greens"
  plotxy time greens
  set-current-plot-pen "blacks"
  plotxy time blacks
end