Source author record

Leonardo Rojas-Nandayapa

Leonardo Rojas-Nandayapa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.PR

Catalog footprint

What is connected

4works

1topics

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.

preprint2016arXiv

Efficient simulation for dependent rare events with applications to extremes

We consider the general problem of estimating probabilities which arise as a union of dependent events. We propose a flexible series of estimators for such probabilities, and describe variance reduction schemes applied to the proposed estimators. We derive efficiency results of the estimators in rare-event settings, in particular those associated with extremes. Finally, we examine the performance of our estimators in a numerical example.

preprint2015arXiv

Approximating the Laplace transform of the sum of dependent lognormals

Let $(X_1, \dots, X_n)$ be multivariate normal, with mean vector $\boldsymbolμ$ and covariance matrix $\boldsymbolΣ$, and $S_n=\mathrm{e}^{X_1}+\cdots+\mathrm{e}^{X_n}$. The Laplace transform ${\cal L}(θ)=\mathbb{E}\mathrm{e}^{-θS_n} \propto \int \exp\{-h_θ(\boldsymbol{x})\} \,\mathrm{d} \boldsymbol{x}$ is represented as $\tilde{\cal L}(θ)I(θ)$, where $\tilde{\cal L}(θ)$ is given in closed-form and $I(θ)$ is the error factor ($\approx 1$). We obtain $\tilde{\cal L}(θ)$ by replacing $h_θ(\boldsymbol{x})$ with a second order Taylor expansion around its minimiser $\boldsymbol{x}^*$. An algorithm for calculating the asymptotic expansion of $\boldsymbol{x}^*$ is presented, and it is shown that $I(θ)\to 1$ as $θ\to\infty$. A variety of numerical methods for evaluating $I(θ)$ are discussed, including Monte Carlo with importance sampling and quasi-Monte Carlo. Numerical examples (including Laplace transform inversion for the density of $S_n$) are also given.

preprint2014arXiv

Exponential Family Techniques for the Lognormal Left Tail

Let $X$ be lognormal$(μ,σ^2)$ with density $f(x)$, let $θ>0$ and define ${L}(θ)=E e^{-θX}$. We study properties of the exponentially tilted density (Esscher transform) $f_θ(x) =e^{-θx}f(x)/{L}(θ)$, in particular its moments, its asymptotic form as $θ\to\infty$ and asymptotics for the Cramér function; the asymptotic formulas involve the Lambert W function. This is used to provide two different numerical methods for evaluating the left tail probability of lognormal sum $S_n=X_1+\cdots+X_n$: a saddlepoint approximation and an exponential twisting importance sampling estimator. For the latter we demonstrate the asymptotic consistency by proving logarithmic efficiency in terms of the mean square error. Numerical examples for the c.d.f.\ $F_n(x)$ and the p.d.f.\ $f_n(x)$ of $S_n$ are given in a range of values of $σ^2,n,x$ motivated from portfolio Value-at-Risk calculations.

preprint2013arXiv

Non-Existence of Stabilizing Policies for the Critical Push-Pull Network and Generalizations

The push-pull queueing network is a simple example in which servers either serve jobs or generate new arrivals. It was previously conjectured that there is no policy that makes the network positive recurrent (stable) in the critical case. We settle this conjecture and devise a general sufficient condition for non-stabilizability of queueing networks which is based on a linear martingale and further applies to generalizations of the push-pull network.