Demographic Modelling of Human Population Growth

This work presents models and methods that have been used in producing forecasts of population growth. The work is intended to emphasize the reliability bounds of the model forecasts. Leslie model and various versions of logistic population models are presented. References to literature and several studies are given. A lot of relevant methodology has been developed in biological sciences. The Leslie modelling approach involves the use of current trends in mortality,fertility, migration and emigration. The model treats population divided in age groups and the model is given as a recursive system. Other group of models is based on straightforward extrapolation of census data. Trajectories of simple exponential growth function and logistic models are used to produce the forecast. The work presents the basics of Leslie type modelling and the logistic models, including multi- parameter logistic functions. The latter model is also analysed from model reliability point of view. Bayesian approach and MCMC method are used to create error bounds of the model predictions.

Étude de l’évolution du tropisme et de la pression sélective exercée sur le gène de l’enveloppe du virus de l’immunodéficience humaine de type 1 au cours de la grossesse et à l’échelle de la population.

Introduction. Le VIH-1 évolue en fonction de la réponse immunitaire spécifique de l’hôte. La pression sélective exercée par la réponse immunitaire VIH-spécifique de l’hôte entraine l’évolution des gènes viraux et à terme détermine l’évolution de la maladie. Cette évolution du virus à l’échelle d’un individu façonne également l’évolution du virus à l’échelle de la population et détermine le devenir de l’épidémie. Le VIH utilise les corécepteurs d’entrée CCR5 (virus R5) et CXCR4 (virus X4) afin d’infecter la cellule cible, et l’évolution du tropisme du virus de R5 vers X4, appelé switch du tropisme, est associé à la progression de la maladie. Les virus R5 sont rencontrés en début d’infection tandis que les virus X4 apparaissent en fin de maladie chez un certain de nombre de patients et sont considérés comme plus virulents. La pression sélective immunitaire exercée sur le gène de l’enveloppe (env) peut donc entrainer l’évolution du tropisme du VIH. La grossesse est un état immunitaire particulier considéré comme étant principalement caractérisé par un biais Th2 nécessaire à l’établissement de la tolérance materno-fétale. Le switch de tropisme de R5 vers X4 en grossesse n’a jamais été documenté, de même que l’évolution des déterminants du tropisme à l’échelle de la population. Hypothèses. Les changements immunitaires associés à l’initiation et la progression de la grossesse engendrent des changements dans la pression immunitaire exercée sur l’enveloppe et peuvent favoriser le switch du tropisme. L’évolution du tropisme du VIH-1 peut être observé à l’échelle de la population au même titre que l’évolution de l’enveloppe virale. Objectifs. Analyser l’évolution du tropisme et décrire la pression sélective sur l’enveloppe des femmes enceintes infectées par le VIH-1. Analyser l’évolution des déterminants du tropisme à l’échelle de la population. Méthodes. Nous avons dans un premier temps analysé l’évolution des déterminants du tropisme et déterminé le génotype et phénotype du VIH-1 chez 19 femmes enceintes issues de la cohorte du centre maternel et infantile sur le SIDA de l’hôpital Sainte-Justine (CMIS). Nous avons ensuite caractérisé et comparé la pression sélective exercée sur env, par une méthode bayésienne, chez 31 femmes enceinte et 29 femmes non-enceintes. Enfin, nous avons analysé et comparé des déterminants du tropisme entre des séquences d’enveloppe contemporaines et anciennes, issues des bases de données du NCBI. Résultats. Nos résultats montrent la présence de virus X4 chez la moitié de notre cohorte, et un switch de tropisme de R5 vers X4 chez 5/19 sujets. Les séquences des femmes enceintes présentaient des taux de substitutions plus élevées que celles des femmes non-enceintes. La pression sélective dans la région C2 était plus faible chez les femmes enceintes que chez les femmes non-enceintes, et différait dans 4 positions entre ces 2 groupes. Cette sélection diminuait au cours de la grossesse chez les patientes traitées. Enfin, une accumulation de mutations X4 a été observée dans les séquences R5 contemporaines par rapport aux séquences R5 anciennes. Conclusion. Les changements immunitaires associés à la grossesse semblent induire des modifications subtiles dans la pression sélective exercée sur env, suffisant à influencer l’évolution du tropisme de R5 vers X4. Un switch du tropisme à l’échelle de la population impliquerait une épidémie évoluant vers une plus grande virulence du virus. Nos résultats sont d’importance en ce qui concerne la prophylaxie antirétrovirale pour la santé de la mère et la prévention de la transmission mère-enfant du VIH-1. Ils sont aussi importants concernant l’avenir de la thérapie antirétrovirale dans le contexte d’une épidémie évoluant vers une plus grande virulence

