In the usual shock models, the shocks arrive from a single source. Bozbulut and Eryilmaz [(2020). Generalized extreme shock models and their applications. Communications in Statistics – Simulation and Computation 49(1): 110–120] introduced two types of extreme shock models when the shocks arrive from one of $m\geq 1$ possible sources. In Model 1, the shocks arrive from different sources over time. In Model 2, initially, the shocks randomly come from one of $m$ sources, and shocks continue to arrive from the same source. In this paper, we prove that the lifetime of Model 1 is less than Model 2 in the usual stochastic ordering. We further show that if the inter-arrival times of shocks have increasing failure rate distributions, then the usual stochastic ordering can be generalized to the hazard rate ordering. We study the stochastic behavior of the lifetime of Model 2 with respect to the severity of shocks using the notion of majorization. We apply the new stochastic ordering results to show that the age replacement policy under Model 1 is more costly than Model 2.

In this paper, we consider a mixed dividend strategy in a dual risk model. The mixed dividend strategy is the combination of a threshold dividend and a Parisian implementation delays dividend under periodic observation. Given a series of discrete observation points, when the surplus level is larger than the predetermined bonus barrier at observation point, the Parisian implementation delays dividend is immediately carried out, and the threshold dividend is performed continuously during the delayed period. We study the Gerber-Shiu expected discounted penalty function and the expected discounted dividend payments before ruin in such a dual risk model. Numerical illustrations are given to study the influence of relevant parameters on the ruin-related quantities and the selection of the optimal dividend barrier for a given initial surplus level.

Due to the importance of generalized order statistics (GOS) in many branches of Statistics, a wide interest has been shown in investigating stochastic comparisons of GOS. In this article, we study the likelihood ratio ordering of $p$-spacings of GOS, establishing some flexible and applicable results. We also settle certain unresolved related problems by providing some useful lemmas. Since we do not impose restrictions on the model parameters (as previous studies did), our findings yield new results for comparison of various useful models of ordered random variables including order statistics, sequential order statistics, $k$-record values, Pfeifer's record values, and progressive Type-II censored order statistics with arbitrary censoring plans. Some results on preservation of logconvexity properties among spacings are provided as well.

Sequences of non-decreasing (non-increasing) lower (upper) bounds for the renewal-type equation as well as for the renewal function which are improvements of the famous corresponding bounds of Marshal [(1973). Linear bounds on the renewal function. SIAM Journal on Applied Mathematics 24(2): 245–250] are given. Also, sequences such bounds converging to the ordinary renewal function are obtained for several reliability classes of the lifetime distributions of the inter-arrival times, which are refinements of all of the existing known corresponding bounds. For the first time, a lower bound for the renewal function with DMRL lifetimes is given. Finally, sequences of such improved bounds are given for the ordinary renewal density as well as for the right-tail of the distribution of the forward recurrence time.

We consider the problem of group testing (pooled testing), first introduced by Dorfman. For nonadaptive testing strategies, we refer to a nondefective item as “intruding” if it only appears in positive tests. Such items cause misclassification errors in the well-known COMP algorithm and can make other algorithms produce an error. It is therefore of interest to understand the distribution of the number of intruding items. We show that, under Bernoulli matrix designs, this distribution is well approximated in a variety of senses by a negative binomial distribution, allowing us to understand the performance of the two-stage conservative group testing algorithm of Aldridge.

In this paper, we extend the optimal dividend and capital injection problem with affine penalty at ruin in (Xu, R. & Woo, J.K. (2020). Insurance: Mathematics and Economics 92: 1–16) to the case with singular dividend payments. The asymptotic relationships between our value function to the one with bounded dividend density are studied, which also help to verify that our value function is a viscosity solution to the associated Hamilton–Jacob–Bellman Quasi-Variational Inequality (HJBQVI). We also show that the value function is the smallest viscosity supersolution within certain functional class. A modified comparison principle is proved to guarantee the uniqueness of the value function as the viscosity solution within the same functional class. Finally, a band-type dividend and capital injection strategy is constructed based on four crucial sets; and the optimality of such band-type strategy is proved by using fixed point argument. Numerical examples of the optimal band-type strategies are provided at the end when the claim size follows exponential and gamma distribution, respectively.

