Predator-prey models: an application for the plankton dynamics of lake Geneva

In this paper we present a hierarchical Bayesian analysis for a predator-prey model applied to ecology considering the use of Markov Chain Monte Carlo methods. We consider the introduction of a random effect in the model and the presence of a covariate vector. An application to ecology is considered using a data set related to the plankton dynamics of lake Geneva for the year 1990. We also discuss some aspects of discrimination of the proposed models.

Reaction-diffusion stochastic lattice model for a predator-prey system

We have the purpose of analyzing the effect of explicit diffusion processes in a predator-prey stochastic lattice model. More precisely we wish to investigate the possible effects due to diffusion upon the thresholds of coexistence of species, i. e., the possible changes in the transition between the active state and the absorbing state devoid of predators. To accomplish this task we have performed time dependent simulations and dynamic mean-field approximations. Our results indicate that the diffusive process can enhance the species coexistence.

CRITICAL BEHAVIOR AND THRESHOLD OF COEXISTENCE OF A PREDATOR-PREY STOCHASTIC MODEL IN A 2D LATTICE

We investigate the critical behavior of a stochastic lattice model describing a predator-prey system. By means of Monte Carlo procedure we simulate the model defined on a regular square lattice and determine the threshold of species coexistence, that is, the critical phase boundaries related to the transition between an active state, where both species coexist and an absorbing state where one of the species is extinct. A finite size scaling analysis is employed to determine the order parameter, order parameter fluctuations, correlation length and the critical exponents. Our numerical results for the critical exponents agree with those of the directed percolation universality class. We also check the validity of the hyperscaling relation and present the data collapse curves.

The stochastic nature of predator-prey cycles

We study by numerical simulations the time correlation function of a stochastic lattice model describing the dynamics of coexistence of two interacting biological species that present time cycles in the number of species individuals. Its asymptotic behavior is shown to decrease in time as a sinusoidal exponential function from which we extract the dominant eigenvalue of the evolution operator related to the stochastic dynamics showing that it is complex with the imaginary part being the frequency of the population cycles. The transition from the oscillatory to the nonoscillatory behavior occurs when the asymptotic behavior of the time correlation function becomes a pure exponential, that is, when the real part of the complex eigenvalue equals a real eigenvalue. We also show that the amplitude of the undamped oscillations increases with the square root of the area of the habitat as ordinary random fluctuations. (C) 2009 Elsevier B.V. All rights reserved.

The role of predator overlap in the robustness and extinction of a four species predator-prey network

