The interpretation of stem diameter–height allometry in trees: biomechanical constraints,neighbour effects,or biased regressions?

Static diameter–height allometry data have been used by many ecologists to demonstrate that diameter should increase at a faster rate than height during tree growth, as predicted by biomechanical models. We review the available evidence and examine the potential problems that arise in the interpretation of this relationship. In particular, we reveal how few studies investigating patterns of diameter–height allometry in trees have adequately controlled for neighbour effects. We also demonstrate how the interpretation of diameter–height allometry has suffered from a lack of uniformity in the selection of regression models, and how the use of least squares regression to estimate allometric scaling exponents can be biased. We conclude that most of the published data on static diameter–height relationships in trees tell us virtually nothing about either age (developmental) effects or neighbour effects; they are completely confounded. Further studies are required to analyse the long-term dynamic growth trajectories of individual trees in relation to local neighbour effects, and greater effort must be made to establish the consistent use of unbiased statistical methods between studies.     

Assessing variation in life‐history tactics within a population using mixture regression models: a practical guide for evolutionary ecologists

Mixed models are now well‐established methods in ecology and evolution because they allow accounting for and quantifying within‐ and between‐individual variation. However, the required normal distribution of the random effects can often be violated by the presence of clusters among subjects, which leads to multi‐modal distributions. In such cases, using what is known as mixture regression models might offer a more appropriate approach. These models are widely used in psychology, sociology, and medicine to describe the diversity of trajectories occurring within a population over time (e.g. psychological development, growth). In ecology and evolution, however, these models are seldom used even though understanding changes in individual trajectories is an active area of research in life‐history studies. Our aim is to demonstrate the value of using mixture models to describe variation in individual life‐history tactics within a population, and hence to promote the use of these models by ecologists and evolutionary ecologists. We first ran a set of simulations to determine whether and when a mixture model allows teasing apart latent clustering, and to contrast the precision and accuracy of estimates obtained from mixture models versus mixed models under a wide range of ecological contexts. We then used empirical data from long‐term studies of large mammals to illustrate the potential of using mixture models for assessing within‐population variation in life‐history tactics. Mixture models performed well in most cases, except for variables following a Bernoulli distribution and when sample size was small. The four selection criteria we evaluated [Akaike information criterion (AIC), Bayesian information criterion (BIC), and two bootstrap methods] performed similarly well, selecting the right number of clusters in most ecological situations. We then showed that the normality of random effects implicitly assumed by evolutionary ecologists when using mixed models was often violated in life‐history data. Mixed models were quite robust to this violation in the sense that fixed effects were unbiased at the population level. However, fixed effects at the cluster level and random effects were better estimated using mixture models. Our empirical analyses demonstrated that using mixture models facilitates the identification of the diversity of growth and reproductive tactics occurring within a population. Therefore, using this modelling framework allows testing for the presence of clusters and, when clusters occur, provides reliable estimates of fixed and random effects for each cluster of the population. In the presence or expectation of clusters, using mixture models offers a suitable extension of mixed models, particularly when evolutionary ecologists aim at identifying how ecological and evolutionary processes change within a population. Mixture regression models therefore provide a valuable addition to the statistical toolbox of evolutionary ecologists. As these models are complex and have their own limitations, we provide recommendations to guide future users.     

Individual variation in the growth of captive infant gorillas

Serial anthropometric data were obtained during the first year of life of six nursery-reared infant gorillas in the Columbus (Ohio) Zoo. Two of the infants are likely to be monozygotic twins as determined by DNA analysis. Growth curves were fitted to serial measures of cephalo-thoracic-abdominal length, arm length, leg length, head circumference, upper arm circumference, and weight from each gorilla, to describe individual patterns of variation in skeletal growth and body composition. Growth in skeletal measures tended to be curvilinear to varying degrees over the first year of life. Body composition varied more than skeletal measures throughout the first year as a consequence of individual health status. Individual growth and body composition variations appear to reflect both genetic and environmental influences in this small sample of captive infant gorillas.     

Methods for analysis of unbalanced,longitudinal, growth data

