Benchmark Model5

"Model 5 = Model4 + metapopulation"
KendrickModel SIRDemography
	attribute: #(status -> S I R);
	parameters: #(beta lambda gamma b mu);
	transitions: #(
		S -- lambda --> I.
		I -- gamma --> R.
		status -- mu --> Empty.
		Empty -- b --> S.
	);
	lambda: #(beta*I/N).
	
KendrickModel AgeStructure
	attribute: #(age -> J A);
	parameters: #(m);
	transitions: #(
		J -- m --> A.
	).
	
"Host: H, Vector: V"
KendrickModel Species
	attribute: #(species -> H V).
	
KendrickModel MetaPopulation
	attribute: #(zone -> 1 to: 20);
	parameters: #(rho);
	connectionTopo: #ring.
	
Composition Model5
	model: 'SIRDemography';
	model: 'AgeStructure';
	model: 'Species';
	model: 'MetaPopulation'.
	
Scenario M5Parameters
	 on: 'Model5';
	 populationSize: 100;
    beta_species_age: #(
		#(#(0.1 0.01) #(0.01 0.01))
		#(#(0.01 0.1) #(0.06 0.06))
		); 
    gamma_species: {1/30. 1/20}; 
	 mu_age: {0.01. 0.0025};
	 m: 0.02;
	 b_species_age: #(#(0.005 0) #(0 0.05));
	 lambda: #(beta*I_species/N sum);
	 rho_species: #(0.02 0.0);
	 S_species: #(3000 900);
	 I_species: #(0 100).
	
Simulation M5Sim rungeKutta
	scenarios: #(M5Parameters);
	from: 0.0; 
	to: 730; 
	step: 1.
	
Visualization M4Vis diagram 
	for: 'M5Sim';
	data: #(I_species);
	xLabel: 'Time (days)';
	exportToPng.