CRAN Package vitality

For a single source for the vitality models, please refer to the CRAN package vitality: Fitting Routines for the Vitality Family of Mortality Models and the companion documentation.

The vitality-based survival model is a parametric model for relating stressors and environmental properties to organism survivorship. Vitality is an abstract property which changes in response to organism experience. An organism's resistance to disease, level of stress, behavior, success and failure in feeding, frequency of predator attack, mating, parental care, and habitat choice all induce incremental changes in vitality. Mortality occurs when an organism's vitality reaches zero. Mortality can also occur independent of vitality through accidental-based mortality; catastrophic events occurring equally within a population and independent of the different past histories of the individuals. The model survival distribution is then the product of the probability of survival according to the organism's vitality and the probability of avoiding accidental mortality.

### Four (4) Parameter Vitality Model

Li, T. and J.J. Anderson. 2009. The vitality model: A way to understand population survival and demographic heterogeneity. *Theoretical Population Biology*. 76(2): 118-131.

#### Abstract

A four-parameter model describing mortality as the first passage of an abstract measure of survival capacity, vitality, is developed and used to explore four classic problems in demography: (1) medfly demographic paradox, (2) effect of diet restriction on longevity, (3) cross-life stage effects on survival curves and (4) mortality plateaus. The model quantifies the sources of mortality in these classical problems into vitality-dependent and independent parts, and characterizes the vitality-dependent part in terms of initial and evolving heterogeneities. Three temporal scales express the balance of these factors: a time scale of death from senescence, a time scale of accidental mortality and a crossover time between evolving vs. initial heterogeneity. The examples demonstrate how the first-passage approach provides a unique and informative perspective into the processes that shape the survival curves of populations.

Anderson, J.J., M.C. Gildea, D.W. Williams, and T. Li. 2008. Linking Growth, Survival, and Heterogeneity through Vitality. *The American Naturalist*. 171: E20-E43.

#### Abstract

We model the cross-stage effect of juvenile growth on future cohort survival with vitality, a single stochastic measure of an organism's survival capacity that results in death when it reaches 0. In this construct, the distribution of vitality at the end of a growth treatment stage, which is a measure of survival capacity heterogeneity, determines a cohort's susceptibility to starvation in a subsequent challenge stage. The model predicts that the treatment-stage duration and mass gain determine the mean and variance of the initial vitality distribution of the challenge stage, which in turn determine the effect of a challenge-stage stressor on survival. Studies linking the effect of juvenile growth on time to starvation for chinook salmon and yellow perch are compared to model predictions. The feasibility of predicting survival and heterogeneity in overwintering fish populations from first-year growth is considered. Some limitations and potential extensions of the model to other scenarios are discussed.

#### 4 Parameter Model R Code and Data, *June 2009*

The supplemental R script file vitality.gaussian.R contains all functions necessary to run the MLE parameter fitting routine for the vitality-based survival model (Li and Anderson 2009) with an initial Gaussian distribution. The functions are written in the R programming language. Most of the functions are very similar to those written by Salinger et al. (2003) in S-plus for fitting the vitality model with a Dirac delta function initial distribution. The fitting routines differ in how they derive initial parameter estimates for the Newton-Ralphson method.

- Vitality model 4 parameter fitting code: vitality.gaussian.R (
*release date*: 29 June 2009) - Instructions for using parameter fitting code: README
- Sample function call and data: Page 6 of README

### Three (3) Parameter Vitality Model

Anderson JJ. Characterizing System Failure Curves with Vitality. Army Science Conference, 26th [Internet]. 2008 . Available from: www.cbr.washington.edu/sites/default/files/papers/anderson/2008_army26thscienceconf.pdf

#### Abstract

Identifying the causes of failure is fundamental to improving the survivability of Army systems whether they are remote sensors, complex vehicles, or networks. The shape of a repair frequency curve of a complex system is determined by variability in the initial quality of its components, variability in the system’s wear and maintenance, and random failures independent of the condition of the components. The contributions of these sources of failure are modeled using vitality theory, which describes population survival curves in terms of the stochastic decline of a hidden Markov process (vitality) to an absorbing boundary representing death. By reinterpreting the vitality model, the population becomes the complex system, the individuals become the system’s components, vitality becomes the remaining amount of component wear prior to failure, and absorption into the boundary represents component failure and replacement. Fitting the model to the repair frequency curve of a 180 ton rear dump truck, the shape of the repair curve is quantitatively partitioned into three factors: the operating environment, quality control in manufacture, and variability in maintenance. This information may be useful in identifying causes of system failure.