We describe an approach to analysis of growth that does not depend on assumptions about the underlying functional growth pattern and that allows for multiple observations arising from individual-specific, irregularly spaced data. We produce estimated growth curves for predefined subject groups by using LOWESS, a nonparametric smoothing algorithm. We describe how statistical significance of curve features may be evaluated by using the “jackknife,” a sample re-use method; this technique can be used to assess differences between subject groups. We then obtain residuals at each data point by reference to the estimated curve. Consistency of residuals is evaluated as a characteristic of individual subjects, and in the presence of individual consistency, relative size-for-age is then scored by the average residual for each individual. This allows study of relationships between relative size and other individual characteristics such as birth order, dominance rank, or age of maturation. Finally, we indicate flexibility of these methods and alternatives, propose uses related to other questions about growth, and suggest potential applications to variables other than body size. Appendices demonstrate application of the LOWESS and jackknife algorithms to the problem of testing sex differences in growth. © 1992 Wiley-Liss, Inc.     

On using integral projection models to generate demographically driven predictions of species' distributions: development and validation using sparse data

Knowledge of species' geographic distributions is critical for understanding and forecasting population dynamics, responses to environmental change, biodiversity patterns, and conservation planning. While many suggestive correlative occurrence models have been used to these ends, progress lies in understanding the underlying population biology that generates patterns of range dynamics. Here, we show how to use a limited quantity of demographic data to produce demographic distribution models (DDMs) using integral projection models for size‐structured populations. By modeling survival, growth, and fecundity using regression, integral projection models can interpolate across missing size data and environmental conditions to compensate for limited data. To accommodate the uncertainty associated with limited data and model assumptions, we use Bayesian models to propagate uncertainty through all stages of model development to predictions. DDMs have a number of strengths: 1) DDMs allow a mechanistic understanding of spatial occurrence patterns; 2) DDMs can predict spatial and temporal variation in local population dynamics; 3) DDMs can facilitate extrapolation under altered environmental conditions because one can evaluate the consequences for individual vital rates. To illustrate these features, we construct DDMs for an overstory perennial shrub in the Proteaceae family in the Cape Floristic Region of South Africa. We find that the species' population growth rate is limited most strongly by adult survival throughout the range and by individual growth in higher rainfall regions. While the models predict higher population growth rates in the core of the range under projected climates for 2050, they also suggest that the species faces a threat along arid range margins from the interaction of more frequent fire and drying climate. The results (and uncertainties) are helpful for prioritizing additional sampling of particular demographic parameters along these gradients to iteratively refine projections. In the appendices, we provide fully functional R code to perform all analyses.     

Mathematical models of bacterial growth,inhibition and death under combined stress conditions

Summary In this report we review the history of growth theories. We show how classical growth models may be derived as special cases of a generic growth rate equation. We show how growth models may be modified to represent survival data. We use linear combinations of growth and survival models to represent complex growth/survival curves and give practical examples utilizing nonlinear regression analysis. We show that traditional methods of estimating D values are inappropriate for complex, multiphasic growth/survival data. We show how such data may be modeled mathematically and illustrate methods for estimating true D values from such data.     

Statistical analysis of test-day milk yields using random regression models for the comparison of feeding groups during the lactation period

Random regression models are widely used in the field of animal breeding for the genetic evaluation of daily milk yields from different test days. These models are capable of handling different environmental effects on the respective test day, and they describe the characteristics of the course of the lactation period by using suitable covariates with fixed and random regression coefficients. As the numerically expensive estimation of parameters is already part of advanced computer software, modifications of random regression models will considerably grow in importance for statistical evaluations of nutrition and behaviour experiments with animals. Random regression models belong to the large class of linear mixed models. Thus, when choosing a model, or more precisely, when selecting a suitable covariance structure of the random effects, the information criteria of Akaike and Schwarz can be used. In this study, the fitting of random regression models for a statistical analysis of a feeding experiment with dairy cows is illustrated under application of the program package SAS. For each of the feeding groups, lactation curves modelled by covariates with fixed regression coefficients are estimated simultaneously. With the help of the fixed regression coefficients, differences between the groups are estimated and then tested for significance. The covariance structure of the random and subject-specific effects and the serial correlation matrix are selected by using information criteria and by estimating correlations between repeated measurements. For the verification of the selected model and the alternative models, mean values and standard deviations estimated with ordinary least square residuals are used.     