Statistical inference of OH concentrations and air mass dilution rates from successive observations of non-methane hydrocarbons in single air masses

Bayesian inference has been used to determine rigorous estimates of hydroxyl radical concentrations () and air mass dilution rates (K) averaged following air masses between linked observations of nonmethane hydrocarbons (NMHCs) spanning the North Atlantic during the Intercontinental Transport and Chemical Transformation (ITCT)-Lagrangian-2K4 experiment. The Bayesian technique obtains a refined (posterior) distribution of a parameter given data related to the parameter through a model and prior beliefs about the parameter distribution. Here, the model describes hydrocarbon loss through OH reaction and mixing with a background concentration at rate K. The Lagrangian experiment provides direct observations of hydrocarbons at two time points, removing assumptions regarding composition or sources upstream of a single observation. The estimates are sharpened by using many hydrocarbons with different reactivities and accounting for their variability and measurement uncertainty. A novel technique is used to construct prior background distributions of many species, described by variation of a single parameter . This exploits the high correlation of species, related by the first principal component of many NMHC samples. The Bayesian method obtains posterior estimates of , K and following each air mass. Median values are typically between 0.5 and 2.0 × 106 molecules cm−3, but are elevated to between 2.5 and 3.5 × 106 molecules cm−3, in low-level pollution. A comparison of estimates from absolute NMHC concentrations and NMHC ratios assuming zero background (the “photochemical clock” method) shows similar distributions but reveals systematic high bias in the estimates from ratios. Estimates of K are ∼0.1 day−1 but show more sensitivity to the prior distribution assumed.

Brief and precarious lives: Infant mortality in contrasting sites from medieval and post-medieval England (AD 850-1859)

This study compares the infant mortality profiles of 128 infants from two urban and two rural cemetery sites in medieval England. The aim of this paper is to assess the impact of urbanization and industrialization in terms of endogenous or exogenous causes of death. In order to undertake this analysis, two different methods of estimating gestational age from long bone lengths were used: a traditional regression method and a Bayesian method. The regression method tended to produce more marked peaks at 38 weeks, while the Bayesian method produced a broader range of ages and were more comparable with the expected "natural" mortality profiles. At all the sites, neonatal mortality (28-40 weeks) outweighed post-neonatal mortality (41-48 weeks) with rural Raunds Furnells in Northamptonshire, showing the highest number of neonatal deaths and post-medieval Spitalfields, London, showing a greater proportion of deaths due to exogenous or environmental factors. Of the four sites under study, Wharram Percy in Yorkshire showed the most convincing "natural" infant mortality profile, suggesting the inclusion of all births (i.e., stillbirths and unbaptised infants).

Sampling uncertainty properties of cloud fraction estimates from random transect observations

[1] Cloud cover is conventionally estimated from satellite images as the observed fraction of cloudy pixels. Active instruments such as radar and Lidar observe in narrow transects that sample only a small percentage of the area over which the cloud fraction is estimated. As a consequence, the fraction estimate has an associated sampling uncertainty, which usually remains unspecified. This paper extends a Bayesian method of cloud fraction estimation, which also provides an analytical estimate of the sampling error. This method is applied to test the sensitivity of this error to sampling characteristics, such as the number of observed transects and the variability of the underlying cloud field. The dependence of the uncertainty on these characteristics is investigated using synthetic data simulated to have properties closely resembling observations of the spaceborne Lidar NASA-LITE mission. Results suggest that the variance of the cloud fraction is greatest for medium cloud cover and least when conditions are mostly cloudy or clear. However, there is a bias in the estimation, which is greatest around 25% and 75% cloud cover. The sampling uncertainty is also affected by the mean lengths of clouds and of clear intervals; shorter lengths decrease uncertainty, primarily because there are more cloud observations in a transect of a given length. Uncertainty also falls with increasing number of transects. Therefore a sampling strategy aimed at minimizing the uncertainty in transect derived cloud fraction will have to take into account both the cloud and clear sky length distributions as well as the cloud fraction of the observed field. These conclusions have implications for the design of future satellite missions. This paper describes the first integrated methodology for the analytical assessment of sampling uncertainty in cloud fraction observations from forthcoming spaceborne radar and Lidar missions such as NASA's Calipso and CloudSat.

