Source author record

Max O. Souza

Max O. Souza appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2016arXiv

Conservative parabolic problems: non-degenerated theory and degenerated examples from population dynamics

We consider partial differential equations (PDE) of drift-diffusion type in the unit interval, supplemented by either two conservation laws or by a conservation law and a further boundary condition. We treat two different cases: (i) uniform parabolic problems; (ii) degenerated problems at the boundaries. The former can be treated in a very general and complete way, much as the traditional boundary value problems. The latter, however, bring new issues, and we restrict our study to a class of forward Kolmogorov equations that arise naturally when the corresponding stochastic process has either one or two absorbing boundaries. These equations are treated by means of a uniform parabolic regularisation, which then yields a measure solution in the vanishing regularisation limit that is unique. Two prototypical problems from population dynamics are treated in detail.

preprint2015arXiv

A Hedged Monte Carlo Approach to Real Option Pricing

In this work we are concerned with valuing optionalities associated to invest or to delay investment in a project when the available information provided to the manager comes from simulated data of cash flows under historical (or subjective) measure in a possibly incomplete market. Our approach is suitable also to incorporating subjective views from management or market experts and to stochastic investment costs. It is based on the Hedged Monte Carlo strategy proposed by Potters et al (2001) where options are priced simultaneously with the determination of the corresponding hedging. The approach is particularly well-suited to the evaluation of commodity related projects whereby the availability of pricing formulae is very rare, the scenario simulations are usually available only in the historical measure, and the cash flows can be highly nonlinear functions of the prices.

preprint2015arXiv

Fixation in large populations: a continuous view of a discrete problem

We study fixation in large, but finite, populations with two types, and dynamics governed by birth-death processes. By considering a restricted class of such processes, we derive a continuous approximation for the probability of fixation that is valid beyond the weak-selection (WS) limit. From the continuous approximations, we then obtain asymptotic approximations for evolutionary dynamics with at most one equilibrium, in the selection-driven regime, that does not preclude a weak-selection regime. As an application, we study the fixation pattern when the infinite population limit has an interior Evolutionary Stable Strategy (ESS): (i) we show that the fixation pattern for the Hawk and Dove game satisfies what we term the one-half law: if the Evolutionary Stable Strategy (ESS) is outside a small interval around $\sfrac{1}{2}$, the fixation is of dominance type; (ii) we also show that, outside of the weak-selection regime, the long-term dynamics of large populations can have very little resemblance to the infinite population case; in addition, we also present results for the case of two equilibria. Finally, we present continuous restatements valid for large populations of two classical concepts naturally defined in the discrete case: (i) the definition of an $\textsc{ESS}_N$ strategy; (ii) the definition of a risk-dominant strategy. We then present two applications of these restatements: (i) we obtain an asymptotic definition valid in the quasi-neutral regime that recovers both the one-third law under linear fitness and the generalised one-third law for $d$-player games; (ii) we extend the ideas behind the (generalised) one-third law outside the quasi-neutral regime and, as a generalisation, we introduce the concept of critical-frequency; (iii) we recover the classification of risk-dominant strategies for $d$-player games.

preprint2013arXiv

Multiscale Analysis for a Vector-Borne Epidemic Model

Traditional studies about disease dynamics have focused on global stability issues, due to their epidemiological importance. We study a classical SIR-SI model for arboviruses in two different directions: we begin by describing an alternative proof of previously known global stability results by using only a Lyapunov approach. In the sequel, we take a different view and we argue that vectors and hosts can have very distinctive intrinsic time-scales, and that such distinctiveness extends to the disease dynamics. Under these hypothesis, we show that two asymptotic regimes naturally appear: the fast host dynamics and the fast vector dynamics. The former regime yields, at leading order, a SIR model for the hosts, but with a rational incidence rate. In this case, the vector disappears from the model, and the dynamics is similar to a directly contagious disease. The latter yields a SI model for the vectors, with the hosts disappearing from the model. Numerical results show the performance of the approximation, and a rigorous proof validates the reduced models.