Statistical analysis of test-day milk yields using random regression models for the comparison of feeding groups during the lactation period

Abstract

Random regression models are widely used in the field of animal breeding for the genetic evaluation of daily milk yields from different test days. These models are capable of handling different environmental effects on the respective test day, and they describe the characteristics of the course of the lactation period by using suitable covariates with fixed and random regression coefficients. As the numerically expensive estimation of parameters is already part of advanced computer software, modifications of random regression models will considerably grow in importance for statistical evaluations of nutrition and behaviour experiments with animals. Random regression models belong to the large class of linear mixed models. Thus, when choosing a model, or more precisely, when selecting a suitable covariance structure of the random effects, the information criteria of Akaike and Schwarz can be used. In this study, the fitting of random regression models for a statistical analysis of a feeding experiment with dairy cows is illustrated under application of the program package SAS. For each of the feeding groups, lactation curves modelled by covariates with fixed regression coefficients are estimated simultaneously. With the help of the fixed regression coefficients, differences between the groups are estimated and then tested for significance. The covariance structure of the random and subject-specific effects and the serial correlation matrix are selected by using information criteria and by estimating correlations between repeated measurements. For the verification of the selected model and the alternative models, mean values and standard deviations estimated with ordinary least square residuals are used.     

Functional linear models for region-based association analysis

Regional association analysis is one of the most powerful tools for gene mapping because instead analysis of individual variants it simultaneously considers all variants in the region. Recent development of the models for regional association analysis involves functional data analysis approach. In the framework of this approach, genotypes of variants within region as well as their effects are described by continuous functions. Such approach allows us to use information about both linkage and linkage disequilibrium and reduce the influence of noise and/or observation errors. Here we define a functional linear mixed model to test association on independent and structured samples. We demonstrate how to test fixed and random effects of a set of genetic variants in the region on quantitative trait. Estimation of statistical properties of new methods shows that type I errors are in accordance with declared values and power is high especially for models with fixed effects of genotypes. We suppose that new functional regression linear models facilitate identification of rare genetic variants controlling complex human and animal traits. New methods are implemented in computer software FREGAT which is available for free download at http://mga.bionet.nsc.ru/soft/FREGAT/.     

Analysis of infant growth curves using multivariate adaptive splines

In this paper, we study the effect of cocaine use by a pregnant woman on the growth of her infant after birth. Using a data set from a retrospective study, we found that cocaine use was a marginally significant contributor to the infant growth as measured by bodyweight. From a statistical point of view, the data represent a common, though complex, structure that has received little attention in the statistical literature. To analyze these data, we adopt and further enhance an approach developed recently called MASAL (multivariate adaptive splines for the analysis of longitudinal data). In addition to the fitting of growth curves, we demonstrate particularly how to explore and estimate the underlying covariance structures for the longitudinal data that were collected from a rather irregular schedule.     

Implementation of Rao's one-sample polynomial growth curve model using SAS

Longitudinal data are frequently treated with the classic analysis of variance and regression models. However, these models assume independence of observations. Hoel (1964) demonstrated that the use of least-squares methods on intercorrelated serial observations results in the rejection of the null hypothesis much too frequently. Although appropriate models for analyzing longitudinal data have been available for quite some time, they have remained inaccessible due to cumbersome matrix manipulations. We implement Rao's (1959) one-sample polynomial growth curve model using the programming capability and matrix language of SAS, which involves testing the goodness-of-fit and calculation of confidence bands for polynomial growth curves fit to data at equally spaced time points. Confidence intervals for the parameters themselves are also computed. The method and program (presented in the Appendix) are illustrated with examples involving mandibular ramus height in 12 young male rhesus monkeys. The data set, which spans a 4 year period (yearly observations), is fit adequately by a quadratic equation. The data spanning a 2 year period (half-year observations) are fit adequately by the linear equation. These examples illustrate the considerable widening of confidence bands that occurs when polynomial equations having more terms than are needed to meet the goodness-of-fit requirement are considered.     

