Kinetic theory of age-structured stochastic birth-death processes

Classical age-structured mass-action models such as the McKendrick-von Foerster equation have been extensively studied but are unable to describe stochastic fluctuations or population-size-dependent birth and death rates. Stochastic theories that treat semi-Markov age-dependent processes using, e.g., the Bellman-Harris equation do not resolve a population's age structure and are unable to quantify population-size dependencies. Conversely, current theories that include size-dependent population dynamics (e.g., mathematical models that include carrying capacity such as the logistic equation) cannot be easily extended to take into account age-dependent birth and death rates. In this paper, we present a systematic derivation of a new, fully stochastic kinetic theory for interacting age-structured populations. By defining multiparticle probability density functions, we derive a hierarchy of kinetic equations for the stochastic evolution of an aging population undergoing birth and death. We show that the fully stochastic age-dependent birth-death process precludes factorization of the corresponding probability densities, which then must be solved by using a Bogoliubov-–Born–-Green–-Kirkwood-–Yvon-like hierarchy. Explicit solutions are derived in three limits: no birth, no death, and steady state. These are then compared with their corresponding mean-field results. Our results generalize both deterministic models and existing master equation approaches by providing an intuitive and efficient way to simultaneously model age- and population-dependent stochastic dynamics applicable to the study of demography, stem cell dynamics, and disease evolution

Inferring transient dynamics of human populations from matrix non-normality

This is the final version of the article. Available from Springer Verlag via the DOI in this record.In our increasingly unstable and unpredictable world, population dynamics rarely settle uniformly to long-term behaviour. However, projecting period-by-period through the preceding fluctuations is more data-intensive and analytically involved than evaluating at equilibrium. To efficiently model populations and best inform policy, we require pragmatic suggestions as to when it is necessary to incorporate short-term transient dynamics and their effect on eventual projected population size. To estimate this need for matrix population modelling, we adopt a linear algebraic quantity known as non-normality. Matrix non-normality is distinct from normality in the Gaussian sense, and indicates the amplificatory potential of the population projection matrix given a particular population vector. In this paper, we compare and contrast three well-regarded metrics of non-normality, which were calculated for over 1000 age-structured human population projection matrices from 42 European countries in the period 1960 to 2014. Non-normality increased over time, mirroring the indices of transient dynamics that peaked around the millennium. By standardising the matrices to focus on transient dynamics and not changes in the asymptotic growth rate, we show that the damping ratio is an uninformative predictor of whether a population is prone to transient booms or busts in its size. These analyses suggest that population ecology approaches to inferring transient dynamics have too often relied on suboptimal analytical tools focussed on an initial population vector rather than the capacity of the life cycle to amplify or dampen transient fluctuations. Finally, we introduce the engineering technique of pseudospectra analysis to population ecology, which, like matrix non-normality, provides a more complete description of the transient fluctuations than the damping ratio. Pseudospectra analysis could further support non-normality assessment to enable a greater understanding of when we might expect transient phases to impact eventual population dynamics.This work was funded by Wellcome Trust New Investigator 103780 to TE, who is also funded by NERC Fellowship NE/J018163/1. JB gratefully acknowledges the ESRC Centre for Population Change ES/K007394/1

Monitoring international migration flows in Europe. Towards a statistical data base combining data from different sources

The paper reviews techniques developed in demography, geography and statistics that are useful for bridging the gap between available data on international migration flows and the information required for policy making and research. The basic idea of the paper is as follows: to establish a coherent and consistent data base that contains sufficiently detailed, up-to-date and accurate information, data from several sources should be combined. That raises issues of definition and measurement, and of how to combine data from different origins properly. The issues may be tackled more easily if the statistics that are being compiled are viewed as different outcomes or manifestations of underlying stochastic processes governing migration. The link between the processes and their outcomes is described by models, the parameters of which must be estimated from the available data. That may be done within the context of socio-demographic accounting. The paper discusses the experience of the U.S. Bureau of the Census in combining migration data from several sources. It also summarizes the many efforts in Europe to establish a coherent and consistent data base on international migration. The paper was written at IIASA. It is part of the Migration Estimation Study, which is a collaborative IIASA-University of Groningen project, funded by the Netherlands Organization for Scientific Research (NWO). The project aims at developing techniques to obtain improved estimates of international migration flows by country of origin and country of destination

All grown up? The fate after 15 years of a quarter of a million UK firms born in 1998

