IMPROVED BAYESIAN ANALYSIS OF METAPOPULATION DATA WITH AN APPLICATION TO A TREE FROG METAPOPULATION

Metapopulation models are important tools to predict whether a species can persist in a landscape consisting of habitat patches. Here a Bayesian method is presented for estimating parameters of such models from data on patch occupancy in one or more years. Earlier methods were either ad hoc, produced only point estimates, or could only use turnover information. The new method is based on the assumption of quasi-stationarity, which enables it to use not only turnover data, but also snapshot data. Being Bayesian, the method produces reliable information about the uncertainty of the parameters and model-based predictions in the form of posterior distributions. It is computationally demanding, but considerably faster than a recently developed Bayesian method extended beyond turnover data. The method is compared with existing methods (placed in a Bayesian framework) by fitting an extended incidence function model to a data set on a tree frog metapopulation with many missing values and by predicting its viability, mean occupancy, and turnover rate after 100 yr. Corresponding Editor: O. N. Bjornstad.

[1]  Paul Opdam,et al.  European Nuthatch Metapopulations in a Fragmented Agricultural Landscape , 1991 .

[2]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[3]  W. Wong,et al.  The calculation of posterior distributions by data augmentation , 1987 .

[4]  F. Gosselin,et al.  Test of mathematical assumptions behind the 'incidence function' estimation process of metapopulations' dynamic parameters. , 1999, Mathematical biosciences.

[5]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[6]  L. Wasserman,et al.  The Selection of Prior Distributions by Formal Rules , 1996 .

[7]  N. Metropolis,et al.  The Monte Carlo method. , 1949 .

[8]  Atte Moilanen,et al.  PATCH OCCUPANCY MODELS OF METAPOPULATION DYNAMICS: EFFICIENT PARAMETER ESTIMATION USING IMPLICIT STATISTICAL INFERENCE , 1999 .

[9]  Otso Ovaskainen,et al.  The metapopulation capacity of a fragmented landscape , 2000, Nature.

[10]  I. Hanski,et al.  Habitat Deterioration, Habitat Destruction, and Metapopulation Persistence in a Heterogenous Landscape , 1997, Theoretical population biology.

[11]  J. Copas.  Regression, Prediction and Shrinkage , 1983 .

[12]  P. Amarasekare,et al.  Allee Effects in Metapopulation Dynamics , 1998, The American Naturalist.

[13]  I. Hanski.  A Practical Model of Metapopulation Dynamics , 1994 .

[14]  Hannu Toivonen,et al.  BAYESIAN ANALYSIS OF METAPOPULATION DATA , 2002 .

[15]  Atte Moilanen,et al.  The Quantitative Incidence Function Model and Persistence of an Endangered Butterfly Metapopulation , 1996 .

[16]  Ilkka Hanski,et al.  Coexistence of Competitors in Patchy Environment , 1983 .

[17]  James H. Brown,et al.  Turnover Rates in Insular Biogeography: Effect of Immigration on Extinction , 1977 .

[18]  Atte Moilanen,et al.  The equilibrium assumption in estimating the parameters of metapopulation models. , 2000 .

[19]  R. Etienne.  Local populations of different sizes, mechanistic rescue effect and patch preference in the Levins metapopulation model , 2000, Bulletin of mathematical biology.

[20]  R. Levins.  Some Demographic and Genetic Consequences of Environmental Heterogeneity for Biological Control , 1969 .

[21]  Claire C. Vos,et al.  Incidence function modelling and conservation of the tree frog Hyla arborea in the Netherlands , 2000 .