Salinger, D.H., J.J. Anderson, and O.S. Hamel. 2003. A parameter estimation routine for the vitality-based survival model. *Ecological Modelling* 166(3): 287-294.

#### Abstract

An algorithm is described which determines model parameters for the vitality-based survival model (Anderson 2000). The algorithm uses an MLE estimator on interval-based survival data and returns the vitality model parameters: the rate of vitality loss, the variability in the rate of vitality loss and the rate of accidental mortality. Standard errors of the parameter estimates and a Goodness-of-fit measure (via a Pearson's C test) are also returned. Issues of right censored data and study start time are addressed.

Anderson, J.J. 2000. A vitality-based model relating stressors and environmental properties to organism survival. *Ecological Monographs* 70(3): 445-470.

#### Abstract

A survivorship curve is shaped by the differential survivability of the organisms within the population, and a change in a survivorship curve with a stressor reflects the differential response of the organisms to the stressor. Quantifying this linkage in a simple, rigorous way is valuable for characterizing the response of populations to stressors and ultimately for understanding the evolutionary selection of individuals exposed to stressors. To quantify this stressor-individual-population linkage with as few parameters as possible, I present a simple mechanistic model describing organism survival in terms of age-dependent and age-independent mortality rates. The age-independent rate is represented by a Poisson process. For the age-dependent rate, a concept of vitality is defined, and mortality occurs when an organism's vitality is exhausted. The loss of vitality over age is represented by a continuous Brownian-motion process, the Weiner process; vitality-related mortality occurs when the random process reaches the boundary of zero vitality. The age at which vitality-related mortality occurs is represented by the Weiner-process probability distribution for first-arrival time. The basic model has three rate parameters: the rate of accidental mortality, the mean rate of vitality loss, and the variability in the rate of vitality loss. These rates are related to body mass, environmental conditions, and xenobiotic stressors, resulting in a model that characterizes intrinsic and extrinsic factors that control a population's survival and the distribution of vitality of its individuals. The model assumes that these factors contribute to the rate parameters additively and linearly.

The model is evaluated with case studies across a range of species exposed to natural and xenobiotic stressors. The mean rate of vitality loss generally is the dominant factor in determining the shape of survival curves under optimal conditions. Xenobiotic stressors add to the mean rate in proportion to the strength of the stressor. The base, or intrinsic, vitality loss rate is proportional to the -1/3 power of adult body mass across a range of iteroparous species. The increase in vitality loss rate with a xenobiotic stressor can be a function of body mass according to the allometric relationship of the organism structures affected by the stressor. The model's applicability to dose–response studies is illustrated with case studies including natural stressors (temperature, feeding interval, and population density) and xenobiotic stressors (organic and inorganic toxicants). The model provides a way to extrapolate the impact of stressors measured in one environment to another environment; by characterizing how stressors alter the vitality probability distribution, it can quantify the degree to which a stressor differentiates members of a population.

#### 3 Parameter Model S-PLUS Code and Data, *March 2003, December 2004 update*

The algorithm determines model parameters for the vitality-based survival model (Anderson 2000). The algorithm uses an MLE estimator on interval-based survival data and returns the vitality model parameters: the rate of vitality loss, the variability in the rate of vitality loss and the rate of accidental mortality. Standard errors of the parameter estimates and a Goodness-of-fit measure (via a Pearson's C test) are also returned. Issues of right censored data and study start time are addressed.

- Vitality model 3 parameter fitting code: VitalityModelFitting.SSC (
*release date*: 3 December 2004) - Instructions for using parameter fitting code: README
- Sample rainbow trout data file: RT.sdd
- Sample function call: RTtestScript.ssc

#### S-PLUS Code Change Log

- 12/3/2004 - Safeguarded against input of incremental mortality data with zero mortality in initial timestep.