NetLogo Models Library:
This model demonstrates some connections between probability and the natural sciences. Specifically, the model uses combinatorial space, sampling, and distribution in a genotype/phenotype analysis of fish procreation. The model allows you to look "under the hood": you can study a Mendel-type visualization of the combinations of dominant and recessive genes that underlie changes and trends in genetic distribution.
Fish vary by body and fin color, each of which can be either green or blue, so there are four different phenotypes. Underlying this variation are dominant and recessive genes, expressed in "4-blocks." A 4-block is a 2-by-2 array of squares, each of which can be either green or blue. Initially, fish are randomly distributed by 4-block genotype. Then, fish seek mates and reproduce, and the genetic distribution changes.
This model is a part of the ProbLab curriculum. For more information about the ProbLab curriculum please refer to http://ccl.northwestern.edu/curriculum/ProbLab/.
Fishes' Mendel-type dominant/recessive genetic combinations are shown as "4-Blocks," such as the following:
[B][G] = Blue in the top left corner, Green in the top right corner
[B][B] = Blue in the bottom left corner, Blue in the bottom right corner
The top row of the 4-block, e.g., "B G," is the fish's genetic code for body color, and the bottom row, e.g., "B B," is for fin color. Green is the dominant gene both for body and fin color, whereas blue is the recessive gene for those attributes. Thus, a green-green top row makes for a green body color and so does green-blue and blue-green. Only a blue-blue top row would give a blue body color. The same applies to the bottom row, with respect to fin color. For example, the fish with the genes in the 4-block above, has a green body and a blue fin.
When two fish mate, the offspring's each parent contributes one genetic "square" from each of its rows. The selection is random. It is therefore possible for two green-bodied fishes to beget a blue-bodied fishlet.
When you add fish to the pool, each fish's genotype is selected randomly from the sample space of all 16 different combinations that a 4-block can take. Fishes' phenotype is determined directly by their genotype. Upon activating the simulation, the fish swim around randomly. If at least two fish are on the same square (a NetLogo "patch"), they might mate and procreate if they are matched according to the pre-set mating rules (see 'THINGS TO NOTICE'). For instance, the fish might need to have the same fin color (see 'HOW TO USE IT').
When the fish mate, the offspring's genotype, both for the body color and the fin color, is determined as a combination of the parents' genetic material (combinations of blocks from their respective 4-Blocks). When mating is selective, certain offspring are less likely to appear, so the population distribution changes, and these changes, in turn, further impact mating chances. The simulation enables you to witness these processes. Interface monitors and graphs keep track of fish distribution both by genotype and phenotype.
Press SETUP, and then add as many fish as you'd like, in increments of 10, by pressing ADD FISH. Change the mating rules in the MATE-WITH chooser and the LIFE-SPAN of a fish. Then, press GO ONCE or GO to see the fish interact. Below are more features that will let you take full advantage of the model's capabilities.
LIFE-SPAN -- sets the number of "years," or time-steps, a fish will live.
SETUP -- initializes variables, re-colors the tiles, and resets the monitors and the graphs.
ADD FISH -- adds 10 fish to the display, updating the 4-BLOCK DISTRIBUTION histogram.
GO ONCE -- activates a single run through the 'go' procedure during which the fish move one step in the direction they were facing; if two fish land on the same patch and they are suited to mate according to the current settings of the mating rules, an offspring of the two fish is born and appears somewhere in the world. A patch turns yellow to indicate that two fish are mating on it.
GO -- forever button that keeps running through the 'go' procedure (GO ONCE runs through the procedure only once).
REVEAL GENES -- when the button is pressed, click on specific fish on the display to view their 4-Block Genotype. If you click on a mating fish (it's on a yellow patch), the genes of both of the parents and the child are revealed at the top of the display, and all other fish are temporarily hidden.
G-BODY G-FIN -- displays the number of fish with a green body and a green fin.
G-BODY B-FIN -- displays the number of fish with a green body and a blue fin.
B-BODY G-FIN -- displays the number of fish with a blue body and a green fin.
B-BODY B-FIN -- displays the number of fish with a blue body and a blue fin.
G-BODY TOTAL -- displays the number of fish with a green body (G-BODY G-FIN + G-BODY B-FIN).
B-BODY TOTAL -- displays the number of fish with a blue body (B-BODY G-FIN + B-BODY B-FIN).
TOTAL G-FIN -- displays the number of fish with a green fin (G-BODY G-FIN + B-BODY G-FIN).
TOTAL B- FIN -- displays the number of fish with a blue fin (G-BODY B-FIN + B-BODY B-FIN).
# of ITERATIONS -- displays the number of time-steps the current experiment had run.
AVE BLOCK -- gives the mean number of green squares in the genetic material (4-Blocks) of the entire fish population.
PERCENT FISH BY PROPERTIES -- keeps track of the trends of the different sub-populations over multiple time-steps.
4-BLOCK DISTRIBUTION -- histogram of fish genotype, grouped by the number of green squares in the 4-Blocks. The vertical red line represents the AVE BLOCK.
MATE-WITH -- Represents the four different mating-rule choices:
Press SETUP and then press ADD FISH 10 times. Watch the group of monitors in the top-right area of the interface. What can you say about the distribution of fish by body and fin color? Why is this distribution consistent when you repeat this little activity? Why do we get this specific initial distribution?
What is the setting of the LIFE-SPAN slider that would keep the population from depleting? How, if at all, does this depend on the setting of the MATE-WITH choice button and on the initial size of the population?
If you set the MATE-WITH choice button to "SAME BODY," fish will only mate with other fish that have the same body color and they do (they could differ in fin color). How do you expect this should affect trends in the distribution of fish by phenotype, as viewed in the PERCENT FISH BY PROPERTIES plot?
In the current version of the model, at each 'go' only two fish per patch are selected randomly for a possible mating even if more than two fish are present on that patch. In fact, two non-matching fish may be selected even there are matches on that patch. This constraint limits the maximum reproduction rate to the number of patches in the graphic window. Thus, the population may reach equilibrium, rather than explode. See 'EXTENDING THE MODEL' for ideas on how to make the model more realistic.
What is the total number of fish when the population reaches equilibrium, and how does this number depend on the LIFE-SPAN?
An experiment may reach a point where the number of blue-bodied fish is 0. However, as the experiment continues, one or more blue-bodied fish may appear. Can you explain this? Set MATE-WITH to "Same Body" for an example of such behavior.
Note that the PERCENT FISH BY PROPERTIES plot tracks the fish population by phenotype, while the 4-BLOCK DISTRIBUTION plot represents the fish population by genotype.
Note that for low values of LIFE-SPAN, your fish population may deplete and for higher values, the population might grow sharply but eventually arrive at equilibrium, at a point that is determined by the life span.
Press SETUP and then press ADD FISH 10 times. Watch the group of monitors in the top-right area of the model. What can you say about the distribution of fish by body and by fin color? Why is this distribution quite consistent each time you repeat this little activity? Why do we get this specific initial distribution? What if we added a 1000 fish --would the distribution be as far off of your expectation (in terms of absolute numbers?; in proportionate terms?).
Run the model with MATE-WITH set at "Any Fish." After the model has been running for a while, will the population make-up stay the same, or will it change? Do you expect any type of fish to dominate? Which one and why? Which 4-Block will dominate?
Run GO with MATE-WITH set at "Same Body." After the model has been running for a while, what do you predict will happen to the fish population?
How about MATE-WITH set at "Same Fin"? Which 4-Block would dominate now? After a while, what is the relationship between the G-G and G-B population? Do these outcomes depend on the starting population? If so, how so?
How do you expect the setting of MATE-WITH to affect trends in the distribution of fish by phenotype, as viewed in the "Percent Fish by Properties" graph?
Run GO with MATE-WITH set at "Same Both." What do you think will happen? Is it just a faster way of getting to the same body/same fin final result? Track the 4-Block population while doing so by slowing down the model. Over repeated trials, you may notice that the population is quite sensitive to the random configuration of fish at the starting point? (See 'EXTENDING THE MODEL').
Is there a connection between the mating rules and the tendency of the population to die out?
Change the mating procedures of the model so that instead of only two fish mating per patch, all fish on a patch get an opportunity to mate. For example, if there are 6 fish on the same patch, a total of three pairs may reproduce. This would make the model more realistic.
The model continuously plots four different fish populations in terms of the number in each. You may want to track different aspects of the population, such as according to the number of green squares in each fish's 4-block. Add a plot and edit the 'update-graphs' procedure to do so.
The current interface of the model allows the user to add exactly 10 fish at a time. Give the user more control over the number of fish to be added. To do so, you can add a 'user-input,' asking for the number of fish to be added, or you may want to add a slider where the user can select the number of fish to be added. A simple way of adding an exact number of fish, x, is to type 'add-fish x' in the Command Center.
What is the connection between the original set of fish, the MATE-WITH choice, and the set of fish that the population ultimately "settles" on? Can you determine any rules? To do this, you would first want to create buttons that let the user add specific types of fish, for instance a blue-bodied/blue-finned fish. You could then explore the following question, for example: if the original population set has a higher presence of B-BODY B-FIN fish, could this phenotype become dominant under any setting of MATE-WITH?
The model utilizes the mouse-clicking capabilities of NetLogo. Note that a forever button must be running in order for the model to detect a mouse click.
Note the white band at the top of the view. The procedures ensure that no fish may enter the area of the white band. The origin is moved so that the white band has positive y-coordinates so bounds checking is easier. The band was designed for this model as a unique form of output to the user, because we required output features that cannot be displayed in a standard output area.
There are 16 different 4-blocks, for each of which a new shape had to be created in order for the fish to be able to REVEAL GENES. While drawing one shape is not a very arduous task in NetLogo's Shapes Editor, drawing 16 or more is. (And imagine the task of creating all 512 shapes of a "9-block"!) To create the blocks, we opened the NetLogo file in WordPad, in which all shapes are represented in text as circles, polygons and lines with coordinates and colors. Using this text-based version of the graphical elements greatly simplified the task of managing and generating the entire set of desired 4-Block shapes -- we simply used a straightforward combinatorial analysis expressed in simple symbols. When we re-loaded the file in NetLogo, the entire set of 4-blocks was there just as if we had made them in the Shapes Editor.
ProbLab Genetics, though it may not seem so at first glance, is very closely related to the 9-Block Stalagmite model of the ProbLab curriculum. In that model, 4-blocks are chosen and then stacked in a "histogram," to show the distribution of the blocks by type. Similarly, we may track the distribution of the 4-blocks, the fishes' genotype, upon adding fish and watching the 4-BLOCK DISTRIBUTION histogram update.
In terms of interactions between breeds, several NetLogo models track the number of specimens in breeds that compete over resources. For an example, see the Wolf-Sheep Predation model.
In Expected Value Advanced, a model of the ProbLab Curriculum, the same metaphor of roaming fish is used. As in this model, where a fish carries the value of the number of green squares in its 4-block, fish in Expected Value Advanced carry an underlying monetary value that may be collected and quantified. Additionally, the fish's appearance is dependent on their underlying value (the larger the value, the lighter the body color).
This model is a part of the ProbLab curriculum. For more information about the ProbLab Curriculum please refer to http://ccl.northwestern.edu/curriculum/ProbLab/.
Thanks to Steve Gorodetskiy for his design and programming work.
If you mention this model or the NetLogo software in a publication, we ask that you include the citations below.
For the model itself:
Please cite the NetLogo software as:
Copyright 2004 Uri Wilensky.
This work is licensed under the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License. To view a copy of this license, visit https://creativecommons.org/licenses/by-nc-sa/3.0/ or send a letter to Creative Commons, 559 Nathan Abbott Way, Stanford, California 94305, USA.
Commercial licenses are also available. To inquire about commercial licenses, please contact Uri Wilensky at email@example.com.
This model was created as part of the projects: PARTICIPATORY SIMULATIONS: NETWORK-BASED DESIGN FOR SYSTEMS LEARNING IN CLASSROOMS and/or INTEGRATED SIMULATION AND MODELING ENVIRONMENT. The project gratefully acknowledges the support of the National Science Foundation (REPP & ROLE programs) -- grant numbers REC #9814682 and REC-0126227.