2016-02-20: Apple VS FBI

2016-02-19: More Zika may be better than less

2016-02-17: Dependent Non-Commuting Random Variable Systems

2016-01-14: Life at the multifurcation

2015-09-28: AI ain't that smart

2015-06-24: MathEpi citation tree

2015-03-31: Too much STEM is bad

2015-03-24: Dawn of the CRISPR age

2015-02-12: A Comment on How Biased Dispersal can Preclude Competitive Exclusion

2015-02-09: Hamilton's selfish-herd paradox

2015-02-08: Risks and values of microparasite research

2014-11-10: Vaccine mandates and bioethics

2014-10-18: Ebola, travel, president

2014-10-17: Ebola comments

2014-10-12: Ebola numbers

2014-09-23: More stochastic than?

2014-08-17: Feynman's missing method for third-orders?

2014-07-31: CIA spies even on congress

2014-07-16: Rehm on vaccines

2014-06-21: Kurtosis, 4th order diffusion, and wave speed

2014-06-20: Random dispersal speeds invasions

2014-05-06: Preservation of information asymetry in Academia

2014-04-16: Dual numbers are really just calculus infinitessimals

2014-04-14: More on fairer markets

2014-03-18: It's a mad mad mad mad prisoner's dilemma

2014-03-05: Integration techniques: Fourier--Laplace Commutation

2014-02-25: Fiber-bundles for root-polishing in two dimensions

2014-02-17: Is life a simulation or a dream?

2014-01-30: PSU should be infosocialist

2014-01-12: The dark house of math

2014-01-11: Inconsistencies hinder pylab adoption

2013-12-24: Cuvier and the birth of extinction

2013-12-17: Risk Resonance

2013-12-15: The cult of the Levy flight

2013-12-09: 2013 Flu Shots at PSU

2013-12-02: Amazon sucker-punches 60 minutes

2013-11-26: Zombies are REAL, Dr. Tyson!

2013-11-22: Crying wolf over synthetic biology?

2013-11-21: Tilting Drake's Equation

2013-11-18: Why $1^\infty != 1$

2013-11-15: Adobe leaks of PSU data + NSA success accounting

2013-11-14: 60 Minutes misreport on Benghazi

2013-11-11: Making fairer trading markets

2013-11-10: L'Hopital's Rule for Multidimensional Systems

2013-11-09: Using infinitessimals in vector calculus

2013-11-08: Functional Calculus

2013-11-03: Elementary mathematical theory of the health poverty trap

2013-11-02: Proof of the area of a circle using elementary methods

Hamilton's selfish-herd paradox

Herd formation is one of the most obvious examples of emergent spatial structure in animal populations. A number of recent theory papers have addressed how behavior contributes to the emergence of animal herds. A number of methods have been used in these papers. The simplest approach is to assign movement rules to individuals without internal states and then compare the emergent spatial structures to observations of real aggregations. For example, Lee 2006 uses molecular dynamics to create snapshots of aggregations comparable to herds. More complex models may assign movement rules that depend on both internal and external states. Rands 2004 use a rule-of-thumb agent model to study emergent patterns of social foragers under predation. In other cases, models with assigned behaviors can be extended to allow behaviors to evolve. To test how selection affects aggregation behaviors, Reluga 2005 uses an agent-based simulation with a simple evolution rule to show that realistic aggregations can be stable under predation pressure. Wood 2007 combined foraging and predation and obtained more evidence that trade-offs between foraging benefits and predation risk can drive the evolution of herd behaviors.

Surprisingly, the game-theoretic foundations of aggregation behavior have not received much attention (perhaps because Hamilton's informal arguments were convincing). While thought experiments and simulation experiments have shown predation can lead to the evolution of herds, few if any mathematically rigorous results seem to have been obtained so far. Rigorous game-theoretic results would be a great addition to selfish-herd theory.

One often-invoked explanation is the Hamilton's selfish-herd hypothesis. Hamilton originally formulated his hypothesis in terms of frogs living on the circular edge of a pond. Suppose we have \(n\) frogs, and each individual \(i\) has a position \(x_i(t) \in [0,2 \pi)\) on the pond's edge at time \(t\) (for convenience, we use polar coordinates). At any time \(t\), a frog may choose to move around the circle in one direction or the other. The velocity chose by individual \(i\) at time \(t\) is \(\phi_i(t) \in [-1,1]\). The dynamics are governed by \[ \dot{x}_i = \phi_i(t), \quad \forall i \in 1 \ldots n \] At some time \(t_f > 0\), a snake emerges from the water and eats one of the frogs near to where it emerges. Thus, there a risk of death that depends on the positions of the frogs. If each frog wishes to minimize it's risk of death, how should the frogs behave?

This is clearly a game-theory problem. Given initial conditions, can we calculate a Nash equilibrium strategy-set for the movement of the frogs? To determine the Nash equilibria, we first have to specify payoffs for individual behaviors in terms of a risk function \(R_i\) for individual \(i\) that depends on the positions of other individuals at the final time \(t_f\). In the standard thought experiment, the risk of death \(R_i\) is the probability that the predator's position pops up closer to individual \(i\) than to any other individual. If individual \(i\) is between individual \(i-1\) and individual \(i+1\), then \[ R_i = \frac{|x_{i+1}(t_f) - x_{i-1}(t_f)|}{4 \pi}. \] This risk function leads to a singular model in the sense that it is independent of individual \(i\)'s position. Individual \(i\) can not change it's risk through movement. This is a subtle but important issue in the mathematical analysis. Small perturbations to the shape of \(R_i\) can lead to very different equilibria.

Thus, we have reached a bit of a paradox. Clearly, herd formation is evolutionarily stable, given the frequency with which we find it in nature. However, in theory, Hamilton's selfish-herd hypothesis -- the most frequently invoked explanation, says individuals can not change their risk of predation by movement!

A couple of caveates are in order at this point. For one, Hamilton used a discrete-time model with sequential updating and unbounded jumps in his original explanation of selfish herds, which avoids the problem we have revealed. However, in this sense, our model is more realistic than his -- animals can not teleport, but rather must move continuously through space from one point to another! A second related caveate is that in most cases, herd behavior is occurring in 2 or 3 dimensions, rather than the 1 dimension of Hamilton's thought experiment. In 2 and 3 dimensions, the size of the domain of danger DOES change according to the individual's position relative to it's neighbors. However, several simple calculations I have done suggest that moving towards a neighbor does NOT reduce one's domain of danger.

It seems that something else is going on, and some further mathematical studies are needed to justify the selfish herd hypothesis from first principles.