Importance sampling for maxima on trees

Abstract

We consider the distributional fixed-point equation: R D= Q ( i=1N Ci Ri ), where the \Ri\ are i.i.d.~copies of R, independent of the vector (Q, N, \Ci\), where N ∈ N, Q, \Ci\ ≥ 0 and P(Q > 0) > 0. By setting W = R, Xi = Ci, Y = Q it is equivalent to the high-order Lindley equation W D= \ Y, \, 1 ≤ i ≤ N (Xi + Wi) \. It is known that under Kesten assumptions, P(W > t) H e-α t, t ∞, where α>0 solves the Cram\'er-Lundberg equation E [ Σj=1N Ci α ] = E[ Σi=1N eα Xi ] = 1. The main goal of this paper is to provide an explicit representation for P(W > t), which can be directly connected to the underlying weighted branching process where W is constructed and that can be used to construct unbiased and strongly efficient estimators for all t. Furthermore, we show how this new representation can be directly analyzed using Alsmeyer's Markov renewal theorem, yielding an alternative representation for the constant H. We provide numerical examples illustrating the use of this new algorithm.

0

Discussion (0)

Sign in to join the discussion.

Loading comments…