BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Zbigniew Palmowski contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We consider a discrete time parallel queue, which is two-queue network, where at each time-slot there is a the same batch arrival to both queues and at each queue there is a random service available. The service law at each time-slot for each queue is different. Let $(Q_n^1, Q_n^2)$ be the queue length at $n$th time-slot. We present several open questions related to the steady-state of this queue.
In this paper we develop the Gerber-Shiu theory for the classic and dual discrete risk processes in a Markovian (regime switching) environment. In particular, by expressing the Gerber-Shiu function in terms of potential measures of an upward (downward) skip-free discrete-time and discrete-space Markov Additive Process (MAP), we derive closed form expressions for the Gerber-Shiu function in terms of the so-called (discrete) $\boldsymbol{W}_v$ and $\boldsymbol{Z}_v$ scale matrices, which were introduced in arXiv:2008.06697. We show that the discrete scale matrices allow for a unified approach for identifying the Gerber-Shiu function as well as the value function of the associated constant dividend barrier problems.
Stochastic fluid-fluid models (SFFMs) offer powerful modeling ability for a wide range of real-life systems of significance. The existing theoretical framework for this class of models is in terms of operator-analytic methods. For the first time, we establish matrix-analytic methods for the efficient analysis of SFFMs. We illustrate the theory with numerical examples.
In this paper we give a solution to the quickest drift change detection problem for a multivariate Lévy process consisting of both continuous (Gaussian) and jump components in the Bayesian approach. We do it for a general 0-modified continuous prior distribution of the change point. Classically, our criterion of optimality is based on a probability of false alarm and an expected delay of the detection, which is then reformulated in terms of a posterior probability of the change point. We find a generator of the posterior probability, which in case of general prior distribution is inhomogeneous in time. The main solving technique uses the optimal stopping theory and is based on solving a certain free-boundary problem. We also construct a Generelized Shiryaev-Roberts statistic, which can be used for applications. The paper is supplemented by two examples, one of which is further used to analyze Polish life tables (after proper calibration) and detect the drift change in the correlated force of mortality of men and women jointly.
We study critical GI/G/1 queues under finite second moment assumptions. We show that the busy period distribution is regularly varying with index half. We also review previously known M/G/1/ and M/M/1 derivations, yielding exact asymptotics as well as a similar derivation for GI/M/1. The busy period asymptotics determine the growth rate of moments of the renewal process counting busy cycles. We further use this to demonstrate a BRAVO phenomenon (Balancing Reduces Asymptotic Variance of Outputs) for the work-output process (namely the busy-time). This yields new insight on the BRAVO effect. A second contribution of the paper is in settling previous conjectured results about GI/G/1 and GI/G/s BRAVO. Previously, infinite buffer BRAVO was generally only settled under fourth-moment assumptions together with an assumption about the tail of the busy-period. In the current paper we strengthen the previous results by reducing to assumptions to existence of $2+ε$ moments.
We consider the sample average of a centered random walk in $\mathbb{R}^d$ with regularly varying step size distribution. For the first exit time from a compact convex set $A$ not containing the origin, we show that its tail is of lognormal type. Moreover, we show that the typical way for a large exit time to occur is by having a number of jumps growing logarithmically in the scaling parameter.
In this paper, we build on the techniques developed in Albrecher et al. (2013), to generate initial-boundary value problems for ruin probabilities of surplus-dependent premium risk processes, under a renewal case scenario, Erlang (2) claim arrivals, and an exponential claims scenario, Erlang (2) claim sizes. Applying the approximation theory of solutions of linear ordinary differential equations developed in Fedoryuk (1993), we derive the asymptotics of the ruin probabilities when the initial reserve tends to infinity. When considering premiums that are {\it linearly} dependent on reserves, representing for instance returns on risk-free investments of the insurance capital, we firstly derive explicit formulas for the ruin probabilities, from which we can easily determine their asymptotics, only to match the ones obtained for general premiums dependent on reserves. We compare them with the asymptotics of the equivalent ruin probabilities when the premium rate is fixed over time, to measure the gain generated by this additional mechanism of binding the premium rates with the amount of reserve own by the insurance company.
In this paper we consider the following optimal stopping problem $$V^ω_{\rm A}(s) = \sup_{τ\in\mathcal{T}} \mathbb{E}_{s}[e^{-\int_0^τω(S_w) dw} g(S_τ)],$$ where the process $S_t$ is a jump-diffusion process, $\mathcal{T}$ is a family of stopping times while $g$ and $ω$ are fixed payoff function and discount function, respectively. In a financial market context, if $g(s)=(K-s)^+$ or $g(s)=(s-K)^+$ and $\mathbb{E}$ is the expectation taken with respect to a martingale measure, $V^ω_{\rm A}(s)$ describes the price of a perpetual American option with a discount rate depending on the value of the asset process $S_t$. If $ω$ is a constant, the above problem produces the standard case of pricing perpetual American options. In the first part of this paper we find sufficient conditions for the convexity of the value function $V^ω_{\rm A}(s)$. This allows us to determine the stopping region as a certain interval and hence we are able to identify the form of $V^ω_{\rm A}(s)$. We also prove a put-call symmetry for American options with asset-dependent discounting. In the case when $S_t$ is a geometric Lévy process we give exact expressions using the so-called omega scale functions introduced in Li and Palmowski (2018). We prove that the analysed value function satisfies the HJB equation and we give sufficient conditions for the smooth fit property as well. Finally, we present a few examples for which we obtain the analytical form of the value function $V^ω_{\rm A}(s)$.
The main objective of this paper is to present an algorithm of pricing perpetual American put options with asset-dependent discounting. The value function of such an instrument can be described as \begin{equation*} V^ω_{\text{A}^{\text{Put}}}(s) = \sup_{τ\in\mathcal{T}} \mathbb{E}_{s}[e^{-\int_0^τω(S_w) dw} (K-S_τ)^{+}], \end{equation*} where $\mathcal{T}$ is a family of stopping times, $ω$ is a discount function and $\mathbb{E}$ is an expectation taken with respect to a martingale measure. Moreover, we assume that the asset price process $S_t$ is a geometric Lévy process with negative exponential jumps, i.e. $S_t = s e^{ζt + σB_t - \sum_{i=1}^{N_t} Y_i}$. The asset-dependent discounting is reflected in the $ω$ function, so this approach is a generalisation of the classic case when $ω$ is constant. It turns out that under certain conditions on the $ω$ function, the value function $V^ω_{\text{A}^{\text{Put}}}(s)$ is convex and can be represented in a closed form; see Al-Hadad and Palmowski (2021). We provide an option pricing algorithm in this scenario and we present exact calculations for the particular choices of $ω$ such that $V^ω_{\text{A}^{\text{Put}}}(s)$ takes a simplified form.
We analyse an additive-increase and multiplicative-decrease (aka growth-collapse) process that grows linearly in time and that experiences downward jumps at Poisson epochs that are (deterministically) proportional to its present position. This process is used for example in modelling of Transmission Control Protocol (TCP) and can be viewed as a particular example of the so-called shot noise model, a basic tool in modeling earthquakes, avalanches and neuron firings. For this process, and also for its reflected versions, we consider one- and two-sided exit problems that concern the identification of the laws of exit times from fixed intervals and half-lines. All proofs are based on a unified first-step analysis approach at the first jump epoch, which allows us to give explicit, yet involved, formulas for their Laplace transforms. All the eight Laplace transforms can be described in terms of two so-called scale functions $Z_{\uparrow}$ and $L_{\uparrow}$. Here $Z_{\uparrow}$ is described in terms of multiple explicit sums, and $L_{\uparrow}$ in terms of an explicit recursion formula. All other Laplace transforms can be obtained from $Z_{\uparrow}$ and $L_{\uparrow}$ by taking limits, derivatives, integrals and combinations of these.
This paper presents an analysis of the stochastic recursion $W_{i+1} = [V_iW_i+Y_i]^+$ that can be interpreted as an autoregressive process of order 1, reflected at 0. We start our exposition by a discussion of the model's stability condition. Writing $Y_i=B_i-A_i$, for independent sequences of non-negative i.i.d.\ random variables $\{A_i\}_{i\in N_0}$ and $\{B_i\}_{i\in N_0}$, and assuming $\{V_i\}_{i\in N_0}$ is an i.i.d. sequence as well (independent of $\{A_i\}_{i\in N_0}$ and $\{B_i\}_{i\in N_0}$), we then consider three special cases: (i) $V_i$ attains negative values only and $B_i$ has a rational LST, (ii) $V_i$ equals a positive value $a$ with certain probability $p\in (0,1)$ and is negative otherwise, and both $A_i$ and $B_i$ have a rational LST, (iii) $V_i$ is uniformly distributed on $[0,1]$, and $A_i$ is exponentially distributed. In all three cases we derive transient and stationary results, where the transient results are in terms of the transform at a geometrically distributed epoch.
This paper describes an application of dynamic programming to determine the optimal strategy for assigning grapes to pressing tanks in one of the largest Portuguese wineries. To date, linear programming has been employed to generate proposed solutions to analogous problems, but this approach lacks robustness and may, in fact, result in severe losses in cases of sudden changes, which frequently occur in weather-dependent wine factories. Hence, we endowed our model with stochasticity, thereby rendering it less vulnerable to such changes. Our analysis, which is based on real-world data, demonstrates that the proposed algorithm is highly efficient and, after calibration, can be used to support winery's decision-making. The solution proposed herein could also be applied in numerous other contexts where production processes rely on outside supplies.
We consider the Lévy model of the perpetual American call and put options with a negative discount rate under Poisson observations. Similar to the continuous observation case as in De Donno et al. [24], the stopping region that characterizes the optimal stopping time is either a half-line or an interval. The objective of this paper is to obtain explicit expressions of the stopping and continuation regions and the value function, focusing on spectrally positive and negative cases. To this end, we compute the identities related to the first Poisson arrival time to an interval via the scale function and then apply those identities to the computation of the optimal strategies. We also discuss the convergence of the optimal solutions to those in the continuous observation case as the rate of observation increases to infinity. Numerical experiments are also provided.
We consider the distributional fixed-point equation: $$R \stackrel{\mathcal{D}}{=} Q \vee \left( \bigvee_{i=1}^N C_i R_i \right),$$ where the $\{R_i\}$ are i.i.d.~copies of $R$, independent of the vector $(Q, N, \{C_i\})$, where $N \in \mathbb{N}$, $Q, \{C_i\} \geq 0$ and $P(Q > 0) > 0$. By setting $W = \log R$, $X_i = \log C_i$, $Y = \log Q$ it is equivalent to the high-order Lindley equation $$W \stackrel{\mathcal{D}}{=} \max\left\{ Y, \, \max_{1 \leq i \leq N} (X_i + W_i) \right\}.$$ It is known that under Kesten assumptions, $$P(W > t) \sim H e^{-αt}, \qquad t \to \infty,$$ where $α>0$ solves the Cramér-Lundberg equation $E \left[ \sum_{j=1}^N C_i ^α\right] = E\left[ \sum_{i=1}^N e^{αX_i} \right] = 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.
In this work, we adapt a Monte Carlo algorithm introduced by Broadie and Glasserman (1997) to price a $π$-option. This method is based on the simulated price tree that comes from discretization and replication of possible trajectories of the underlying asset's price. As a result this algorithm produces the lower and the upper bounds that converge to the true price with the increasing depth of the tree. Under specific parametrization, this $π$-option is related to relative maximum drawdown and can be used in the real-market environment to protect a portfolio against volatile and unexpected price drops. We also provide some numerical analysis.
This paper considers an optimal dividend distribution problem for an insurance company where the dividends are paid in a foreign currency. In the absence of dividend payments, our risk process follows a spectrally negative Lévy process. We assume that the exchange rate is described by a an exponentially Lévy process, possibly containing the same risk sources like the surplus of the insurance company under consideration. The control mechanism chooses the amount of dividend payments. The objective is to maximise the expected dividend payments received until the time of ruin and a penalty payment at the time of ruin, which is an increasing function of the size of the shortfall at ruin. A complete solution is presented to the corresponding stochastic control problem. Via the corresponding Hamilton--Jacobi--Bellman equation we find the necessary and sufficient conditions for optimality of a single dividend barrier strategy. A number of numerical examples illustrate the theoretical analysis.
This paper discusses the valuation of credit default swaps, where default is announced when the reference asset price has gone below certain level from the last record maximum, also known as the high-water mark or drawdown. We assume that the protection buyer pays premium at fixed rate when the asset price is above a pre-specified level and continuously pays whenever the price increases. This payment scheme is in favour of the buyer as she only pays the premium when the market is in good condition for the protection against financial downturn. Under this framework, we look at an embedded option which gives the issuer an opportunity to call back the contract to a new one with reduced premium payment rate and slightly lower default coverage subject to paying a certain cost. We assume that the buyer is risk neutral investor trying to maximize the expected monetary value of the option over a class of stopping time. We discuss optimal solution to the stopping problem when the source of uncertainty of the asset price is modelled by Lévy process with only downward jumps. Using recent development in excursion theory of Lévy process, the results are given explicitly in terms of scale function of the Lévy process. Furthermore, the value function of the stopping problem is shown to satisfy continuous and smooth pasting conditions regardless of regularity of the sample paths of the Lévy process. Optimality and uniqueness of the solution are established using martingale approach for drawdown process and convexity of the scale function under Esscher transform of measure. Some numerical examples are discussed to illustrate the main results.
We consider in this paper a risk reserve process where the claims and gains arrive according to two independent Poisson processes. While the gain sizes are phase-type distributed, we assume instead that the claim sizes are phase-type perturbed by a heavy-tailed component; that is, the claim size distribution is formally chosen to be phase-type with large probability $1-ε$ and heavy-tailed with small probability $ε$. We analyze the seminal Gerber-Shiu function coding the joint distribution of the time to ruin, the surplus immediately before ruin, and the deficit at ruin. We derive its value as an expansion with respect to powers of $ε$ with known coefficients and we construct approximations from the first two terms of the aforementioned series. The main idea is based on the so-called fluid embedding that allows to put the considered risk process into the framework of spectrally negative Markov-additive processes and use its fluctuation theory developed in Ivanovs and Palmowski (2012).
We revisit the optimal capital structure model with endogenous bankruptcy first studied by Leland \cite{Leland94} and Leland and Toft \cite{Leland96}. Differently from the standard case, where shareholders observe continuously the asset value and bankruptcy is executed instantaneously without delay, we assume that the information of the asset value is updated only at intervals, modeled by the jump times of an independent Poisson process. Under the spectrally negative Lévy model, we obtain the optimal bankruptcy strategy and the corresponding capital structure. A series of numerical studies are given to analyze the sensitivity of observation frequency on the optimal solutions, the optimal leverage and the credit spreads.
In this paper we provide the analysis of the limiting conditional distribution (Yaglom limit) for stochastic fluid models (SFMs), a key class of models in the theory of matrix-analytic methods. So far, transient and stationary analyses of the SFMs have been only considered in the literature. The limiting conditional distribution gives useful insights into what happens when the process has been evolving for a long time, given its busy period has not ended yet. We derive expressions for the Yaglom limit in terms of the singularity $s^*$ such that the key matrix of the SFM, ${\bfΨ}(s)$, is finite (exists) for all $s\geq s^*$ and infinite for $s<s^*$. We show the uniqueness of the Yaglom limit and illustrate the application of the theory with simple examples.