Multivariate meta‐analysis with an increasing number of parameters

Meta‐analysis can average estimates of multiple parameters, such as a treatment's effect on multiple outcomes, across studies. Univariate meta‐analysis (UVMA) considers each parameter individually, while multivariate meta‐analysis (MVMA) considers the parameters jointly and accounts for the correlation between their estimates. The performance of MVMA and UVMA has been extensively compared in scenarios with two parameters. Our objective is to compare the performance of MVMA and UVMA as the number of parameters, p, increases. Specifically, we show that (i) for fixed‐effect (FE) meta‐analysis, the benefit from using MVMA can substantially increase as p increases; (ii) for random effects (RE) meta‐analysis, the benefit from MVMA can increase as p increases, but the potential improvement is modest in the presence of high between‐study variability and the actual improvement is further reduced by the need to estimate an increasingly large between study covariance matrix; and (iii) when there is little to no between‐study variability, the loss of efficiency due to choosing RE MVMA over FE MVMA increases as p increases. We demonstrate these three features through theory, simulation, and a meta‐analysis of risk factors for non‐Hodgkin lymphoma.     

Comparison of Linear,Nonlinear and Semiparametric Mixed‐effects Models for Estimating HIV Dynamic Parameters

The potency of antiretroviral agents in AIDS clinical trials can be assessed on the basis of an early viral response such as viral decay rate or change in viral load (number of copies of HIV RNA) of the plasma. Linear, parametric nonlinear, and semiparametric nonlinear mixed‐effects models have been proposed to estimate viral decay rates in viral dynamic models. However, before applying these models to clinical data, a critical question that remains to be addressed is whether these models produce coherent estimates of viral decay rates, and if not, which model is appropriate and should be used in practice. In this paper, we applied these models to data from an AIDS clinical trial of potent antiviral treatments and found significant incongruity in the estimated rates of reduction in viral load. Simulation studies indicated that reliable estimates of viral decay rate were obtained by using the parametric and semiparametric nonlinear mixed‐effects models. Our analysis also indicated that the decay rates estimated by using linear mixed‐effects models should be interpreted differently from those estimated by using nonlinear mixed‐effects models. The semiparametric nonlinear mixed‐effects model is preferred to other models because arbitrary data truncation is not needed. Based on real data analysis and simulation studies, we provide guidelines for estimating viral decay rates from clinical data. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Small area estimation for semicontinuous skewed spatial data: An application to the grape wine production in Tuscany

Linear‐mixed models are frequently used to obtain model‐based estimators in small area estimation (SAE) problems. Such models, however, are not suitable when the target variable exhibits a point mass at zero, a highly skewed distribution of the nonzero values and a strong spatial structure. In this paper, a SAE approach for dealing with such variables is suggested. We propose a two‐part random effects SAE model that includes a correlation structure on the area random effects that appears in the two parts and incorporates a bivariate smooth function of the geographical coordinates of units. To account for the skewness of the distribution of the positive values of the response variable, a Gamma model is adopted. To fit the model, to get small area estimates and to evaluate their precision, a hierarchical Bayesian approach is used. The study is motivated by a real SAE problem. We focus on estimation of the per‐farm average grape wine production in Tuscany, at subregional level, using the Farm Structure Survey data. Results from this real data application and those obtained by a model‐based simulation experiment show a satisfactory performance of the suggested SAE approach.     

Using functional traits to predict species growth trajectories,and cross‐validation to evaluate these models for ecological prediction

Modeling plant growth using functional traits is important for understanding the mechanisms that underpin growth and for predicting new situations. We use three data sets on plant height over time and two validation methods—in‐sample model fit and leave‐one‐species‐out cross‐validation—to evaluate non‐linear growth model predictive performance based on functional traits. In‐sample measures of model fit differed substantially from out‐of‐sample model predictive performance; the best fitting models were rarely the best predictive models. Careful selection of predictor variables reduced the bias in parameter estimates, and there was no single best model across our three data sets. Testing and comparing multiple model forms is important. We developed an R package with a formula interface for straightforward fitting and validation of hierarchical, non‐linear growth models. Our intent is to encourage thorough testing of multiple growth model forms and an increased emphasis on assessing model fit relative to a model's purpose.     