Local environment and density-dependent feedbacks determine population growth in a forest herb

Linking spatial variation in environmental factors to variation in demographic rates is essential for a mechanistic understanding of the dynamics of populations. However, we still know relatively little about such links, partly because feedbacks via intraspecific density make them difficult to observe in natural populations. We conducted a detailed field study and investigated simultaneous effects of environmental factors and the intraspecific density of individuals on the demography of the herb Lathyrus vernus. In regression models of vital rates we identified effects associated with spring shade on survival and growth, while density was negatively correlated with these vital rates. Density was also negatively correlated with average individual size in the study plots, which is consistent with self-thinning. In addition, average plant sizes were larger than predicted by density in plots that were less shaded by the tree canopy, indicating an environmentally determined carrying capacity. A size-structured integral projection model based on the vital rate regressions revealed that the identified effects of shade and density were strong enough to produce differences in stable population sizes similar to those observed in the field. The results illustrate how the local environment can determine dynamics of populations and that intraspecific density may have to be more carefully considered in studies of plant demography and population viability analyses of threatened species. We conclude that demographic approaches incorporating information about both density and key environmental factors are powerful tools for understanding the processes that interact to determine population dynamics and abundances.     

Estimating polymorphic growth curve sets with nonchronological data

1. When we collect the growth curves of many individuals, orderly variation in the curves is often observed rather than a completely random mixture of various curves. Small individuals may exhibit similar growth curves, but the curves differ from those of large individuals, whereby the curves gradually vary from small to large individuals. It has been recognized that after standardization with the asymptotes, if all the growth curves are the same (anamorphic growth curve set), the growth curve sets can be estimated using nonchronological data; otherwise, that is, if the growth curves are not identical after standardization with the asymptotes (polymorphic growth curve set), this estimation is not feasible. However, because a given set of growth curves determines the variation in the observed data, it may be possible to estimate polymorphic growth curve sets using nonchronological data.
2. In this study, we developed an estimation method by deriving the likelihood function for polymorphic growth curve sets. The method involves simple maximum likelihood estimation. The weighted nonlinear regression and least‐squares method after the log‐transform of the anamorphic growth curve sets were included as special cases.
3. The growth curve sets of the height of cypress (Chamaecyparis obtusa) and larch (Larix kaempferi) trees were estimated. With the model selection process using the AIC and likelihood ratio test, the growth curve set for cypress was found to be polymorphic, whereas that for larch was found to be anamorphic. Improved fitting using the polymorphic model for cypress is due to resolving underdispersion (less dispersion in real data than model prediction).
4. The likelihood function for model estimation depends not only on the distribution type of asymptotes, but the definition of the growth curve set as well. Consideration of these factors may be necessary, even if environmental explanatory variables and random effects are introduced.

     

Lifetime growth in wild meerkats: incorporating life history and environmental factors into a standard growth model

Lifetime records of changes in individual size or mass in wild animals are scarce and, as such, few studies have attempted to model variation in these traits across the lifespan or to assess the factors that affect them. However, quantifying lifetime growth is essential for understanding trade-offs between growth and other life history parameters, such as reproductive performance or survival. Here, we used model selection based on information theory to measure changes in body mass over the lifespan of wild meerkats, and compared the relative fits of several standard growth models (monomolecular, von Bertalanffy, Gompertz, logistic and Richards). We found that meerkats exhibit monomolecular growth, with the best model incorporating separate growth rates before and after nutritional independence, as well as effects of season and total rainfall in the previous nine months. Our study demonstrates how simple growth curves may be improved by considering life history and environmental factors, which may be particularly relevant when quantifying growth patterns in wild populations.     

Random regression analyses using B-spline functions to model growth of Nellore cattle

