PARENT SESSION
Tuesday, August 8, 8:00-11:30 am
OOS 4 - The modern paradigm in population ecology: stochastic, statistical, and inferential
Ballroom D, Ballroom Level, Cook Convention Center
Organized by: EE Holmes (eli.holmes@noaa.gov), C Jordan, and B Dennis

This session features contemporary research on stochastic ecological dynamics and estimation that is fundamentally changing the way we think about and make inferences about ecological processes.

Estimating density dependence, process noise, and observation error.

Dennis, Brian^*,1, ¹ University of Idaho, Moscow, ID

ABSTRACT- In this presentation I describe a discrete time, stochastic population model with density dependence, environmental-type process noise and lognormal observation or sampling error. The model, a stochastic version of the Gompertz model, can be transformed into a linear Gaussian state space model (Kalman filter) for convenient fitting to time series data. The model has a multivariate normal likelihood function and is simple enough for a variety of uses ranging from theoretical study of parameter estimation issues to routine data analyses in population monitoring. A special case of the model is the discrete time, stochastic exponential growth model (density independence) with environmental-type process error and lognormal observation error. Two methods for estimating parameters in the Gompertz state space model, maximum likelihood based on time series observations and restricted maximum likelihood based on first differences, are compared with computer simulations. Both offer adequate statistical properties. Because the likelihood function is identical to a repeated measures analysis of variance model with a random time effect, parameter estimates can be calculated using PROC MIXED of SAS. Data sets from the Breeding Bird Survey provide illustrative analyses. For one data set, the fitted model suggests that over 70% of the noise in the population's growth rate is due to observation error. The model describes the autocovariance properties of the data especially well. While observation error and process noise variance parameters can both be estimated from one time series, multimodal likelihood functions can and do occur. For data arising from the model, the statistically consistent parameter estimates do not necessarily correspond to the global maximum in the likelihood function. Maximization, simulation, and bootstrapping programs must accomodate the phenomenon of multimodal likelihood functions to produce statistically valid results.

Key words: state space model, density dependence, statistical ecology