Source author record

Henrik Hult

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

ResearcherUnclaimed source record

math.PR Computation math.AP physics.soc-ph Populations and Evolution q-fin.PR

Catalog footprint

What is connected

8works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Asymptotic behaviour of sampling and transition probabilities in coalescent models under selection and parent dependent mutations

The results in this paper provide new information on asymptotic properties of classical models: the neutral Kingman coalescent under a general finite-alleles, parent-dependent mutation mechanism, and its generalisation, the ancestral selection graph. Several relevant quantities related to these fundamental models are not explicitly known when mutations are parent dependent. Examples include the probability that a sample taken from a population has a certain type configuration, and the transition probabilities of their block counting jump chains. In this paper, asymptotic results are derived for these quantities, as the sample size goes to infinity. It is shown that the sampling probabilities decay polynomially in the sample size with multiplying constant depending on the stationary density of the Wright-Fisher diffusion and that the transition probabilities converge to the limit of frequencies of types in the sample.

preprint2020arXiv

Estimates of the proportion of SARS-CoV-2 infected individuals in Sweden

In this paper a Bayesian SEIR model is studied to estimate the proportion of the population infected with SARS-CoV-2, the virus responsible for COVID-19. To capture heterogeneity in the population and the effect of interventions to reduce the rate of epidemic spread, the model uses a time-varying contact rate, whose logarithm has a Gaussian process prior. A Poisson point process is used to model the occurrence of deaths due to COVID-19 and the model is calibrated using data of daily death counts in combination with a snapshot of the the proportion of individuals with an active infection, performed in Stockholm in late March. The methodology is applied to regions in Sweden. The results show that the estimated proportion of the population who has been infected is around 13.5% in Stockholm, by 2020-05-15, and ranges between 2.5% - 15.6% in the other investigated regions. In Stockholm where the peak of daily death counts is likely behind us, parameter uncertainty does not heavily influence the expected daily number of deaths, nor the expected cumulative number of deaths. It does, however, impact the estimated cumulative number of infected individuals. In the other regions, where random sampling of the number of active infections is not available, parameter sharing is used to improve estimates, but the parameter uncertainty remains substantial.

preprint2020arXiv

Exact simulation of coupled Wright-Fisher diffusions

In this paper an exact rejection algorithm for simulating paths of the coupled Wright-Fisher diffusion is introduced. The coupled Wright-Fisher diffusion is a family of multidimensional Wright-Fisher diffusions that have drifts depending on each other through a coupling term and that find applications in the study of interacting genes' networks as those encountered in studies of antibiotic resistance. Our algorithm uses independent neutral Wright-Fisher diffusions as candidate proposals, which can be sampled exactly by means of existing algorithms and are only needed at a finite number of points. Once a candidate is accepted, the remaining of the path can be recovered by sampling from a neutral multivariate Wright-Fisher bridge, for which we also provide an exact sampling strategy. The technique relies on a modification of the alternating series method and extends existing algorithms that are currently available for the one-dimensional case. Finally, the algorithm's complexity is derived and its performance demonstrated in a simulation study.

preprint2016arXiv

Exact and efficient simulation of tail probabilities of heavy-tailed infinite series

We develop an efficient simulation algorithm for computing the tail probabilities of the infinite series $S = \sum_{n \geq 1} a_n X_n$ when random variables $X_n$ are heavy-tailed. As $S$ is the sum of infinitely many random variables, any simulation algorithm that stops after simulating only fixed, finitely many random variables is likely to introduce a bias. We overcome this challenge by rewriting the tail probability of interest as a sum of a random number of telescoping terms, and subsequently developing conditional Monte Carlo based low variance simulation estimators for each telescoping term. The resulting algorithm is proved to result in estimators that a) have no bias, and b) require only a fixed, finite number of replications irrespective of how rare the tail probability of interest is. Thus, by combining a traditional variance reduction technique such as conditional Monte Carlo with more recent use of auxiliary randomization to remove bias in a multi-level type representation, we develop an efficient and unbiased simulation algorithm for tail probabilities of $S$. These have many applications including in analysis of financial time-series and stochastic recurrence equations arising in models in actuarial risk and population biology.

preprint2014arXiv

Large deviations for weighted empirical measures arising in importance sampling

Importance sampling is a popular method for efficient computation of various properties of a distribution such as probabilities, expectations, quantiles etc. The output of an importance sampling algorithm can be represented as a weighted empirical measure, where the weights are given by the likelihood ratio between the original distribution and the sampling distribution. In this paper the efficiency of an importance sampling algorithm is studied by means of large deviations for the weighted empirical measure. The main result, which is stated as a Laplace principle for the weighted empirical measure arising in importance sampling, can be viewed as a weighted version of Sanov's theorem. The main theorem is applied to quantify the performance of an importance sampling algorithm over a collection of subsets of a given target set as well as quantile estimates. The analysis yields an estimate of the sample size needed to reach a desired precision as well as of the reduction in cost for importance sampling compared to standard Monte Carlo.

preprint2014arXiv

Min-max representations of viscosity solutions of Hamilton-Jacobi equations and applications in rare-event simulation

In this paper a duality relation between the Mañé potential and Mather's action functional is derived in the context of convex and state-dependent Hamiltonians. The duality relation is used to obtain min-max representations of viscosity solutions of first order Hamilton-Jacobi equations. These min-max representations naturally suggest class\-es of subsolutions of Hamilton-Jacobi equations that arise in the theory of large deviations. The subsolutions, in turn, are good candidates for designing efficient rare-event simulation algorithms.

preprint2013arXiv

A simple time-consistent model for the forward density process

In this paper a simple model for the evolution of the forward density of the future value of an asset is proposed. The model allows for a straightforward initial calibration to option prices and has dynamics that are consistent with empirical findings from option price data. The model is constructed with the aim of being both simple and realistic, and avoid the need for frequent re-calibration. The model prices of $n$ options and a forward contract are expressed as time-varying functions of an $(n+1)$-dimensional Brownian motion and it is investigated how the Brownian trajectory can be determined from the trajectories of the price processes. An approach based on particle filtering is presented for determining the location of the driving Brownian motion from option prices observed in discrete time. A simulation study and an empirical study of call options on the S&P 500 index illustrates that the model provides a good fit to option price data.

preprint2012arXiv

Markov chain Monte Carlo for computing rare-event probabilities for a heavy-tailed random walk

In this paper a method based on a Markov chain Monte Carlo (MCMC) algorithm is proposed to compute the probability of a rare event. The conditional distribution of the underlying process given that the rare event occurs has the probability of the rare event as its normalizing constant. Using the MCMC methodology a Markov chain is simulated, with that conditional distribution as its invariant distribution, and information about the normalizing constant is extracted from its trajectory. The algorithm is described in full generality and applied to the problem of computing the probability that a heavy-tailed random walk exceeds a high threshold. An unbiased estimator of the reciprocal probability is constructed whose normalized variance vanishes asymptotically. The algorithm is extended to random sums and its performance is illustrated numerically and compared to existing importance sampling algorithms.