preprint2013arXiv

The frequency-dependent Wright-Fisher model: diffusive and non-diffusive approximations

We study a class of processes that are akin to the Wright-Fisher model, with transition probabilities weighted in terms of the frequency-dependent fitness of the population types. By considering an approximate weak formulation of the discrete problem, we are able to derive a corresponding continuous weak formulation for the probability density. Therefore, we obtain a family of partial differential equations (PDE) for the evolution of the probability density, and which will be an approximation of the discrete process in the joint large population, small time-steps and weak selection limit. If the fitness functions are sufficiently regular, we can recast the weak formulation in a more standard formulation, without any boundary conditions, but supplemented by a number of conservation laws. The equations in this family can be purely diffusive, purely hyperbolic or of convection-diffusion type, with frequency dependent convection. The particular outcome will depend on the assumed scalings. The diffusive equations are of the degenerate type; using a duality approach, we also obtain a frequency dependent version of the Kimura equation without any further assumptions. We also show that the convective approximation is related to the replicator dynamics and provide some estimate of how accurate is the convective approximation, with respect to the convective-diffusion approximation. In particular, we show that the mode, but not the expected value, of the probability distribution is modelled by the replicator dynamics. Some numerical simulations that illustrate the results are also presented.

preprint2009arXiv

A non-standard evolution problem arising in population genetics

We study the evolution of the probability density of an asexual, one locus population under natural selection and random evolution. This evolution is governed by a Fokker-Planck equation with degenerate coefficients on the boundaries, supplemented by a pair of conservation laws. It is readily shown that no classical or standard weak solution definition yields solvability of the problem. We provide an appropriate definition of weak solution for the problem, for which we show existence and uniqueness. The solution displays a very distinctive structure and, for large time, we show convergence to a unique stationary solution that turns out to be a singular measure supported at the endpoints. An exponential rate of convergence to this steady state is also proved.

preprint2009arXiv

Global stability for a class of virus models with CTL immune response and antigenic variation

We study the global stability of a class of models for in-vivo virus dynamics, that take into account the CTL immune response and display antigenic variation. This class includes a number of models that have been extensively used to model HIV dynamics. We show that models in this class are globally asymptotically stable, under mild hypothesis, by using appropriate Lyapunov functions. We also characterise the stable equilibrium points for the entire biologically relevant parameter range. As a byproduct, we are able to determine what is the diversity of the persistent strains.

preprint2009arXiv

The SIR epidemic model from a PDE point of view

We present a derivation of the classical SIR model through a mean-field approximation from a discrete version of SIR. We then obtain a hyperbolic forward Kolmogorov equation, and show that its projected characteristics recover the standard SIR model. Moreover, we show that the long time limit of the evolution will be a Dirac measure. The exact position will depend on the well-know $R_0$ parameter, and it will be supported on the corresponding stable SIR equilibrium.

preprint2008arXiv

From discrete to continuous evolution models: a unifying approach to drift-diffusion and replicator dynamics

We study the large population limit of the Moran process, assuming weak-selection, and for different scalings. Depending on the particular choice of scalings, we obtain a continuous model that may highlight the genetic-drift (neutral evolution) or natural selection; for one precise scaling, both effects are present. For the scalings that take the genetic-drift into account, the continuous model is given by a singular diffusion equation, together with two conservation laws that are already present at the discrete level. For scalings that take into account only natural selection, we obtain a hyperbolic singular equation that embeds the Replicator Dynamics and satisfies only one conservation law. The derivation is made in two steps: a formal one, where the candidate limit model is obtained, and a rigorous one, where convergence of the probability density is proved. Additional results on the fixation probabilities are also presented.

preprint2006arXiv

Discrete versus continuous models in evolutionary dynamics: from simple to simpler -- and even simpler -- models

There are many different models--both continuous and discrete--used to describe gene mutation fixation. In particular, the Moran process, the Kimura equation and the replicator dynamics are all well known models, that might lead to different conclusions. We present a discussion of a unified framework to embrace all these models, in the large population regime.