A company with $n$ geographically widely dispersed sites seeks insurance that pays off if $m$ out of the $n$ sites experience rarely occurring catastrophes (e.g., earthquakes) during a year. This study describes an adaptive dynamic strategy that enables an insurance company to offer the policy with smaller loss probability than more conventional static policies induce, but at a comparable or lower premium. The strategy accomplishes this by periodically purchasing reinsurance on individual sites. Exploiting rarity, the policy induces zero loss with probability one if no more than one quake occurs during any review interval. The policy also may induce a profit if $m$ or more quakes occur in an interval if no quakes have occurred in previous intervals. The study also examines the benefit of more than one reinsurance policy per active site. The study relies on expected utility to determine indifference premiums and derives an upper bound on loss probability independent of premium.

This paper studies the open-loop equilibrium strategies for a class of non-zero-sum reinsurance–investment stochastic differential games between two insurers with a state-dependent mean expectation in the incomplete market. Both insurers are able to purchase proportional reinsurance contracts and invest their wealth in a risk-free asset and a risky asset whose price is modeled by a general stochastic volatility model. The surplus processes of two insurers are driven by two standard Brownian motions. The objective for each insurer is to find the equilibrium investment and reinsurance strategies to balance the expected return and variance of relative terminal wealth. Incorporating the forward backward stochastic differential equations (FBSDEs), we derive the sufficient conditions and obtain the general solutions of equilibrium controls for two insurers. Furthermore, we apply our theoretical results to two special stochastic volatility models (Hull–White model and Heston model). Numerical examples are also provided to illustrate our results.

In the literature, some stochastic orders have been extended to the higher orders in different scenarios. In this paper, inspired by interesting properties of the excess wealth order and its wide range application particularly in comparing the tail variability of risks, we consider the second-order excess wealth order and study its main properties. We obtain two results characterizing the proposed order. We also investigate its relationship with other well-known variability orders and criteria to compare risks. An application of the results in comparing the epoch times of two nonhomogeneous poisson processes is also given.

As the queue becomes exhausted, different maintenance tasks can be performed according to the fatigue load and wear degree of the service equipment. At the same time, considering the customer's sensitivity to time delay, the service facility will not completely remain inactive during the maintenance period. To describe this objectively existing phenomenon arising in the waiting line system, we consider a hyper-exponential working vacation queue with a batch renewal arrival process. Through the calculation of the well-structured roots of the associated characteristic equation, the shift operator method in the theory of difference equations and the supplementary variable technique for stochastic modeling plays a central role in the queue-length distribution analysis. Comparison with other ways to analyze queueing models, the advantage of our approach is that we can avoid deriving the complex transition probability matrix of the queue-length process embedded at input points. The feasibility of this approach is verified by extensive numerical examples.

In this article, we provide a comprehensive analyses of two continuous review lost sales inventory system based on different replenishment policies, namely $(s,S)$ and $(s,Q)$. We assume that the arrival times of demands form a Poisson process and that the demand sizes have i.i.d. exponential distribution. We assume that the items in stock may obsolete after an exponential time. The lead time for replenishment is exponential. We also assume that the excess demands and the demands that occurred during stock out periods are lost. Using the system point method of level crossing and integral equation method, we derive the steady-state probability distribution of inventory level explicitly. After deriving some system performance measures, we computed the total expected cost rate. We also provide numerical examples of sensitivity analyses involving different parameters and costs. As a result of our numerical analysis, we provide several insights on the optimal $(s,S)$ and $(s,Q)$ policies for inventory systems of obsolescence items with positive lead times. The better policy for maintaining inventory can be quantified numerically.

In this paper, the multi-state survival signature is first redefined for multi-state coherent or mixed systems with independent and identically distributed (i.i.d.) multi-state components. With the assumption of independence of component lifetimes at different state levels, transformation formulas of multi-state survival signatures of different sizes are established through the use of equivalent systems and a generalized triangle rule for order statistics from several independent and non-identical distributions. The results obtained facilitate stochastic comparisons of multi-state coherent or mixed systems with different numbers of i.i.d. multi-state components. Specific examples are finally presented to illustrate the transformation formulas established here, and also their use in comparing systems of different sizes.