Identifying adaptive genetic divergence among populations from genome scans

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of F-ST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of F-ST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of F-ST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

The demand for a healthy diet: estimating the almost ideal demand system with infrequency of purchase

A Bayesian method of estimating multivariate sample selection models is introduced and applied to the estimation of a demand system for food in the UK to account for censoring arising from infrequency of purchase. We show how it is possible to impose identifying restrictions on the sample selection equations and that, unlike a maximum likelihood framework, the imposition of adding up at both latent and observed levels is straightforward. Our results emphasise the role played by low incomes and socio-economic circumstances in leading to poor diets and also indicate that the presence of children in a household has a negative impact on dietary quality.

Sea surface temperature estimation from the Geostationary Operational Environmental Satellite 12 (GOES-12)

This paper describes the techniques used to obtain sea surface temperature (SST) retrievals from the Geostationary Operational Environmental Satellite 12 (GOES-12) at the National Oceanic and Atmospheric Administration’s Office of Satellite Data Processing and Distribution. Previous SST retrieval techniques relying on channels at 11 and 12 μm are not applicable because GOES-12 lacks the latter channel. Cloud detection is performed using a Bayesian method exploiting fast-forward modeling of prior clear-sky radiances using numerical weather predictions. The basic retrieval algorithm used at nighttime is based on a linear combination of brightness temperatures at 3.9 and 11 μm. In comparison with traditional split window SSTs (using 11- and 12-μm channels), simulations show that this combination has maximum scatter when observing drier colder scenes, with a comparable overall performance. For daytime retrieval, the same algorithm is applied after estimating and removing the contribution to brightness temperature in the 3.9-μm channel from solar irradiance. The correction is based on radiative transfer simulations and comprises a parameterization for atmospheric scattering and a calculation of ocean surface reflected radiance. Potential use of the 13-μm channel for SST is shown in a simulation study: in conjunction with the 3.9-μm channel, it can reduce the retrieval error by 30%. Some validation results are shown while a companion paper by Maturi et al. shows a detailed analysis of the validation results for the operational algorithms described in this present article.

Fatores genéticos e ambientais sobre a probabilidade de prenhez precoce em bovinos Caracu

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Efeitos genéticos e ambientais sobre o perímetro escrotal de animais da raça Caracu

Neste trabalho foram avaliados os efeitos de idade (IDS) e peso aos 15 meses (P15) sobre medidas de perímetro escrotal ao sobreano de 1.892 machos da raça Caracu e estimada a herdabilidade desta característica. Utilizaram-se dois modelos de análise: em um dos modelos, foram incluídos o efeito de grupo de contemporâneos (GC) e, como covariável, a idade ao sobreano (efeitos linear e quadrático); e, no outro, o efeito de GC e, como covariável, o peso ao sobreano (efeitos linear e quadrático). A idade não teve efeito significativo no perímetro escrotal, o que está relacionado à homogeneidade dos animais nos grupos de contemporâneos, enquanto os efeitos linear e quadrático de P15 foram significativos sobre o perímetro escrotal, indicando grande influência do peso sobre a variação desta característica. A herdabilidade do perímetro escrotal foi estimada pelo método Bayesiano utilizando-se um modelo animal. O modelo incluiu os efeitos de GC e P15 (linear e quadrático) e os efeitos genéticos aditivo direto e residual. O valor médio estimado de herdabilidade do PE ao sobreano foi 0,38 e comprova que esta característica pode ser utilizada como critério de seleção para precocidade sexual em programas de melhoramento genético de animais da raça Caracu.

Herdabilidades e correlações genéticas de características de machos e fêmeas, em um rebanho bovino da raça Canchim

