Search

Dependence among cyber risks has been an essential and challenging component of risk management. The current study characterizes cyber dependence from both qualitative and quantitative perspectives based on L-hop propagation model. From the qualitative side, it is shown that cyber risks always possess positive association based on the proposed risk propagation model. From the quantitative side, an explicit formula for computing the fundamental dependence measure of covariance is provided for an arbitrary network. In particular, we study the impacts of factors—especially external and internal compromise probabilities, propagation depth, and network topologies—on dependence among cyber risks. We conclude by presenting some examples and applications.

Modeling cyber threats, such as the computer malicious software (malware) propagation dynamics in cyberspace, is an important research problem because models can deepen our understanding of dynamical cyber threats. In this paper, we study the statistical modeling of the macro-level evolution of dynamical cyber attacks. Specifically, we propose a Bayesian structural time series approach for modeling the computer malware propagation dynamics in cyberspace. Our model not only possesses the parsimony property (i.e. using few model parameters) but also can provide the predictive distribution of the dynamics by accommodating uncertainty. Our simulation study shows that the proposed model can fit and predict the computer malware propagation dynamics accurately, without requiring to know the information about the underlying attack-defense interaction mechanism and the underlying network topology. We use the model to study the propagation of two particular kinds of computer malware, namely the Conficker and Code Red worms, and show that our model has very satisfactory fitting and prediction accuracies.

Data breaches cause millions of dollars in financial losses each year. The insurance industry has been exploring the ways to transfer such extreme risk. In this work, we investigate data breach cat...

Malicious hacking data breaches cause millions of dollars in financial losses each year, and more companies are seeking cyber insurance coverage. The lack of suitable statistical approaches for sco...

Cybersecurity risk has attracted considerable attention in recent decades. However, the modeling of cybersecurity risk is still in its infancy, mainly because of its unique characteristics. In this...

Modeling cybersecurity risks is an important, yet challenging, problem. In this paper, we initiate the study of modeling multivariate cybersecurity risks. We develop the first statistical a...

Cyber attacks have become a problem that is threatening the economy, human privacy, and even national security. Before we can adequately address the problem, we need to have a crystal clear understanding about cyber attacks from various perspectives. This is a challenge because the Internet is a large-scale complex system with humans in the loop. In this paper, we investigate a particular perspective of the problem, namely the extreme value phenomenon that is exhibited by cyber attack rates, which are the numbers of attacks against a system of interest per time unit. It is important to explore this perspective because understanding the statistical properties of extreme cyber attack rates will pave the way for cost-effective, if not optimal, allocation of resources in real-life cyber defense operations. Specifically, we propose modeling and predicting extreme cyber attack rates via marked point processes, while using the Value-at-Risk as a natural measure of intense cyber attacks. The point process...

In this paper, we propose a novel method for constructing upper bounds of the quasi-stationary distribution of SIS processes. Using this method, we obtain an upper bound that is better than the state-of-the-art upper bound. Moreover, we prove that the fixed point map Φ [7] actually preserves the equilibrium reversed hazard rate order under a certain condition. This allows us to further improve the upper bound. Some numerical results are presented to illustrate the results.

Summary Discrete power-law distributions have significant consequences for understanding many phenomena in practice, and have attracted much attention in recent decades. However, in many practical applications, there exists a natural upper bound for the probability tail. In this paper, we develop maximum likelihood estimates for truncated discrete power-law distributions based on the upper order statistics, and large sample properties are mentioned as well. Monte Carlo simulation is carried out to examine the finite sample performance of the estimates. Applications in real cyber attack data and peak gamma-ray intensity of solar flares are highlighted.

This paper studies capital allocation problems with the aggregate risk exceeding a certain threshold. We propose a novel capital allocation rule based on the Tail Mean–Variance principle. General formulas for the optimal capital allocations are proposed. Explicit formulas for optimal capital allocations are derived for multivariate elliptical distributions. Moreover, we give asymptotic allocation formulas for multivariate regular variation variables. Various numerical examples are given to illustrate the results, and real insurance data is discussed as well.