The objective of this study was to estimate (co)variance components using random regression on B-spline functions to weight records obtained from birth to adulthood. A total of 82 064 weight records of 8145 females obtained from the data bank of the Nellore Breeding Program (PMGRN/Nellore Brazil) which started in 1987, were used. The models included direct additive and maternal genetic effects and animal and maternal permanent environmental effects as random. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of age (cubic regression) were considered as random covariate. The random effects were modeled using B-spline functions considering linear, quadratic and cubic polynomials for each individual segment. Residual variances were grouped in five age classes. Direct additive genetic and animal permanent environmental effects were modeled using up to seven knots (six segments). A single segment with two knots at the end points of the curve was used for the estimation of maternal genetic and maternal permanent environmental effects. A total of 15 models were studied, with the number of parameters ranging from 17 to 81. The models that used B-splines were compared with multi-trait analyses with nine weight traits and to a random regression model that used orthogonal Legendre polynomials. A model fitting quadratic B-splines, with four knots or three segments for direct additive genetic effect and animal permanent environmental effect and two knots for maternal additive genetic effect and maternal permanent environmental effect, was the most appropriate and parsimonious model to describe the covariance structure of the data. Selection for higher weight, such as at young ages, should be performed taking into account an increase in mature cow weight. Particularly, this is important in most of Nellore beef cattle production systems, where the cow herd is maintained on range conditions. There is limited modification of the growth curve of Nellore cattle with respect to the aim of selecting them for rapid growth at young ages while maintaining constant adult weight.     

Is my study system good enough? A case study for identifying maternal effects

In this paper, we demonstrate how simulation studies can be used to answer questions about identifiability and consequences of omitting effects from a model. The methodology is presented through a case study where identifiability of genetic and/or individual (environmental) maternal effects is explored. Our study system is a wild house sparrow (Passer domesticus) population with known pedigree. We fit pedigree‐based (generalized) linear mixed models (animal models), with and without additive genetic and individual maternal effects, and use deviance information criterion (DIC) for choosing between these models. Pedigree and R‐code for simulations are available. For this study system, the simulation studies show that only large maternal effects can be identified. The genetic maternal effect (and similar for individual maternal effect) has to be at least half of the total genetic variance to be identified. The consequences of omitting a maternal effect when it is present are explored. Our results indicate that the total (genetic and individual) variance are accounted for. When an individual (environmental) maternal effect is omitted from the model, this only influences the estimated (direct) individual (environmental) variance. When a genetic maternal effect is omitted from the model, both (direct) genetic and (direct) individual variance estimates are overestimated.     

The home-range concept: are traditional estimators still relevant with modern telemetry technology?

Recent advances in animal tracking and telemetry technology have allowed the collection of location data at an ever-increasing rate and accuracy, and these advances have been accompanied by the development of new methods of data analysis for portraying space use, home ranges and utilization distributions. New statistical approaches include data-intensive techniques such as kriging and nonlinear generalized regression models for habitat use. In addition, mechanistic home-range models, derived from models of animal movement behaviour, promise to offer new insights into how home ranges emerge as the result of specific patterns of movements by individuals in response to their environment. Traditional methods such as kernel density estimators are likely to remain popular because of their ease of use. Large datasets make it possible to apply these methods over relatively short periods of time such as weeks or months, and these estimates may be analysed using mixed effects models, offering another approach to studying temporal variation in space-use patterns. Although new technologies open new avenues in ecological research, our knowledge of why animals use space in the ways we observe will only advance by researchers using these new technologies and asking new and innovative questions about the empirical patterns they observe.     

Breeding value accuracy estimates for growth traits using random regression and multi-trait models in Nelore cattle