The theory of firm growth is in a rather unsatisfactory state. However, the analysis of large firm-level datasets which have become available in recent years allows us to begin building an evidence base which can, in turn, be used to underpin the development of more satisfactory theory. Here we study the 239 thousand UK private sector firms born in 1998 over their first 15 years of life. A first, and quite striking, finding is the extraordinary force of mortality. By age 15, 90% of the UK firms born in 1998 are dead, and, for those surviving to age 15, the hazard of death is still about 10% a year. The chance of death is related to the size and growth of firms in an interesting way. Whilst the hazard rate after 15 years is largely independent of size at birth, it is strongly affected by the current (age 14) size. In particular, firms with more than five employees are half as likely to die in the next year as firms with less than five employees. A second important finding is that most firms, even those which survive to age 15, do not grow very much. By age 15 more than half the 26,000 survivors still have less than five jobs. In other words, the growth paths – what we call the ‘growth trajectories’ – of most of the 26,000 survivors are pretty flat. However, of the firms that do grow, firms born smaller grow faster than those born larger. Another striking finding is that growth is heavily concentrated in the first five years. Whilst growth does continue, even up to age 15, each year after age five it involves only a relatively small proportion of firms. Finally, there are two groups of survivors which contribute importantly to job creation. Some are those born relatively large (with more than 20 jobs) although their growth rate is quite modest. More striking though, is a very small group of firms born very small with less than five jobs (about 5% of all survivors) which contribute a substantial proportion (more than one third) of the jobs added to the cohort total by age 15

The Leverage of Demographic Dynamics on Carbon Dioxide Emissions: Does Age Structure Matter?

This article provides a methodological contribution to the study of the effect of changes in population age structure on carbon dioxide (CO2) emissions. First, I propose a generalization of the IPAT equation to a multisector economy with an age-structured population and discuss the insights that can be obtained in the context of stable population theory. Second, I suggest a statistical model of household consumption as a function of household size and age structure to quantitatively evaluate the extent of economies of scale in consumption of energy-intensive goods, and to estimate age-specific profiles of consumption of energy-intensive goods and of CO2 emissions. Third, I offer an illustration of the methodologies using data for the United States. The analysis shows that per-capita CO2 emissions increase with age until the individual is in his or her 60s, and then emissions tend to decrease. Holding everything else constant, the expected change in U.S. population age distribution during the next four decades is likely to have a small, but noticeable, positive impact on CO2 emissions

A Multigenerational View of Inequality

The study of intergenerational mobility and most population research are governed by a two-generation (parent-to-offspring) view of intergenerational influence, to the neglect of the effects of grandparents and other ancestors and nonresident contemporary kin. While appropriate for some populations in some periods, this perspective may omit important sources of intergenerational continuity of family-based social inequality. Social institutions, which transcend individual lives, help support multigenerational influence, particularly at the extreme top and bottom of the social hierarchy, but to some extent in the middle as well. Multigenerational influence also works through demographic processes because families influence subsequent generations through differential fertility and survival, migration, and marriage patterns, as well as through direct transmission of socioeconomic rewards, statuses, and positions. Future research should attend more closely to multigenerational effects; to the tandem nature of demographic and socioeconomic reproduction; and to data, measures, and models that transcend coresident nuclear families

Real-Time Epidemic Monitoring and Forecasting of H1N1-2009 Using Influenza-Like Illness from General Practice and Family Doctor Clinics in Singapore

10.1371/journal.pone.0010036PLoS ONE54

An exactly solvable, spatial model of mutation accumulation in cancer

One of the hallmarks of cancer is the accumulation of driver mutations which increase the net reproductive rate of cancer cells and allow them to spread. This process has been studied in mathematical models of well mixed populations, and in computer simulations of three-dimensional spatial models. But the computational complexity of these more realistic, spatial models makes it difficult to simulate realistically large and clinically detectable solid tumours. Here we describe an exactly solvable mathematical model of a tumour featuring replication, mutation and local migration of cancer cells. The model predicts a quasi-exponential growth of large tumours, even if different fragments of the tumour grow sub-exponentially due to nutrient and space limitations. The model reproduces clinically observed tumour growth times using biologically plausible rates for cell birth, death, and migration rates. We also show that the expected number of accumulated driver mutations increases exponentially in time if the average fitness gain per driver is constant, and that it reaches a plateau if the gains decrease over time. We discuss the realism of the underlying assumptions and possible extensions of the model