Front page Wednesday, March 28, 2007 Will an infection become an epidemic? A simple model   posted by p-ter @ 3/28/2007 03:17:00 PM One subject I've developed a mild interest in is the modeling of epidemics. So in that vein, here is a post (the first of a possible series, if I feel really inspired) on mathematical epidemiology:Let's start with a simple discrete-time model for the evolution of an infection: imagine an infinite population of people, of which one has become infected with a disease. This one person will infect each of the other susceptible members of the population with some probability p. In the next time unit, the infected individual will recover and become immune to further infection (i.e. he either has developed resistance to the disease or is dead). Each of the newly infected individuals then independently can infect each of the remaining susceptible individuals with probability p once more (so note that the probability of being infected is p if there is one infected individual in the population, 2p if there are two, 3p if there are three, etc), and so on and so forth. The basic model is shown in the figure-- asusceptible individual (S) can be infected, becoming an infected individual (I). Infected individuals then always recover and are removed from the population (R).One of the first (and most mathematically tractable) questions to ask of this model is: will an infection become an epidemic?To answer this question, let's think of the first stages of the infection of the population as a Galton-Watson process (see this wikipedia link if you're wondering what role Francis Galton played in outlining this process). That is, each infected individual infects a random number of susceptible individuals in the next time point according to some probability distribution (called the "offspring distribution") X. If the probability of infecting each individual person is sufficiently low, we can say this distribution is Poisson(lambda). The quantity of interest, then, is the probability that this Galton-Watson process goes extinct; that is, that there exists a time point in which zero people are infected. Let's call this probability E. At a given point in time, if there is one infected individual, the probablity of excinction is E. If there are two infected individuals, the probability of extinction is E^2. Three infected individuals: E^3. And so on. Given that we start with one infected individual, we can than write:E = P(X=0) + P(X=1)E + P(X=2)E^2 + P(X=3)E^3 + P(X=4)E^4 ...Mathematically oriented readers will notice that the right-hand side of this equation is the probability generating function of the offspring distribution, evaluated at E. Indeed, the extinction probability of a Galton-Watson process is the minumum solution to the equation E = A(E) [1]. For our Poisson distribution, this has simple consequences: if the mean of the offspring distribution is less than 1, the infection will not become an epidemic (the only solution to the equation is E = 1). If the mean of the distribution is greater than or equal to one, there is some probability that the Galton-Watson process will never die (in a real population, of course, the epidemic will eventually run out of gas, but we're not interested in those dynamics here). This may seem rather intuitive (if on average, each infected individual infected less than one individual, the infection will die out quickly), but the mathematical formalism allows one to construct more complex models. Plus, reducing lambda (the mean of the Poission distribution), of course, is then the goal of nearly all anti-epidemic precautions-- quarantines, hand-washing, etc.; all can be though of as practical ways to reduce a mathematical parameter. [1] The probability generation function A(z) of a distribution X is the sum over all possible values of X of P(X=x)z^x. Labels: Epidemiology input{width:180px; border: thin #000000 inset;} .sub{width:180px; border: thin #000000 inset;} Razib's Home Page GNXP Archives Interviews Blogroll Principles of Population Genetics Genetics of Populations Molecular Evolution Quantitative Genetics Evolutionary Quantitative Genetics Evolutionary Genetics Evolution Molecular Markers, Natural History, and Evolution The Genetics of Human Populations Genetics and Analysis of Quantitative Traits Epistasis and Evolutionary Process Evolutionary Human Genetics Biometry Mathematical Models in Biology Speciation Evolutionary Genetics: Case Studies and Concepts Narrow Roads of Gene Land 1 Narrow Roads of Gene Land 2 Narrow Roads of Gene Land 3 Statistical Methods in Molecular Evolution The History and Geography of Human Genes Population Genetics and Microevolutionary Theory Population Genetics, Molecular Evolution, and the Neutral Theory Genetical Theory of Natural Selection Evolution and the Genetics of Populations Genetics and Origins of Species Tempo and Mode in Evolution Causes of Evolution Evolution The Great Human Diasporas Bones, Stones and Molecules Natural Selection and Social Theory Journey of Man Mapping Human History The Seven Daughters of Eve Evolution for Everyone Why Sex Matters Mother Nature Grooming, Gossip, and the Evolution of Language Genome R.A. Fisher, the Life of a Scientist Sewall Wright and Evolutionary Biology Origins of Theoretical Population Genetics A Reason for Everything The Ancestor's Tale Dragon Bone Hill Endless Forms Most Beautiful The Selfish Gene Adaptation and Natural Selection Nature via Nurture The Symbolic Species The Imitation Factor The Red Queen Out of Thin Air Mutants Evolutionary Dynamics The Origin of Species The Descent of Man Age of Abundance The Darwin Wars The Evolutionists The Creationists Of Moths and Men The Language Instinct How We Decide Predictably Irrational The Black Swan Fooled By Randomness Descartes' Baby Religion Explained In Gods We Trust Darwin's Cathedral A Theory of Religion The Meme Machine Synaptic Self The Mating Mind A Separate Creation The Number Sense The 10,000 Year Explosion The Math Gene Explaining Culture Origin and Evolution of Cultures Dawn of Human Culture The Origins of Virtue Prehistory of the Mind The Nurture Assumption The Moral Animal Born That Way No Two Alike Sociobiology Survival of the Prettiest The Blank Slate The g Factor The Origin Of The Mind Unto Others Defenders of the Truth The Cultural Origins of Human Cognition Before the Dawn Behavioral Genetics in the Postgenomic Era The Essential Difference Geography of Thought The Classical World The Fall of the Roman Empire The Fall of Rome History of Rome How Rome Fell The Making of a Christian Aristoracy The Rise of Western Christendom Keepers of the Keys of Heaven A History of the Byzantine State and Society Europe After Rome The Germanization of Early Medieval Christianity The Barbarian Conversion A History of Christianity God's War Infidels Fourth Crusade and the Sack of Constantinople The Sacred Chain Divided by the Faith Europe The Reformation Pursuit of Glory Albion's Seed 1848 Postwar From Plato to Nato China: A New History China in World History Genghis Khan and the Making of the Modern World Children of the Revolution When Baghdad Ruled the Muslim World The Great Arab Conquests After Tamerlane A History of Iran The Horse, the Wheel, and Language A World History Guns, Germs, and Steel The Human Web Plagues and Peoples 1491 A Concise Economic History of the World Power and Plenty A Splendid Exchange Contours of the World Economy 1-2030 AD Knowledge and the Wealth of Nations A Farewell to Alms The Ascent of Money The Great Divergence Clash of Extremes War and Peace and War Historical Dynamics The Age of Lincoln The Great Upheaval What Hath God Wrought Freedom Just Around the Corner Throes of Democracy Grand New Party A Beautiful Math When Genius Failed Catholicism and Freedom American Judaism Search this site: (2005 - present, Blogger Archives) Hello Movable Type archives August 11,2002 August 18,2002 August 25,2002 September 01,2002 September 15,2002 October 20,2002 December 08,2002 December 22,2002 December 29,2002 January 05,2003 January 12,2003 January 19,2003 January 26,2003 February 02,2003 February 09,2003 February 16,2003 February 23,2003 March 02,2003 March 09,2003 March 16,2003 March 23,2003 March 30,2003 April 06,2003 April 13,2003 April 20,2003 April 27,2003 May 04,2003 May 11,2003 May 18,2003 May 25,2003 June 01,2003 June 08,2003 June 15,2003 June 22,2003 June 29,2003 July 06,2003 July 13,2003 July 20,2003 July 27,2003 August 03,2003 August 10,2003 August 17,2003 August 24,2003 August 31,2003 September 07,2003 September 14,2003 September 21,2003 September 28,2003 October 05,2003 October 12,2003 October 19,2003 October 26,2003 November 02,2003 November 09,2003 November 16,2003 November 23,2003 November 30,2003 December 07,2003 December 14,2003 December 21,2003 December 28,2003 January 04,2004 January 11,2004 January 18,2004 January 25,2004 February 01,2004 February 08,2004 February 15,2004 February 22,2004 February 29,2004 March 07,2004 March 14,2004 March 21,2004 March 28,2004 April 04,2004 April 11,2004 April 18,2004 April 25,2004 May 02,2004 May 09,2004 May 16,2004 May 23,2004 May 30,2004 June 06,2004 June 13,2004 June 20,2004 June 27,2004 July 04,2004 July 11,2004 July 18,2004 July 25,2004 August 01,2004 August 08,2004 August 15,2004 August 22,2004 August 29,2004 September 05,2004 September 12,2004 September 19,2004 September 26,2004 October 03,2004 October 10,2004 October 17,2004 October 24,2004 October 31,2004 November 07,2004 November 14,2004 November 21,2004 November 28,2004 December 05,2004 December 12,2004 December 19,2004 December 26,2004 January 02,2005 January 09,2005 January 16,2005 January 23,2005 January 30,2005 February 06,2005 February 13,2005 February 20,2005 February 27,2005 March 06,2005 March 13,2005 March 20,2005 March 27,2005 April 03,2005 April 10,2005 April 17,2005 April 24,2005 May 01,2005 May 08,2005 May 15,2005 May 22,2005 May 29,2005 June 05,2005 June 12,2005 June 19,2005 June 26,2005 July 03,2005 July 17,2005 August 07,2005 Blogspot archives June 2002 July 2002 August 2002 September 2002 October 2002 November 2002 December 2002