A coefficient of determination (R2) for generalized linear mixed models

Extensions of linear models are very commonly used in the analysis of biological data. Whereas goodness of fit measures such as the coefficient of determination (R²) or the adjusted R² are well established for linear models, it is not obvious how such measures should be defined for generalized linear and mixed models. There are by now several proposals but no consensus has yet emerged as to the best unified approach in these settings. In particular, it is an open question how to best account for heteroscedasticity and for covariance among observations present in residual error or induced by random effects. This paper proposes a new approach that addresses this issue and is universally applicable for arbitrary variance‐covariance structures including spatial models and repeated measures. It is exemplified using three biological examples.     

Competing Risks and Time‐Dependent Covariates

Time‐dependent covariates are frequently encountered in regression analysis for event history data and competing risks. They are often essential predictors, which cannot be substituted by time‐fixed covariates. This study briefly recalls the different types of time‐dependent covariates, as classified by Kalbfleisch and Prentice [The Statistical Analysis of Failure Time Data, Wiley, New York, 2002] with the intent of clarifying their role and emphasizing the limitations in standard survival models and in the competing risks setting. If random (internal) time‐dependent covariates are to be included in the modeling process, then it is still possible to estimate cause‐specific hazards but prediction of the cumulative incidences and survival probabilities based on these is no longer feasible. This article aims at providing some possible strategies for dealing with these prediction problems. In a multi‐state framework, a first approach uses internal covariates to define additional (intermediate) transient states in the competing risks model. Another approach is to apply the landmark analysis as described by van Houwelingen [Scandinavian Journal of Statistics 2007, 34 , 70–85] in order to study cumulative incidences at different subintervals of the entire study period. The final strategy is to extend the competing risks model by considering all the possible combinations between internal covariate levels and cause‐specific events as final states. In all of those proposals, it is possible to estimate the changes/differences of the cumulative risks associated with simple internal covariates. An illustrative example based on bone marrow transplant data is presented in order to compare the different methods.     

Linking demographic and food‐web models to understand management trade‐offs

Alternatives in ecosystem‐based management often differ with respect to trade‐offs between ecosystem values. Ecosystem or food‐web models and demographic models are typically employed to evaluate alternatives, but the approaches are rarely integrated to uncover conflicts between values. We applied multistate models to a capture–recapture dataset on common guillemots Uria aalge breeding in the Baltic Sea to identify factors influencing survival. The estimated relationships were employed together with Ecopath‐with‐Ecosim food‐web model simulations to project guillemot survival under six future scenarios incorporating climate change. The scenarios were based on management alternatives for eutrophication and cod fisheries, issues considered top priority for regional management, but without known direct effects on the guillemot population. Our demographic models identified prey quantity (abundance and biomass of sprat Sprattus sprattus) as the main factor influencing guillemot survival. Most scenarios resulted in projections of increased survival, in the near (2016–2040) and distant (2060–2085) future. However, in the scenario of reduced nutrient input and precautionary cod fishing, guillemot survival was projected to be lower in both future periods due to lower sprat stocks. Matrix population models suggested a substantial decline of the guillemot population in the near future, 24% per 10 years, and a smaller reduction, 1.1% per 10 years, in the distant future. To date, many stakeholders and Baltic Sea governments have supported reduced nutrient input and precautionary cod fishing and implementation is underway. Negative effects on nonfocal species have previously not been uncovered, but our results show that the scenario is likely to negatively impact the guillemot population. Linking model results allowed identifying trade‐offs associated with management alternatives. This information is critical to thorough evaluation by decision‐makers, but not easily obtained by food‐web models or demographic models in isolation. Appropriate datasets are often available, making it feasible to apply a linked approach for better‐informed decisions in ecosystem‐based management.     

Evidence of reduced individual heterogeneity in adult survival of long‐lived species

