Things we need to do:
1) Make our set of randomized environment files:
These files will consist of a set of values
where sum(abcdefg)=7, 0=h=i
REACTION NOT not process:value=a:type=add requisite:max_count=1
REACTION NAND nand process:value=b:type=add requisite:max_count=1
REACTION AND and process:value=c:type=add requisite:max_count=1
REACTION ORN orn process:value=d:type=add requisite:max_count=1
REACTION OR or process:value=e:type=add requisite:max_count=1
REACTION ANDN andn process:value=f:type=add requisite:max_count=1
REACTION NOR nor process:value=g:type=add requisite:max_count=1
REACTION XOR xor process:value=h:type=add requisite:max_count=1
REACTION EQU equ process:value=i:type=add requisite:max_count=1
2) Set up 4 events files (python script to generate events files)
1-(Control A) Constant environment where 1=a=b=c=d=e=f, 0=h=i, every 500 updates kill off 95% of population
2-(Control B) Constant environment where 1=a=b=c=d=e=f, 0=h=i, every 5,000 updates kill off 95% of population
(these second two need to be x50, with different random environments)
3-(Treatment A) Changing environment, every 500 updates, change to random environment.cfg, kill off 95% of population
4-(Treatment B) Changing environment, every 5000 updates, change to random environment.cfg, kill off 95% of population
IN ALL EVENTS FILES at 250,000 updates switch to environment file: where 0=a=b=c=d=e=f=g=h, 10=i .
allow to run for 50,000 more updates, log time EQU evolves
3) Set up avida.cfg
100 x 100 grid
Decrease mutation rate to 0.02
run for 300,000 updates
4) Run all 4 treatments x50.
5) Analysis!!