Professors Balakrishnan and Zhao have written an excellent survey on the recent developments of stochastic comparisons of order statistics, which cover almost every aspect of ordering properties of order statistics from both continuous and discrete heterogeneous populations. My discussion will be limited to the skewness of order statistics and order statistics from heterogeneous populations with different shape parameters.

In this paper, the sample range from a heterogeneous exponential sample is shown to be larger than that from a homogeneous exponential sample in the sense of the star ordering. Then, by using this result, some equivalent characterizations of stochastic comparisons of sample ranges with respect to various stochastic orders are established. In this process, two open problems mentioned in Mao and Hu (2010) [16] are solved. The main results established here extend and strengthen several known results in the literature including those of Khaledi and Kochar (2000) [8] , Zhao and Li (2009) [22] and Genest et al. (2009) [7] .

In this paper, a new sufficient condition for comparing linear combinations of independent gamma random variables according to star ordering is given. This unifies some of the newly proved results on this problem. Equivalent characterizations between various stochastic orders are established by utilizing the new condition. The main results in this paper generalize and unify several results in the literature including those of Amiri, Khaledi, and Samaniego [2], Zhao [18], and Kochar and Xu [9].

This paper studies capital allocation problems using a general loss function. Stochastic comparisons are conducted for general loss functions in several scenarios: independent and identically distributed risks; independent but non-identically distributed risks; comonotonic risks. Applications in optimal capital allocations and policy limits allocations are discussed as well.

In this paper we provide some sufficient conditions to stochastically compare linear combinations of independent random variables. The main results extend those given in Proschan (1965), Ma (1998), Zhao et al. (2011), and Yu (2011). In particular, we propose a new sufficient condition to compare the peakedness of linear combinations of independent random variables which may have heavy-tailed properties.

In this article, we study the order statistics from heterogenous negative binomial random variables. Sufficient conditions are provided for comparing the extreme order statistics according to the usual stochastic order. For the special case of geometric distribution, a sufficient condition is established for comparing order statistics in the sense of multivariate stochastic order. Applications in the Poisson–Gamma shock model and redundant systems have been described as well.

Two sufficient conditions for comparing convolutions of heterogeneous gamma random variables in terms of star order are established. It is further shown that if the scale parameters of heterogeneous gamma random variables are more dispersed in the sense of majorization, then the convolutions are more dispersed according to the right spread order, which generalizes and strengthens the results in Diaconis and Perlman (1987), Kochar and Xu (2010) and Zhao and Balakrishnan (2009).

A sufficient condition for comparing convolutions of heterogeneous exponential random variables in terms of right spread order is established. As a consequence, it is shown that a convolution of heterogeneous independent exponential random variables is more skewed than that of homogeneous exponential random variables in the sense of NBUE order. This gives a new insight into the distribution theory of convolutions of independent random variables. A sufficient condition is also derived for comparing such convolutions in terms of Lorenz order.

In this article, mixture representations of survival functions of residual lifetimes of k-out-of-n systems are obtained when the components are independent but not necessarily identically distributed. Then we stochastically compare the residual lifetimes of k-out-of-n systems in one- and two-sample problems. In particular, the results extend some results in Li and Zhao [14], Khaledi and Shaked [13], Sadegh [17], Gurler and Bairamov [7] and Navarro, Balakrishnan, and Samaniego [16]. Applications in the proportional hazard rates model are presented as well.

Let R"n be the range of a random sample X"1,...,X"n of exponential random variables with hazard rate @l. Let S"n be the range of another collection Y"1,...,Y"n of mutually independent exponential random variables with hazard rates @l"1,...,@l"n whose average is @l. Finally, let r and s denote the reversed hazard rates of R"n and S"n, respectively. It is shown here that the mapping [email protected]?s(t)/r(t) is increasing on (0,~) and that as a result, R"n=X"("n")-X"("1") is smaller than S"n=Y"("n")-Y"("1") in the likelihood ratio ordering as well as in the dispersive ordering. As a further consequence of this fact, X"("n") is seen to be more stochastically increasing in X"("1") than Y"("n") is in Y"("1"). In other words, the pair (X"("1"),X"("n")) is more dependent than the pair (Y"("1"),Y"("n")) in the monotone regression dependence ordering. The latter finding extends readily to the more general context where X"1,...,X"n form a random sample from a continuous distribution while Y"1,...,Y"n are mutually independent lifetimes with proportional hazard rates.