We quantified the potential increase in accuracy of expected breeding value for weights of Nelore cattle, from birth to mature age, using multi-trait and random regression models on Legendre polynomials and B-spline functions. A total of 87,712 weight records from 8144 females were used, recorded every three months from birth to mature age from the Nelore Brazil Program. For random regression analyses, all female weight records from birth to eight years of age (data set I) were considered. From this general data set, a subset was created (data set II), which included only nine weight records: at birth, weaning, 365 and 550 days of age, and 2, 3, 4, 5, and 6 years of age. Data set II was analyzed using random regression and multi-trait models. The model of analysis included the contemporary group as fixed effects and age of dam as a linear and quadratic covariable. In the random regression analyses, average growth trends were modeled using a cubic regression on orthogonal polynomials of age. Residual variances were modeled by a step function with five classes. Legendre polynomials of fourth and sixth order were utilized to model the direct genetic and animal permanent environmental effects, respectively, while third-order Legendre polynomials were considered for maternal genetic and maternal permanent environmental effects. Quadratic polynomials were applied to model all random effects in random regression models on B-spline functions. Direct genetic and animal permanent environmental effects were modeled using three segments or five coefficients, and genetic maternal and maternal permanent environmental effects were modeled with one segment or three coefficients in the random regression models on B-spline functions. For both data sets (I and II), animals ranked differently according to expected breeding value obtained by random regression or multi-trait models. With random regression models, the highest gains in accuracy were obtained at ages with a low number of weight records. The results indicate that random regression models provide more accurate expected breeding values than the traditionally finite multi-trait models. Thus, higher genetic responses are expected for beef cattle growth traits by replacing a multi-trait model with random regression models for genetic evaluation. B-spline functions could be applied as an alternative to Legendre polynomials to model covariance functions for weights from birth to mature age.     

Marginal Models Via GLS: A Convenient Yet Neglected Tool for the Analysis of Correlated Data in the Behavioural Sciences

Behavioural research often produces data that have a complicated structure. For instance, data can represent repeated observations of the same individual and suffer from heteroscedasticity as well as other technical snags. The regression analysis of such data is often complicated by the fact that the observations (response variables) are mutually correlated. The correlation structure can be quite complex and might or might not be of direct interest to the user. In any case, one needs to take correlations into account (e.g. by means of random‐effect specification) in order to arrive at correct statistical inference (e.g. for construction of the appropriate test or confidence intervals). Over the last decade, such data have been more and more frequently analysed using repeated‐measures ANOVA and mixed‐effects models. Some researchers invoke the heavy machinery of mixed‐effects modelling to obtain the desired population‐level (marginal) inference, which can be achieved by using simpler tools – namely marginal models. This paper highlights marginal modelling (using generalized least squares [GLS] regression) as an alternative method. In various concrete situations, such marginal models can be based on fewer assumptions and directly generate estimates (population‐level parameters) which are of immediate interest to the behavioural researcher (such as population mean). Sometimes, they might be not only easier to interpret but also easier to specify than their competitors (e.g. mixed‐effects models). Using five examples from behavioural research, we demonstrate the use, advantages, limits and pitfalls of marginal and mixed‐effects models implemented within the functions of the ‘nlme’ package in R.     

Estimation of fitness from energetics and life‐history data: An example using mussels

Changing environments have the potential to alter the fitness of organisms through effects on components of fitness such as energy acquisition, metabolic cost, growth rate, survivorship, and reproductive output. Organisms, on the other hand, can alter aspects of their physiology and life histories through phenotypic plasticity as well as through genetic change in populations (selection). Researchers examining the effects of environmental variables frequently concentrate on individual components of fitness, although methods exist to combine these into a population level estimate of average fitness, as the per capita rate of population growth for a set of identical individuals with a particular set of traits. Recent advances in energetic modeling have provided excellent data on energy intake and costs leading to growth, reproduction, and other life‐history parameters; these in turn have consequences for survivorship at all life‐history stages, and thus for fitness. Components of fitness alone (performance measures) are useful in determining organism response to changing conditions, but are often not good predictors of fitness; they can differ in both form and magnitude, as demonstrated in our model. Here, we combine an energetics model for growth and allocation with a matrix model that calculates population growth rate for a group of individuals with a particular set of traits. We use intertidal mussels as an example, because data exist for some of the important energetic and life‐history parameters, and because there is a hypothesized energetic trade‐off between byssus production (affecting survivorship), and energy used for growth and reproduction. The model shows exactly how strong this trade‐off is in terms of overall fitness, and it illustrates conditions where fitness components are good predictors of actual fitness, and cases where they are not. In addition, the model is used to examine the effects of environmental change on this trade‐off and on both fitness and on individual fitness components.