  As corporates and governments become more digital, they become vulnerable to
various forms of cyber attack. Cyber insurance products have been used as risk
management tools, yet their pricing does not reflect actual risk, including
that of multiple, catastrophic and contagious losses. For the modelling of
aggregate losses from cyber events, in this paper we introduce a bivariate
compound dynamic contagion process, where the bivariate dynamic contagion
process is a point process that includes both externally excited joint jumps,
which are distributed according to a shot noise Cox process and two separate
self-excited jumps, which are distributed according to the branching structure
of a Hawkes process with an exponential fertility rate, respectively. We
analyse the theoretical distributional properties for these processes
systematically, based on the piecewise deterministic Markov process developed
by Davis (1984) and the univariate dynamic contagion process theory developed
by Dassios and Zhao (2011). The analytic expression of the Laplace transform of
the compound process and its moments are presented, which have the potential to
be applicable to a variety of problems in credit, insurance, market and other
operational risks. As an application of this process, we provide insurance
premium calculations based on its moments. Numerical examples show that this
compound process can be used for the modelling of aggregate losses from cyber
events. We also provide the simulation algorithm for statistical analysis,
further business applications and research.