A parallel system with heterogeneous exponential component lifetimes is shown to be more skewed (according to the convex transform order) than the system with independent and identically distributed exponential components. As a consequence, equivalent conditions for comparing the variabilities of the largest order statistics from heterogeneous and homogeneous exponential samples in the sense of the dispersive order and the right-spread order are established. A sufficient condition is also given for the proportional hazard rate model.

A parallel system with heterogeneous exponential component lifetimes is shown to be more skewed (according to the convex transform order) than the system with independent and identically distributed exponential components. As a consequence, equivalent conditions for comparing the variabilities of the largest order statistics from heterogeneous and homogeneous exponential samples in the sense of the dispersive order and the right-spread order are established. A sufficient condition is also given for the proportional hazard rate model.

In this paper, we introduce a new copula-based dependence order to compare the relative degree of dependence between two pairs of random variables. Relationship of the new order to the existing dependence orders is investigated. In particular, the new ordering is stronger than the partial ordering, more monotone regression dependence as developed by Avérous et al. [J. Avérous, C. Genest, S.C. Kochar, On dependence structure of order statistics, Journal of Multivariate Analysis 94 (2005) 159–171]. Applications of this partial order to order statistics, k-record values and frailty models are given.

This note investigates the negative dependence in frailty models. First, we show that the frailty variable and the overall population variable are negatively likelihood ratio dependent and derive an upper bound for the survival function of the population with higher frailty. Secondly, we prove that the DFR property and the logconvex hazard rate of the baseline variable imply the DLR property of the population variable. Finally, we further prove that the likelihood ratio order, hazard rate order and reversed hazard rate order between two frailty variables imply the likelihood ratio order, reversed hazard rate order, and hazard rate order between the corresponding overall population variables, respectively.

This paper investigates the stochastic comparison of spacings of record values. The likelihood ratio order, the hazard rate order and the mean residual life order between spacings of record values from one sample sequence are developed. Furthermore, the increasing convex order and the mean residual life order between simple spacings of record values from two sample sequences are built as well.

Let X1, … , Xn be independent random variables with Xi having survival function Fλi, i = 1, … , n, and let Y1, … ,Yn be a random sample with common population survival distribution F, where c = ∑i=1nλi/n. Let Xn:n and Yn:n denote the lifetimes of the parallel systems consisting of these components, respectively. It is shown that Xn:n is greater than Yn:n in terms of likelihood ratio order. It is also proved that the sample range Xn:n − X1:n is larger than Yn:n − Y1:n according to reverse hazard rate ordering. These two results strengthen and generalize the results in Dykstra, Kochar, and Rojo [6] and Kochar and Rojo [11], respectively.

This paper deals with preservation of the reversed hazard rate order between equilibrium random variables under formations of some reliability structures. We further investigate a new aging notion based upon the reversed hazard rate order between a random life and its equilibrium version. A nonparametric method is developed to test the exponentiality against such a strict aging property, some numerical results are presented as well.

For two non-negative random variables X and Y, it is shown that if X is smaller than Y in the sense of the likelihood ratio order and either X or Y is of decreasing likelihood ratio, then the m-spacing V k : n ( m ) of an X-sample is smaller than W k : n ( m ) , the m-spacing of a Y-sample, in terms of the likelihood ratio order. It is also proved that a similar result for the upshifted likelihood ratio order is valid. Finally, it is proved that the decreasing failure rate in average property of a population X is preserved by sample spacing V k : n ( 1 ) .

We investigate some new properties of mean inactivity time (MIT) order and increasing MIT (IMIT) class of life distributions. The preservation property of MIT order under increasing and concave transformations, reversed preservation properties of MIT order, and IMIT class of life distributions under the taking of maximum are developed. Based on the residual life at a random time and the excess lifetime in a renewal process, stochastic comparisons of both IMIT and decreasing mean residual life distributions are conducted as well.

In our article [3] we have found a gap in the middle of the proof of Theorem 3.2. Therefore, we do not know whether Theorem 3.2 is true for the reverse hazard rate order. However, we could prove the following weaker result for the stochastic order.