Stage duration distributions in matrix population models.
Ontology highlight
ABSTRACT:
SUBMITTER: Okuyama T
PROVIDER: S-EPMC6144974 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A matrix population model is a convenient tool for summarizing per capita survival and reproduction rates (collectively vital rates) of a population and can be used for calculating an asymptotic finite population growth rate (λ) and generation time. These two pieces of information can be used for determining the status of a threatened species. The use of stage-structured population models has increased in recent years, and the vital rates in such models are often estimated using a life table analysis. However, potential bias introduced when converting age-structured vital rates estimated from a life table into parameters for a stage-structured population model has not been assessed comprehensively. The objective of this study was to investigate the performance of methods for such conversions using simulated life histories of organisms. The underlying models incorporate various types of life history and true population growth rates of varying levels. The performance was measured by comparing differences in λ and the generation time calculated using the Euler-Lotka equation, age-structured population matrices, and several stage-structured population matrices that were obtained by applying different conversion methods. The results show that the discretization of age introduces only small bias in λ or generation time. Similarly, assuming a fixed age of maturation at the mean age of maturation does not introduce much bias. However, aggregating age-specific survival rates into a stage-specific survival rate and estimating a stage-transition rate can introduce substantial bias depending on the organism's life history type and the true values of λ. In order to aggregate survival rates, the use of the weighted arithmetic mean was the most robust method for estimating λ. Here, the weights are given by survivorship curve after discounting with λ. To estimate a stage-transition rate, matching the proportion of individuals transitioning, with λ used for discounting the rate, was the best approach. However, stage-structured models performed poorly in estimating generation time, regardless of the methods used for constructing the models. Based on the results, we recommend using an age-structured matrix population model or the Euler-Lotka equation for calculating λ and generation time when life table data are available. Then, these age-structured vital rates can be converted into a stage-structured model for further analyses.
| S-EPMC5660883 | biostudies-literature
Project description:In learning language, children must discover how to interpret the linguistic significance of phonetic variation. On some accounts, receptive phonology is grounded in perceptual learning of phonetic categories from phonetic distributions drawn over the infant's sample of speech. On other accounts, receptive phonology is instead based on phonetic generalizations over the words in the lexicon. Tests of these hypotheses have been rare and indirect, usually making use of idealized estimates of phonetic variation. Here we evaluated these hypotheses, using as our test case English and Dutch toddlers' different interpretation of the lexical significance of vowel duration. Analysis of thousands of vowels of one Dutch and three English mothers' speech suggests that children's language-specific differences in interpretation of vowel duration are likely due to detection of lexically specific patterns, rather than bimodality in raw phonetic distributions.
| S-EPMC6788768 | biostudies-literature
Project description:Deep generative models of molecules have grown immensely in popularity, trained on relevant datasets, these models are used to search through chemical space. The downstream utility of generative models for the inverse design of novel functional compounds, depends on their ability to learn a training distribution of molecules. The most simple example is a language model that takes the form of a recurrent neural network and generates molecules using a string representation. Since their initial use, subsequent work has shown that language models are very capable, in particular, recent research has demonstrated their utility in the low data regime. In this work, we investigate the capacity of simple language models to learn more complex distributions of molecules. For this purpose, we introduce several challenging generative modeling tasks by compiling larger, more complex distributions of molecules and we evaluate the ability of language models on each task. The results demonstrate that language models are powerful generative models, capable of adeptly learning complex molecular distributions. Language models can accurately generate: distributions of the highest scoring penalized LogP molecules in ZINC15, multi-modal molecular distributions as well as the largest molecules in PubChem. The results highlight the limitations of some of the most popular and recent graph generative models- many of which cannot scale to these molecular distributions.
| S-EPMC9174447 | biostudies-literature
Project description:History matters when individual prior conditions contain important information about the fate of individuals. We present a general framework for demographic models which incorporates the effects of history on population dynamics. The framework incorporates prior condition into the i-state variable and includes an algorithm for constructing the population projection matrix from information on current state dynamics as a function of prior condition. Three biologically motivated classes of prior condition are included: prior stages, linear functions of current and prior stages, and equivalence classes of prior stages. Taking advantage of the matrix formulation of the model, we show how to calculate sensitivity and elasticity of any demographic outcome. Prior condition effects are a source of inter-individual variation in vital rates, i.e., individual heterogeneity. As an example, we construct and analyze a second-order model of Lathyrus vernus, a long-lived herb. We present population growth rate, the stable population distribution, the reproductive value vector, and the elasticity of λ to changes in the second-order transition rates. We quantify the contribution of prior conditions to the total heterogeneity in the stable population of Lathyrus using the entropy of the stable distribution.
| S-EPMC6445492 | biostudies-literature
Project description:The placental weight ratio (PWR) is a health indicator that reflects the balance between fetal and placental growth. The PWR is defined as the placental weight divided by the birth weight, and it changes across gestation. Its ranges are not well established. We aimed to establish PWR distributions by gestational age and to investigate whether the PWR distributions vary by fetal growth adequacy, small, average, and large for gestational age (SGA, AGA, and LGA). The data came from a hospital based retrospective cohort, using all births at two London, Ontario hospitals in the past 10 years. All women who delivered a live singleton infant between 22 and 42 weeks of gestation were included (n = 41441). Nonparametric quantile regression was used to fit the curves. The results demonstrate decreasing PWR and dispersion, with increasing gestational age. A higher proportion of SGA infants have extreme PWRs than AGA and LGA, especially at lower gestational ages. On average, SGA infants had higher PWRs than AGA and LGA infants. The overall curves offer population standards for use in research studies. The curves stratified by fetal growth adequacy are the first of their kind, and they demonstrate that PWR differs for SGA and LGA infants.
| S-EPMC4033358 | biostudies-literature
Project description:Many human diseases are characterized by multiple stages of progression. While the typical sequence of disease progression can be identified, there may be large individual variations among patients. Identifying mean stage durations and their variations is critical for statistical hypothesis testing needed to determine if treatment is having a significant effect on the progression, or if a new therapy is showing a delay of progression through a multistage disease. In this paper we focus on two methods for extracting stage duration statistics from longitudinal datasets: an extension of the linear regression technique, and a counting algorithm. Both are non-iterative, non-parametric and computationally cheap methods, which makes them invaluable tools for studying the epidemiology of diseases, with a goal of identifying different patterns of progression by using bioinformatics methodologies. Here we show that the regression method performs well for calculating the mean stage durations under a wide variety of assumptions, however, its generalization to variance calculations fails under realistic assumptions about the data collection procedure. On the other hand, the counting method yields reliable estimations for both means and variances of stage durations. Applications to Alzheimer disease progression are discussed.
| S-EPMC3233570 | biostudies-literature
Project description:Models of molecular evolution tend to be overly simplistic caricatures of biology that are prone to assigning high probabilities to biologically implausible DNA or protein sequences. Here, we explore how to construct time-reversible evolutionary models that yield stationary distributions of sequences that match given target distributions. By adopting comparatively realistic target distributions,evolutionary models can be improved. Instead of focusing on estimating parameters, we concentrate on the population genetic implications of these models. Specifically, we obtain estimates of the product of effective population size and relative fitness difference of alleles. The approach is illustrated with two applications to protein-coding DNA. In the first, a codon-based evolutionary model yields a stationary distribution of sequences, which, when the sequences are translated,matches a variable-length Markov model trained on human proteins. In the second, we introduce an insertion-deletion model that describes selectively neutral evolutionary changes to DNA. We then show how to modify the neutral model so that its stationary distribution at the amino acid level can match a profile hidden Markov model, such as the one associated with the Pfam database.
| S-EPMC2607412 | biostudies-literature
Project description:The measurement of phytoplankton distributions in ocean ecosystems provides the basis for elucidating the influences of physical processes on plankton dynamics. Technological advances allow for measurement of phytoplankton data to greater resolution, displaying high spatial variability. In conventional mathematical models, the mean value of the measured variable is approximated to compare with the model output, which may misinterpret the reality of planktonic ecosystems, especially at the microscale level. To consider intermittency of variables, in this work, a new modelling approach to the planktonic ecosystem is applied, called the closure approach. Using this approach for a simple nutrient-phytoplankton model, we have shown how consideration of the fluctuating parts of model variables can affect system dynamics. Also, we have found a critical value of variance of overall fluctuating terms below which the conventional non-closure model and the mean value from the closure model exhibit the same result. This analysis gives an idea about the importance of the fluctuating parts of model variables and about when to use the closure approach. Comparisons of plot of mean versus standard deviation of phytoplankton at different depths, obtained using this new approach with real observations, give this approach good conformity.
| S-EPMC4008380 | biostudies-literature
Project description:BACKGROUND: Matrix models are widely used to study the dynamics and demography of populations. An important but overlooked issue is how the number of individuals sampled influences estimates of the population growth rate (lambda) calculated with matrix models. Even unbiased estimates of vital rates do not ensure unbiased estimates of lambda-Jensen's Inequality implies that even when the estimates of the vital rates are accurate, small sample sizes lead to biased estimates of lambda due to increased sampling variance. We investigated if sampling variability and the distribution of sampling effort among size classes lead to biases in estimates of lambda. METHODOLOGY/PRINCIPAL FINDINGS: Using data from a long-term field study of plant demography, we simulated the effects of sampling variance by drawing vital rates and calculating lambda for increasingly larger populations drawn from a total population of 3842 plants. We then compared these estimates of lambda with those based on the entire population and calculated the resulting bias. Finally, we conducted a review of the literature to determine the sample sizes typically used when parameterizing matrix models used to study plant demography. CONCLUSIONS/SIGNIFICANCE: We found significant bias at small sample sizes when survival was low (survival = 0.5), and that sampling with a more-realistic inverse J-shaped population structure exacerbated this bias. However our simulations also demonstrate that these biases rapidly become negligible with increasing sample sizes or as survival increases. For many of the sample sizes used in demographic studies, matrix models are probably robust to the biases resulting from sampling variance of vital rates. However, this conclusion may depend on the structure of populations or the distribution of sampling effort in ways that are unexplored. We suggest more intensive sampling of populations when individual survival is low and greater sampling of stages with high elasticities.
| S-EPMC2518208 | biostudies-literature
Project description:Actomyosin stress fibers (SFs) enable cells to exert traction on planar extracellular matrices (ECMs) by tensing focal adhesions (FAs) at the cell-ECM interface. Although it is widely appreciated that the spatial and temporal distribution of these tensile forces play key roles in polarity, motility, fate choice, and other defining cell behaviors, virtually nothing is known about how an individual SF quantitatively contributes to tensile loads borne by specific molecules within associated FAs. We address this key open question by using femtosecond laser ablation to sever single SFs in cells while tracking tension across vinculin using a molecular optical sensor. We show that disruption of a single SF reduces tension across vinculin in FAs located throughout the cell, with enriched vinculin tension reduction in FAs oriented parallel to the targeted SF. Remarkably, however, some subpopulations of FAs exhibit enhanced vinculin tension upon SF irradiation and undergo dramatic, unexpected transitions between tension-enhanced and tension-reduced states. These changes depend strongly on the location of the severed SF, consistent with our earlier finding that different SF pools are regulated by distinct myosin activators. We critically discuss the extent to which these measurements can be interpreted in terms of whole-FA tension and traction and propose a model that relates SF tension to adhesive loads and cell shape stability. These studies represent the most direct and high-resolution intracellular measurements of SF contributions to tension on specific FA proteins to date and offer a new paradigm for investigating regulation of adhesive complexes by cytoskeletal force.
| S-EPMC3711198 | biostudies-literature