The canalization hypothesis postulates that the rate at which trait variation generates variation in the average individual fitness in a population determines how buffered traits are against environmental and genetic factors. The ranking of a species on the slow‐fast continuum – the covariation among life‐history traits describing species‐specific life cycles along a gradient going from a long life, slow maturity, and low annual reproductive output, to a short life, fast maturity, and high annual reproductive output – strongly correlates with the relative fitness impact of a given amount of variation in adult survival. Under the canalization hypothesis, long‐lived species are thus expected to display less individual heterogeneity in survival at the onset of adulthood, when reproductive values peak, than short‐lived species. We tested this life‐history prediction by analysing long‐term time series of individual‐based data in nine species of birds and mammals using capture‐recapture models. We found that individual heterogeneity in survival was higher in species with short‐generation time (< 3 years) than in species with long generation time (> 4 years). Our findings provide the first piece of empirical evidence for the canalization hypothesis at the individual level from the wild.     

A two-part mixed-effects pattern-mixture model to handle zero-inflation and incompleteness in a longitudinal setting

Two-part regression models are frequently used to analyze longitudinal count data with excess zeros, where the same set of subjects is repeatedly observed over time. In this context, several sources of heterogeneity may arise at individual level that affect the observed process. Further, longitudinal studies often suffer from missing values: individuals dropout of the study before its completion, and thus present incomplete data records. In this paper, we propose a finite mixture of hurdle models to face the heterogeneity problem, which is handled by introducing random effects with a discrete distribution; a pattern-mixture approach is specified to deal with non-ignorable missing values. This approach helps us to consider overdispersed counts, while allowing for association between the two parts of the model, and for non-ignorable dropouts. The effectiveness of the proposal is tested through a simulation study. Finally, an application to real data on skin cancer is provided.     

Darwin was right: inbreeding depression on male fertility in the Darwin family

Charles Darwin, who was married to his first cousin Emma Wedgwood, was the first experimentalist to demonstrate the adverse effects of inbreeding. He documented the deleterious consequences of self‐fertilization on progeny in numerous plant species, and this research led him to suspect that the health problems of his 10 children, who were very often ill, might have been a consequence of his marriage to his first cousin. Because Darwin's concerns regarding the consequences of cousin marriage on his children even nowadays are considered controversial, we analyzed the potential effects of inbreeding on fertility in 30 marriages of the Darwin–Wedgwood dynasty, including the marriages of Darwin's children, which correspond to the offspring of four cousin marriages and three marriages between unrelated individuals. Analysis of the number of children per woman through zero‐inflated regression models showed a significantly adverse effect of the husband inbreeding coefficient on family size. Furthermore, a statistically significant adverse effect of the husband inbreeding coefficient on reproductive period duration was also detected. To our knowledge, this is the first time that inbreeding depression on male fertility has been detected in humans. Because Darwin's sons had fewer children in comparison to non‐inbred men of the dynasty, our findings give empirical support to Darwin's concerns on the consequences of consanguineous marriage in his own progeny. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 474–483.     

Fibril elongation mechanisms of HET-s prion-forming domain: topological evidence for growth polarity

The prion-forming C-terminal domain of the fungal prion HET-s forms infectious amyloid fibrils at physiological pH. The conformational switch from the nonprion soluble form to the prion fibrillar form is believed to have a functional role, as HET-s in its prion form participates in a recognition process of different fungal strains. On the basis of the knowledge of the high-resolution structure of the prion forming domain HET-s(218-289) in its fibrillar form, we here present a numerical simulation of the fibril growth process, which emphasizes the role of the topological properties of the fibrillar structure. An accurate thermodynamic analysis of the way an intervening HET-s chain is recruited to the tip of the growing fibril suggests that elongation proceeds through a dock and lock mechanism. First, the chain docks onto the fibril by forming the longest β-strands. Then, the re-arrangement in the fibrillar form of all the rest of the molecule takes place. Interestingly, we also predict that one side of the HET-s fibril is more suitable for sustaining its growth with respect to the other. The resulting strong polarity of fibril growth is a consequence of the complex topology of HET-s fibrillar structure, as the central loop of the intervening chain plays a crucially different role in favoring or not the attachment of the C-terminus tail to the fibril, depending on the growth side.     

Mixed‐Effects State‐Space Models for Analysis of Longitudinal Dynamic Systems

