Effects of Sample Size on Estimates of Population Growth Rates Calculated with Matrix Models

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

Egg distributions of insect parasitoids: modelling and analysis of temporal data with host density dependence

A simple numerical procedure is presented for the problem of estimating the parameters of models for the distribution of eggs oviposited in a host. The modelling is extended to incorporate both host density and time dependence to produce a remarkably parsimonious structure with only seven parameters to describe a data set of over 3,000 observations. This is further refined using a mixed model to accommodate several large outliers. Both models show that the level of superparasitism declines with increasing host density, and the rate declines over time. It is proposed that the differing behaviours represented by the mixed model may reflect a balance between behavioural strategies of different selective benefit