O objetivo deste trabalho foi estimar as herdabilidades e as correlações genéticas do peso (P12) e do perímetro escrotal (PE12) de machos aos 12 meses de idade, da idade de descarte (TPR, tempo de permanência no rebanho) de fêmeas e do número (ND10) e de quilogramas (QD10) de bezerros desmamados pelas fêmeas em até dez anos de idade, em um rebanho da raça Canchim. Foram utilizadas 1.370, 826, 826, 2.726 e 1.051 observações de TPR, ND10, QD10, P12 e PE12, respectivamente. As estimativas dos componentes de (co)variância foram obtidas pelo método bayesiano, para todas as características em questão, P12, PE12, TPR, ND10 e QD10. O modelo estatístico incluiu, além dos efeitos aleatórios genético aditivo direto e residual, os efeitos fixos de ano de nascimento do animal para todas as características, de mês de nascimento para P12 e PE12 e da covariável idade do animal para PE12. As estimativas de herdabilidade, obtidas pelas análises unicaráter foram iguais a 0,38; 0,52; 0,24; 0,33 e 0,34 para P12, PE12, TPR, ND10 e QD10, respectivamente, indicando que as características possuem variação genética aditiva suficiente para apresentar boa resposta à seleção. As correlações genéticas de TPR (0,33 e 0,33, respectivamente), ND10 (0,38 e 0,30, respectivamente) e QD10 (0,61 e 0,41, respectivamente) com P12 e PE12, obtidas pelas análises bicaráter, sugerem que a seleção com base no peso e no perímetro escrotal dos machos não deve resultar em decréscimo no tempo de permanência das fêmeas no rebanho e no número e quilogramas de bezerros produzidos pelas fêmeas em até dez anos de idade.

Estudo da interação genótipo × ambiente sobre características de crescimento de bovinos de corte utilizando-se inferência bayesiana

The objective of this study was to evaluate the genotype x environment interaction for weaning and yearling weights, daily weight gain from weaning to 12 months of age and the growth performance in Canchim (5/8 Charolais + 3/8 Zebu) beef cattle estimated by a principal components analysis including those three traits. The environment was defined by season of birth (first and second semesters of the year). Genetic parameters were estimated by bayesian method with the Gibbs sampler using bivariate analyses (considering the trait in each of the two seasons as a different one) and models that included the fixed effects of year and month of birth, sex and age of cow (linear and quadratic) and the random effects of animal and residual. The results suggested that genetic evaluation and selection in Canchim beef cattle for the traits studied should consider the genotype and season of birth interaction.

Planning accelerated life tests under Exponentiated-Weibull-Arrhenius model

Purpose - The purpose of this paper is to present designs for an accelerated life test (ALT). Design/methodology/approach - Bayesian methods and simulation Monte Carlo Markov Chain (MCMC) methods were used. Findings - In the paper a Bayesian method based on MCMC for ALT under EW distribution (for life time) and Arrhenius models (relating the stress variable and parameters) was proposed. The paper can conclude that it is a reasonable alternative to the classical statistical methods since the implementation of the proposed method is simple, not requiring advanced computational understanding and inferences on the parameters can be made easily. By the predictive density of a future observation, a procedure was developed to plan ALT and also to verify if the conformance fraction of the manufactured process reaches some desired level of quality. This procedure is useful for statistical process control in many industrial applications. Research limitations/implications - The results may be applied in a semiconductor manufacturer. Originality/value - The Exponentiated-Weibull-Arrhenius model has never before been used to plan an ALT. © Emerald Group Publishing Limited.

Molecular phylogenetic relationships and phenotypic diversity in miniaturized toadlets, genus Brachycephalus (Amphibia: Anura: Brachycephalidae)

Toadlets of the genus Brachycephalus are endemic to the Atlantic rainforests of southeastern and southern Brazil. The 14 species currently described have snout-vent lengths less than 18. mm and are thought to have evolved through miniaturization: an evolutionary process leading to an extremely small adult body size. Here, we present the first comprehensive phylogenetic analysis for Brachycephalus, using a multilocus approach based on two nuclear (Rag-1 and Tyr) and three mitochondrial (Cyt b, 12S, and 16S rRNA) gene regions. Phylogenetic relationships were inferred using a partitioned Bayesian analysis of concatenated sequences and the hierarchical Bayesian method (BEST) that estimates species trees based on the multispecies coalescent model. Individual gene trees showed conflict and also varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. The concatenated gene tree was completely resolved and is identical in topology and degree of statistical support to the individual mtDNA gene tree. On the other hand, the BEST species tree showed reduced significant node support relative to the concatenate tree and recovered a basal trichotomy, although some bipartitions were significantly supported at the tips of the species tree. Comparison of the log likelihoods for the concatenated and BEST trees suggests that the method implemented in BEST explains the multilocus data for Brachycephalus better than the Bayesian analysis of concatenated data. Landmark-based geometric morphometrics revealed marked variation in cranial shape between the species of Brachycephalus. In addition, a statistically significant association was demonstrated between variation in cranial shape and genetic distances estimated from the mtDNA and nuclear loci. Notably, B. ephippium and B. garbeana that are predicted to be sister-species in the individual and concatenated gene trees and the BEST species tree share an evolutionary novelty, the hyperossified dorsal plate. © 2011 Elsevier Inc.