Copula is one of the widely used techniques to describe the dependency structure between components of a system. Among all existing copulas, the family of Archimedean copulas is the popular one due to its wide range of capturing the dependency structures. In this paper, we consider the systems that are formed by dependent and identically distributed components, where the dependency structures are described by Archimedean copulas. We study some stochastic comparisons results for series, parallel, and general $r$-out-of-$n$ systems. Furthermore, we investigate whether a system of used components performs better than a used system with respect to different stochastic orders. Furthermore, some aging properties of these systems have been studied. Finally, some numerical examples are given to illustrate the proposed results.

This paper studies the optimal allocation policy of a coherent system with independent heterogeneous components and dependent subsystems, the systems are assumed to consist of two groups of components whose lifetimes follow proportional hazard (PH) or proportional reversed hazard (PRH) models. We investigate the optimal allocation strategy by finding out the number $k$ of components coming from Group A in the up-series system. First, some sufficient conditions are provided in the sense of the usual stochastic order to compare the lifetimes of two-parallel–series systems with dependent subsystems, and we obtain the hazard rate and reversed hazard rate orders when two subsystems have independent lifetimes. Second, similar results are also obtained for two-series–parallel systems under certain conditions. Finally, we generalize the corresponding results to parallel–series and series–parallel systems with multiple subsystems in the viewpoint of the minimal path and the minimal cut sets, respectively. Some numerical examples are presented to illustrate the theoretical findings.

The star-shaped ordering between probability distributions is a common way to express aging properties. A well-known criterion was proposed by Saunders and Moran [(1978). On the quantiles of the gamma and F distributions. Journal of Applied Probability 15(2): 426–432], to order families of distributions depending on one real parameter. However, the lifetime of complex systems usually depends on several parameters, especially when considering heterogeneous components. We extend the Saunders and Moran criterion characterizing the star-shaped order when the multidimensional parameter moves along a given direction. A few applications to the lifetime of complex models, namely parallel and series models assuming different individual components behavior, are discussed.

By now there is a solid theory for Polya urns. Finding the covariances is somewhat laborious. While these papers are “structural,” our purpose here is “computational.” We propose a practicable method for building the asymptotic covariance matrix in tenable balanced urn schemes, whereupon the asymptotic covariance matrix is obtained by solving a linear system of equations. We demonstrate the use of the method in growing tenable balanced irreducible schemes with a small index and in critical urns. In the critical case, the solution to the linear system of equations is explicit in terms of an eigenvector of the scheme.

In this paper, we identify some conditions to compare the largest order statistics from resilience-scale models with reduced scale parameters in the sense of mean residual life order. As an example of the established result, the exponentiated generalized gamma distribution is examined. Also, for the special case of the scale model, power-generalized Weibull and half-normal distributions are investigated.

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.

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.

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.

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.

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.

A retrial queue with classical retrial policy, where each blocked customer in the orbit retries for service, and general retrial times is modeled by a piecewise deterministic Markov process (PDMP). From the extended generator of the PDMP of the retrial queue, we derive the associated martingales. These results are used to derive the conditional expected number of customers in the orbit in the transient regime.

Innovations in medicine provide us longer and healthier life, leading lower mortality. Sooner rather than later, much greater longevity would be possible for us due to artificial intelligence advances in health care. Similarly, Advanced Driver Assistance Systems (ADAS) in highly automated vehicles may reduce or even eventually eliminate accidents by perceiving dangerous situations, which would minimize the number of accidents and lead to fewer loss claims for insurance companies. To model the survivor function capturing greater longevity as well as the number of claims reflecting less accidents in the long run, in this paper, we study a Cox process whose intensity process is piecewise-constant and decreasing. We derive its ultimate distributional properties, such as the Laplace transform of intensity integral process, the probability generating function of point process, their associated moments and cumulants, and the probability of no more claims for a given time point. In general, this simple model may be applicable in many other areas for modeling the evolution of gradually disappearing events, such as corporate defaults, dividend payments, trade arrivals, employment of a certain job type (e.g., typists) in the labor market, and release of particles. In particular, we discuss some potential applications to insurance.

This paper is concerned with developing low variance simulation estimators of probabilities related to the sum of Bernoulli random variables. It shows how to utilize an identity used in the Chen-Stein approach to bounding Poisson approximations to obtain low variance estimators. Applications and numerical examples in such areas as pattern occurrences, generalized coupon collecting, system reliability, and multivariate normals are presented. We also consider the problem of estimating the probability that a positive linear combination of Bernoulli random variables is greater than some specified value, and present a simulation estimator that is always less than the Markov inequality bound on that probability.

This paper proposes a novel stochastic volatility model with a flexible jump structure. This model allows both contemporaneous and independent arrival of jumps in return and volatility. Moreover, time-varying jump intensities are used to capture jump clustering. In the proposed framework, we provide a semi-analytical solution for the pricing problem of VIX futures and options. Through numerical experiments, we verify the accuracy of our pricing formula and explore the impact of the jump structure on the pricing of VIX derivatives. We find that the correct identification of the market jump structure is crucial for pricing VIX derivatives, and misspecified model setting can yield large errors in pricing.

In this paper, we study the credit default swap (CDS) pricing with counterparty risk in a reduced form model. The default jump intensities of the reference firm and counterparty are both assumed to follow the mean-reverting CIR processes with independent jumps respectively and a common jump. The approximate closed-form solutions of the joint survival probability density and the probability density of the first default can be obtained by using the PDE method. Then with the expressions of the probability densities, we can get the formula for the CDS price with counterparty risk in a reduced form model with a common jump. In the numerical analysis part, we find that the default of the reference asset has a greater impact on the CDS price than that of the default of counterparty after introducing the common jump process.

This paper concentrates on the fundamental concepts of entropy, information and divergence to the case where the distribution function and the respective survival function play the central role in their definition. The main aim is to provide an overview of these three categories of measures of information and their cumulative and survival counterparts. It also aims to introduce and discuss Csiszár's type cumulative and survival divergences and the analogous Fisher's type information on the basis of cumulative and survival functions.