Predators and preys often form species networks with asymmetric patterns of interaction. We study the dynamics of a four species network consisting of two weakly connected predator-prey pairs. We focus our analysis on the effects of the cross interaction between the predator of the first pair and the prey of the second pair. This is an example where the predator overlap, which is the proportion of predators that a given prey shares with other preys, is not uniform across the network due to asymmetries in patterns of interaction. We explore the behavior of the system under different interaction strengths and study the dynamics of survival and extinction. In particular, we consider situations in which the four species have initial populations lower than their long-term equilibrium, simulating catastrophic situations in which their abundances are reduced due to human action or environmental change. We show that, under these reduced initial conditions, and depending on the strength of the cross interaction, the populations tend to oscillate before re-equilibrating, disturbing the community equilibrium and sometimes reaching values that are only a small fraction of the equilibrium population, potentially leading to their extinction. We predict that, contrary to one`s intuition, the most likely scenario is the extinction of the less predated preys. (C) 2010 Elsevier B.V. All rights reserved.

Role of noise in population dynamics cycles

Noise is an intrinsic feature of population dynamics and plays a crucial role in oscillations called phase-forgetting quasicycles by converting damped into sustained oscillations. This function of noise becomes evident when considering Langevin equations whose deterministic part yields only damped oscillations. We formulate here a consistent and systematic approach to population dynamics, leading to a Fokker-Planck equation and the associate Langevin equations in accordance with this conceptual framework, founded on stochastic lattice-gas models that describe spatially structured predator-prey systems. Langevin equations in the population densities and predator-prey pair density are derived in two stages. First, a birth-and-death stochastic process in the space of prey and predator numbers and predator-prey pair number is obtained by a contraction method that reduces the degrees of freedom. Second, a van Kampen expansion in the inverse of system size is then performed to get the Fokker-Planck equation. We also study the time correlation function, the asymptotic behavior of which is used to characterize the transition from the cyclic coexistence of species to the ordinary coexistence.

Time correlation function in systems with two coexisting biological species

We study a stochastic lattice model describing the dynamics of coexistence of two interacting biological species. The model comprehends the local processes of birth, death, and diffusion of individuals of each species and is grounded on interaction of the predator-prey type. The species coexistence can be of two types: With self-sustained coupled time oscillations of population densities and without oscillations. We perform numerical simulations of the model on a square lattice and analyze the temporal behavior of each species by computing the time correlation functions as well as the spectral densities. This analysis provides an appropriate characterization of the different types of coexistence. It is also used to examine linked population cycles in nature and in experiment.

Effect of crab size and habitat type on the locomotory activity of juvenile shore crabs, Carcinus maenas

Post-settlement processes are a major focus in the study of the dynamics of marine populations and communities. Post-settlement movement of juveniles is an important, but often ignored, process which affects local predator-prey and competitive interactions. We used benthic suction sampling and pitfall traps to examine density and locomotory activity of Carcinus maenas juveniles in different intertidal habitat types in the Rio Mira Estuary, Portugal, to better understand intra-specific interactions in a system where density-dependent processes are known to regulate population dynamics. As expected, significantly higher densities of juvenile shore crabs were found from bare mud compared to densely vegetated habitats. At the time of sampling, small and intermediate stages together outnumbered by far the larger juveniles. Conversely, larger crabs were much more frequent than smaller ones in traps. A locomotory index (LI), i.e. the ratio between crab catch in pitfall traps and their density within their moving range, is proposed as a measure of movement. LI analyses indicated that: (1) movement is an order of magnitude higher in large than small juveniles and much higher in sparse than dense vegetation cover; (2) activity of small juveniles is mostly crepuscular, regardless of vegetation cover; and (3) movement of large juveniles is very limited in dense Zostera patches, but very high in sparsely vegetated areas, during the day and night. These results suggest that small juveniles are relatively protected under dense vegetation cover due to lower mobility of larger crabs, and provide evidence of temporal segregation of activity windows between juvenile crabs of different sizes, which may be a key mechanism to reduce cannibalism and therefore increase the carrying capacity of nursery habitats. (C) 2008 Elsevier Ltd. All rights reserved.

Olfactory eavesdropping between two competing stingless bee species

Foragers can improve search efficiency, and ultimately fitness, by using social information: cues and signals produced by other animals that indicate food location or quality. Social information use has been well studied in predator-prey systems, but its functioning within a trophic level remains poorly understood. Eavesdropping, use of signals by unintended recipients, is of particular interest because eavesdroppers may exert selective pressure on signaling systems. We provide the most complete study to date of eavesdropping between two competing social insect species by determining the glandular source and composition of a recruitment pheromone, and by examining reciprocal heterospecific responses to this signal. We tested eavesdropping between Trigona hyalinata and Trigona spinipes, two stingless bee species that compete for floral resources, exhibit a clear dominance hierarchy and recruit nestmates to high-quality food sources via pheromone trails. Gas chromatography-mass spectrometry of T. hyalinata recruitment pheromone revealed six carboxylic esters, the most common of which is octyl octanoate, the major component of T. spinipes recruitment pheromone. We demonstrate heterospecific detection of recruitment pheromones, which can influence heterospecific and conspecific scout orientation. Unexpectedly, the dominant T. hyalinata avoided T. spinipes pheromone in preference tests, while the subordinate T. spinipes showed neither attraction to nor avoidance of T. hyalinata pheromone. We suggest that stingless bees may seek to avoid conflict through their eavesdropping behavior, incorporating expected costs associated with a choice into the decision-making process.

Improvement of Insulin Resistance and Reduction of Cardiovascular Risk Among Obese Patients with Type 2 Diabetes with the Duodenojejunal Bypass Liner

This study aims to evaluate the effectiveness of the duodenojejunal bypass liner (DJBL) in the improvement of insulin resistance and reduction of cardiovascular risk among morbidly obese patients with type 2 diabetes mellitus, using the triglyceride/high-density lipoprotein (HDL) cholesterol ratio, percentage of weight loss, and glycemic control. We used the TG/HDL ratio with a cutoff value of 3.5 to identify patients with insulin resistance. The value of the initial ratio was compared with the ratio obtained 6 months after implantation to evaluate whether an improvement in insulin resistance occurred. We also evaluated the improvement of glycated hemoglobin levels and the weight loss resulted from the use of the device and correlated that with the improvement of the TG/HDL ratio. All patients implanted with the device presented a statistically significant reduction of the HbA1c levels, with most patients (70.3%) obtaining diabetes control with HbA1c levels lower than 7% at the end of the study. All patients also presented a significant weight reduction, with an average loss of 12.6% of their initial weight. We observed an important improvement in insulin resistance and metabolic syndrome, with a significant reduction of the TG/HDL ratio from 5.75 to 4.36 (p < 0.001) and 42.6% of the patients presenting a TG/HDL ratio lower than 3.5 at the end of the study. The DJBL, when used for a period of 6 months, is effective in the control of diabetes, weight loss, improvement of insulin resistance, and decrease of cardiovascular risk among morbidly obese patients with type 2 diabetes mellitus.

Aging and fluctuation-dissipation ratio in a nonequilibrium q-state lattice model

A generalized version of the nonequilibrium linear Glauber model with q states in d dimensions is introduced and analyzed. The model is fully symmetric, its dynamics being invariant under all permutations of the q states. Exact expressions for the two-time autocorrelation and response functions on a d-dimensional lattice are obtained. In the stationary regime, the fluctuation-dissipation theorem holds, while in the transient the aging is observed with the fluctuation-dissipation ratio leading to the value predicted for the linear Glauber model.

A new scale-invariant ratio and finite-size scaling for the stochastic susceptible-infected-recovered model

The critical behavior of the stochastic susceptible-infected-recovered model on a square lattice is obtained by numerical simulations and finite-size scaling. The order parameter as well as the distribution in the number of recovered individuals is determined as a function of the infection rate for several values of the system size. The analysis around criticality is obtained by exploring the close relationship between the present model and standard percolation theory. The quantity UP, equal to the ratio U between the second moment and the squared first moment of the size distribution multiplied by the order parameter P, is shown to have, for a square system, a universal value 1.0167(1) that is the same for site and bond percolation, confirming further that the SIR model is also in the percolation class.

Stochastic lattice gas model describing the dynamics of the SIRS epidemic process

We study a stochastic process describing the onset of spreading dynamics of an epidemic in a population composed of individuals of three classes: susceptible (S), infected (I), and recovered (R). The stochastic process is defined by local rules and involves the following cyclic process: S -> I -> R -> S (SIRS). The open process S -> I -> R (SIR) is studied as a particular case of the SIRS process. The epidemic process is analyzed at different levels of description: by a stochastic lattice gas model and by a birth and death process. By means of Monte Carlo simulations and dynamical mean-field approximations we show that the SIRS stochastic lattice gas model exhibit a line of critical points separating the two phases: an absorbing phase where the lattice is completely full of S individuals and an active phase where S, I and R individuals coexist, which may or may not present population cycles. The critical line, that corresponds to the onset of epidemic spreading, is shown to belong in the directed percolation universality class. By considering the birth and death process we analyze the role of noise in stabilizing the oscillations. (C) 2009 Elsevier B.V. All rights reserved.

Signed likelihood ratio tests in the Birnbaum-Saunders regression model

The Birnbaum-Saunders regression model is becoming increasingly popular in lifetime analyses and reliability studies. In this model, the signed likelihood ratio statistic provides the basis for testing inference and construction of confidence limits for a single parameter of interest. We focus on the small sample case, where the standard normal distribution gives a poor approximation to the true distribution of the statistic. We derive three adjusted signed likelihood ratio statistics that lead to very accurate inference even for very small samples. Two empirical applications are presented. (C) 2010 Elsevier B.V. All rights reserved.

Validation of an experimental polyurethane model for biomechanical studies on implant supported prosthesis - tension tests

OBJECTIVES: The complexity and heterogeneity of human bone, as well as ethical issues, frequently hinder the development of clinical trials. The purpose of this in vitro study was to determine the modulus of elasticity of a polyurethane isotropic experimental model via tension tests, comparing the results to those reported in the literature for mandibular bone, in order to validate the use of such a model in lieu of mandibular bone in biomechanical studies. MATERIAL AND METHODS: Forty-five polyurethane test specimens were divided into 3 groups of 15 specimens each, according to the ratio (A/B) of polyurethane reagents (PU-1: 1/0.5, PU-2: 1/1, PU-3: 1/1.5). RESULTS: Tension tests were performed in each experimental group and the modulus of elasticity values found were 192.98 MPa (SD=57.20) for PU-1, 347.90 MPa (SD=109.54) for PU-2 and 304.64 MPa (SD=25.48) for PU-3. CONCLUSION: The concentration of choice for building the experimental model was 1/1.