Summary The rapid development of new biotechnologies allows us to deeply understand biomedical dynamic systems in more detail and at a cellular level. Many of the subject‐specific biomedical systems can be described by a set of differential or difference equations that are similar to engineering dynamic systems. In this article, motivated by HIV dynamic studies, we propose a class of mixed‐effects state‐space models based on the longitudinal feature of dynamic systems. State‐space models with mixed‐effects components are very flexible in modeling the serial correlation of within‐subject observations and between‐subject variations. The Bayesian approach and the maximum likelihood method for standard mixed‐effects models and state‐space models are modified and investigated for estimating unknown parameters in the proposed models. In the Bayesian approach, full conditional distributions are derived and the Gibbs sampler is constructed to explore the posterior distributions. For the maximum likelihood method, we develop a Monte Carlo EM algorithm with a Gibbs sampler step to approximate the conditional expectations in the E‐step. Simulation studies are conducted to compare the two proposed methods. We apply the mixed‐effects state‐space model to a data set from an AIDS clinical trial to illustrate the proposed methodologies. The proposed models and methods may also have potential applications in other biomedical system analyses such as tumor dynamics in cancer research and genetic regulatory network modeling.     

Testing models of trophic dynamics: The problem of translating from model to nature

Abstract Mathematical models of interacting populations have a prominent position in population and community ecology, but are often criticized for not being testable. The authors reviewed tests of a particular model, the exploitation ecosystem hypothesis as it was formulated in Oksanen et al. (1981), in order to study problems that may be encountered when testing models. A general problem is how to determine if an experimental system should be regarded as within the model's theoretical domain or not. The theoretical domain defines the type of system the model is meant to apply to. It is noted that both liberal and strict domain definitions can be problematic. Most important is that a too liberal domain definition can result in false understanding (i.e. that it is falsely concluded that the processes included in the model are controlling the study system). Other problems encountered were more system‐specific. Equilibrium predictions were tested in experiments that were too short to reach steady state and in several studies ambiguous definitions and measurements of model variables were found such as productivity, biomass and the number trophic levels. It is concluded that a major obstacle when performing tests is the conceptual and methodological problems encountered when translating model abstractions into an empirical reality.     

Incorporating Gravity Model Principles into Disease Mapping

The space‐time variation of disease risk, considering the closeness and availability of care units, is studied by a hierarchical Bayesian model in a gravitational formulation. The pattern of mortality for Hodgkin's lymphoma over the last twenty‐five years in the Tuscany region and the evaluation of the consequences of distance from the nearest radiation‐therapy center are described. Results show a decrease in mortality for Hodgkin's lymphoma and the attraction exercised by each care units.     

Assessing seasonal demographic covariation to understand environmental‐change impacts on a hibernating mammal

Natural populations are exposed to seasonal variation in environmental factors that simultaneously affect several demographic rates (survival, development and reproduction). The resulting covariation in these rates determines population dynamics, but accounting for its numerous biotic and abiotic drivers is a significant challenge. Here, we use a factor‐analytic approach to capture partially unobserved drivers of seasonal population dynamics. We use 40 years of individual‐based demography from yellow‐bellied marmots (Marmota flaviventer) to fit and project population models that account for seasonal demographic covariation using a latent variable. We show that this latent variable, by producing positive covariation among winter demographic rates, depicts a measure of environmental quality. Simultaneously, negative responses of winter survival and reproductive‐status change to declining environmental quality result in a higher risk of population quasi‐extinction, regardless of summer demography where recruitment takes place. We demonstrate how complex environmental processes can be summarized to understand population persistence in seasonal environments.     

A New Method for Estimating Growth Transition Matrices

Summary The vast majority of population models work using age or stage not length but there are many cases where animals cannot be aged sensibly or accurately. For these cases length‐based models form the logical alternative but there has been little work done to develop and compare different methods of estimating growth transition matrices to be used in such models. This article demonstrates how a consistent Bayesian framework for estimating growth parameters and a novel method for constructing length transition matrices accounts for variation in growth in a clear and consistent manner and avoids potential subjective choices required using more established methods. The inclusion of the resultant growth uncertainty in population assessment models and the potential impact on management decisions is also addressed.