We introduce the generalized join the shortest queue model with retrials and two infinite capacity orbit queues. Three independent Poisson streams of jobs, namely a smart, and two dedicated streams, flow into a single-server system, which can hold at most one job. Arriving jobs that find the server occupied are routed to the orbits as follows: Blocked jobs from the smart stream are routed to the shortest orbit queue, and in case of a tie, they choose an orbit randomly. Blocked jobs from the dedicated streams are routed directly to their orbits. Orbiting jobs retry to connect with the server at different retrial rates, i.e., heterogeneous orbit queues. Applications of such a system are found in the modeling of wireless cooperative networks. We are interested in the asymptotic behavior of the stationary distribution of this model, provided that the system is stable. More precisely, we investigate the conditions under which the tail asymptotic of the minimum orbit queue length is exactly geometric. Moreover, we apply a heuristic asymptotic approach to obtain approximations of the steady-state joint orbit queue-length distribution. Useful numerical examples are presented and shown that the results obtained through the asymptotic analysis and the heuristic approach are agreed.

This paper studies a Pareto-optimal reinsurance problem when the contract is subject to default of the reinsurer. We assume that the reinsurer can invest a share of its wealth in a risky asset and default occurs when the reinsurer's end-of-period wealth is insufficient to cover the indemnity. We show that without the solvency regulation, the optimal indemnity function is of excess-of-loss form, regardless of the investment decision. Under the solvency regulation constraint, by assuming the investment decision remains unchanged, the optimal indemnity function is characterized element-wisely. Partial results are derived when both the indemnity function and investment decision are impacted by the solvency regulation. Numerical examples are provided to illustrate the implications of our results and the sensitivity of solution to the model parameters.

Let $V_{(r,n,\tilde {m}_n,k)}^{(p)}$ and $W_{(r,n,\tilde {m}_n,k)}^{(p)}$ be the $p$-spacings of generalized order statistics based on absolutely continuous distribution functions $F$ and $G$, respectively. Imposing some conditions on $F$ and $G$ and assuming that $m_1=\cdots =m_{n-1}$, Hu and Zhuang (2006. Stochastic orderings between p-spacings of generalized order statistics from two samples. Probability in the Engineering and Informational Sciences 20: 475) established $V_{(r,n,\tilde {m}_n,k)}^{(p)} \leq _{{\rm hr}} W_{(r,n,\tilde {m}_n,k)}^{(p)}$ for $p=1$ and left the case $p\geq 2$ as an open problem. In this article, we not only resolve it but also give the result for unequal $m_i$'s. It is worth mentioning that this problem has not been proved even for ordinary order statistics so far.

This paper examines the value-at-risk (VaR) implications of mean-variance hedging. We derive an equivalence between the VaR-based hedge and the mean-variance hedging. This method transfers the investor's subjective risk-aversion coefficient into the estimated VaR measure. As a result, we characterize the collapse probability bounds under which the VaR-based hedge could be insignificantly different from the minimum-variance hedge in the presence of estimation risk. The results indicate that the squared information ratio of futures returns is the primary factor determining the difference between the minimum-variance and VaR-based hedges.

Zhu and He [(2018). A new closed-form formula for pricing European options under a skew Brownian motion. The European Journal of Finance 24(12): 1063–1074] provided an innovative closed-form solution by replacing the standard Brownian motion in the Black–Scholes framework using a particular skew Brownian motion. Their formula involves numerically integrating the product of the Guassian density and corresponding distribution function. Being different from their pricing formula, we derive a much simpler formula that only involves the Gaussian distribution function and Owen's $T$ function.

We obtain here sufficient conditions for increasing concave order and location independent more riskier order of lower record values based on stochastic comparisons of minimum order statistics. We further discuss stochastic orderings of lower record spacings. In particular, we show that increasing convex order of adjacent spacings between minimum order statistics is a sufficient condition for increasing convex order of adjacent spacings of their lower records.

Layer reinsurance treaty is a common form obtained in the problem of optimal reinsurance design. In this paper, we study allocations of policy limits in layer reinsurance treaties with dependent risks. We investigate the effects of orderings and heterogeneity among policy limits on the expected utility functions of the terminal wealth from the viewpoint of risk-averse insurers faced with right tail weakly stochastic arrangement increasing losses. Orderings on optimal allocations are presented for normal layer reinsurance contracts under certain conditions. Parallel studies are also conducted for randomized layer reinsurance contracts. As a special case, the worst allocations of policy limits are also identified when the exact dependence structure among the losses is unknown. Numerical examples are presented to shed light on the theoretical findings.