In [58]:

    

path = joinpath(Pkg.dir(),"Sugarscape", "src", "chapterII.jl")


    

    
Out[58]:





"C:\\Users\\Ken\\.julia\\v0.3\\Sugarscape\\src\\chapterII.jl"

In [59]:

    

include(path)


    

    




Warning: replacing module ChapterII

In [60]:

    

sc = ChapterII


    

    
Out[60]:





ChapterII

In [61]:

    

sc.run(1)
@time scape = sc.run(500)


    

    




elapsed time: 0.21361148 seconds (34653772 bytes allocated, 13.46% gc time)






    
Out[61]:





Sugarscape with 400 agents

In [62]:

    

sc.plot(scape)


    

    













    
Out[62]:





PyObject <matplotlib.text.Text object at 0x000000003A5E15C0>

In [63]:

    

function run_carrypop(N; kwargs...)
    carrypop = Int[]
    scape = sc.init_scape(sc.init_capacity(); kwargs...)    
    push!(carrypop, sum(map(sc.alive,scape.agents)))
    
    for cnt = 1:N
        sc.timestep(scape; kwargs...)
        push!(carrypop, sum(map(sc.alive,scape.agents)))
    end
    scape,carrypop
end


    

    
Out[63]:





run_carrypop (generic function with 1 method)

In [64]:

    

scape, carrypop = run_carrypop(500)
PyPlot.plot(carrypop)
PyPlot.title("Carrying capacity")
PyPlot.xlabel("time steps")


    

    













    
Out[64]:





PyObject <matplotlib.text.Text object at 0x0000000039E780B8>

In [65]:

    

using DataFrames


    

In [69]:

    

function carrying(;sim_steps=3, time_steps=500)
    carry = DataFrame(vision=Int[], metabolism=Int[], population = Int[])
    
    for vision = 1:10, metabolism = 1:3, sims = 1:sim_steps
        scape = sc.run(time_steps; init_vision=()->rand(1:(2vision-1)), init_metabolism=()->rand(1:(2metabolism-1)))
        push!(carry, [vision, metabolism, sum(map(sc.alive, scape.agents))])
    end
    carry
end


    

    
Out[69]:





carrying (generic function with 1 method)

In [70]:

    

carry = carrying();


    

    
Out[70]:




vision
metabolism
population
1
1
1
359
2
1
1
360
3
1
1
364
4
1
2
184
5
1
2
189
6
1
2
183
7
1
3
104
8
1
3
138
9
1
3
138
10
2
1
378
11
2
1
380
12
2
1
375
13
2
2
217
14
2
2
216
15
2
2
220
16
2
3
159
17
2
3
146
18
2
3
149
19
3
1
388
20
3
1
381
21
3
1
380
22
3
2
227
23
3
2
243
24
3
2
226
25
3
3
153
26
3
3
158
27
3
3
149
28
4
1
387
29
4
1
384
30
4
1
386
&vellip
&vellip
&vellip
&vellip

In [92]:

    

using Gadfly


    

In [104]:

    

carry_aggr = by(carry, [:vision, :metabolism],df -> DataFrame(pop = mean(df[:population])))
plot(carry_aggr, x=:vision, y=:pop, color=:metabolism, Geom.line)


    

    
Out[104]:




vision
metabolism
pop
1
1
1
361.0
2
1
2
185.33333333333334
3
1
3
126.66666666666667
4
2
1
377.6666666666667
5
2
2
217.66666666666666
6
2
3
151.33333333333334
7
3
1
383.0
8
3
2
232.0
9
3
3
153.33333333333334
10
4
1
385.6666666666667
11
4
2
233.33333333333334
12
4
3
157.33333333333334
13
5
1
391.0
14
5
2
239.0
15
5
3
154.0
16
6
1
393.6666666666667
17
6
2
259.3333333333333
18
6
3
156.33333333333334
19
7
1
393.0
20
7
2
247.0
21
7
3
159.66666666666666
22
8
1
395.0
23
8
2
254.0
24
8
3
159.0
25
9
1
395.6666666666667
26
9
2
249.66666666666666
27
9
3
166.66666666666666
28
10
1
397.6666666666667
29
10
2
248.66666666666666
30
10
3
156.0
&vellip
&vellip
&vellip
&vellip

In [108]:

    

plot(carry_aggr, x=:vision, y=:pop, color=:metabolism, Geom.line)


    

    
Out[108]:

In [112]:

    

import StatsBase


    

In [115]:

    

StatsBase.hist(scape.agents,10)


    

    




`histrange` has no method matching histrange(::Array{AgentInfo,1}, ::Int64)
while loading In[115], in expression starting on line 1

 in hist at statistics.jl:598