
In this article, we introduce the fundamentals of the theory of demicontractive mappings in metric spaces and expose the key concepts and tools for building a constructive approach to approximating the fixed points of demicontractive mappings in this setting. By using an appropriate geodesic averaged perturbation technique, we obtained strong convergence and Δ-convergence theorems for a Krasnoselskij-Mann type iterative algorithm to approximate the fixed points of quasi-nonexpansive mappings within the framework of CAT(0) spaces. The main results obtained in this nonlinear setting are natural extensions of the classical results from linear settings (Hilbert and Banach spaces) for both quasi-nonexpansive mappings and demicontractive mappings. We applied our results to solving an equilibrium problem in CAT(0) spaces and showed how we can approximate the equilibrium points by using our fixed point results. In this context we also provided a numerical example in the case of a demicontractive mapping that is not a quasi-nonexpansive mapping and highlighted the convergence pattern of the algorithm in
Citation: Sani Salisu, Vasile Berinde, Songpon Sriwongsa, Poom Kumam. Approximating fixed points of demicontractive mappings in metric spaces by geodesic averaged perturbation techniques[J]. AIMS Mathematics, 2023, 8(12): 28582-28600. doi: 10.3934/math.20231463
[1] | Jin Liang, Hong-Ming Yin, Xinfu Chen, Yuan Wu . On a Corporate Bond Pricing Model with Credit Rating Migration Risksand Stochastic Interest Rate. Quantitative Finance and Economics, 2017, 1(3): 300-319. doi: 10.3934/QFE.2017.3.300 |
[2] | Sang Phu Nguyen, Toan Luu Duc Huynh . Portfolio optimization from a Copulas-GJR-GARCH-EVT-CVAR model: Empirical evidence from ASEAN stock indexes. Quantitative Finance and Economics, 2019, 3(3): 562-585. doi: 10.3934/QFE.2019.3.562 |
[3] | Lianzhang Bao, Guangliang Zhao, Zhuo Jin . A new equilibrium trading model with asymmetric information. Quantitative Finance and Economics, 2018, 2(1): 217-229. doi: 10.3934/QFE.2018.1.217 |
[4] | Manuela Larguinho, José Carlos Dias, Carlos A. Braumann . Pricing and hedging bond options and sinking-fund bonds under the CIR model. Quantitative Finance and Economics, 2022, 6(1): 1-34. doi: 10.3934/QFE.2022001 |
[5] | Gabriel Frahm . How often is the financial market going to collapse?. Quantitative Finance and Economics, 2018, 2(3): 590-614. doi: 10.3934/QFE.2018.3.590 |
[6] | Joseph Atta-Mensah . Commodity-linked bonds as an innovative financing instrument for African countries to build back better. Quantitative Finance and Economics, 2021, 5(3): 516-541. doi: 10.3934/QFE.2021023 |
[7] | Richard J. Cebula, Don Capener . The impact of federal income tax rate cuts on the municipal bond market in the U.S.: A brief exploratory empirical note. Quantitative Finance and Economics, 2018, 2(2): 407-412. doi: 10.3934/QFE.2018.2.407 |
[8] | Ke Liu, Changqing Luo, Zhao Li . Investigating the risk spillover from crude oil market to BRICS stock markets based on Copula-POT-CoVaR models. Quantitative Finance and Economics, 2019, 3(4): 754-771. doi: 10.3934/QFE.2019.4.754 |
[9] | Jie Xiong, Geng Deng, Xindong Wang . Extension of SABR Libor Market Model to handle negative interest rates. Quantitative Finance and Economics, 2020, 4(1): 148-171. doi: 10.3934/QFE.2020007 |
[10] | Kazuo Sano . New concept for the value function of prospect theory. Quantitative Finance and Economics, 2024, 8(4): 733-756. doi: 10.3934/QFE.2024028 |
In this article, we introduce the fundamentals of the theory of demicontractive mappings in metric spaces and expose the key concepts and tools for building a constructive approach to approximating the fixed points of demicontractive mappings in this setting. By using an appropriate geodesic averaged perturbation technique, we obtained strong convergence and Δ-convergence theorems for a Krasnoselskij-Mann type iterative algorithm to approximate the fixed points of quasi-nonexpansive mappings within the framework of CAT(0) spaces. The main results obtained in this nonlinear setting are natural extensions of the classical results from linear settings (Hilbert and Banach spaces) for both quasi-nonexpansive mappings and demicontractive mappings. We applied our results to solving an equilibrium problem in CAT(0) spaces and showed how we can approximate the equilibrium points by using our fixed point results. In this context we also provided a numerical example in the case of a demicontractive mapping that is not a quasi-nonexpansive mapping and highlighted the convergence pattern of the algorithm in
The year 2020 will be remembered for the global financial crisis triggered by COVID-19. Against the backdrop of disruption and upheaval caused by the pandemic, millions of people also experienced severe natural and man-made catastrophes, such as bush fires in Australia, derecho storms, Hurricane Laura and a Beirut explosion. In inflation-adjusted terms, global economic losses from natural and man-made catastrophes were USD 202 billion in 2020. The insurance industry has played an essential role, which covered USD 89 billion of the absolute economic costs incurred in 2020 on account of disaster events. Notwithstanding, the global catastrophe protection gap was around USD 113 billion (Swiss Re, 2020). Insurance-linked securities (ILSs), which connect both the capital market and the insurance market, provide an effective channel for the insurance industry to raise funds from the capital market. Catastrophe bonds (CAT bonds) are the most active ILS to date; they not only enhance the catastrophe underwriting capacity of the insurance industry, but also enrich the investment portfolios for investors.
The successful issuance of CAT bonds is inseparable from reasonable pricing. The default of the CAT bonds' principal or remaining interest payment is triggered if a catastrophic event exceeds a certain pre-suppositional condition (Smack, 2016). Clearly, defining the default trigger event plays an important role in pricing CAT bonds. In general, there are three trigger types: indemnity triggers, index triggers and hybrid triggers (Ma and Ma, 2013). The indicator of indemnity triggers is the level of actual economic loss suffered by the sponsor. There are three broad types of index triggers: industry loss indices, modeled loss indices and parametric indices. Industry loss index-triggered CAT bonds take the estimated industry-wide losses from a catastrophic event as the trigger indicator. Published data for industry loss index triggers are obtained from an index maintained by an independent third party, such as the Sigma index and Property Claim Service (PCS) index. The modeled loss index is calculated using the model provided by a catastrophe modeling agency, such as RMS and EQECAT (Ma and Ma, 2013). A parameter trigger mechanism refers to taking the physical parameters of catastrophic events as trigger indicators, such as the earthquake magnitude, typhoon wind speed and flood flow level. The moral hazard faced by bond investors and the basis risk borne by bond sponsors are two important factors to be considered in determining the trigger type for CAT bonds. When an indemnity trigger is used to insure the sponsor's actual loss, it results in the highest degree of moral hazard. The application of the three types of index triggers can eliminate the moral hazard and enhance the investment attraction for investors, but the basic risk is inevitable. It can be seen that there are obvious deficiencies in CAT bonds with a single trigger type, no matter whether it is an indemnity trigger or index trigger. A hybrid trigger mechanism refers to the simultaneous application of two or more types of trigger indicators in the contract, which can be a hybrid of different index triggers, or a hybrid of index triggers and indemnity triggers (Bouriaux and MacMinn, 2009). In practice, it is more advantageous to combine indemnity triggers and index triggers, because this combination can achieve the purpose of balancing moral hazard and basic risk.
In the last two decades, several research studies have been conducted on CAT bonds. Mainly, these studies have focused on the pricing of indemnity-triggered CAT bonds (e.g., Cox and Pedersen, 2000; Lee and Yu, 2002; Nowak and Romaniuk, 2013; Chen et al., 2013; Romaniuk, 2017; Zhang and Tsai, 2018; Kurniawan et al., 2021) and index-triggered CAT bonds (e.g., Zimbidis et al., 2007; Ma and Ma, 2013; Shao et al., 2015; Ma et al., 2017; Mousavi et al., 2019). Nevertheless, studies on the pricing of hybrid-triggered CAT bonds are not very numerous. Taking terrorism risk CAT bonds as an example, Woo (2004) proposed the concept of bonds triggered by multiple events. Reshetar (2008) developed a framework for the pricing of multiple-event CAT bonds that provides a theoretical basis for the construction of the hybrid-triggered CAT bonds pricing model. In this study, the payoffs on the bond are linked to catastrophic property losses and daily deaths, in which the attachment point for property losses is given as a threshold total annual loss, and the attachment point for death is defined as a threshold number of deaths per day (Reshetar, 2008). It should be noted that most insurance companies expect to be supplemented with capital in a timely manner after a single disaster with extreme losses to avoid bankruptcy due to substantial claims. Therefore, this study takes the property loss and magnitude of a single earthquake disaster as the hybrid trigger indicators.
CAT bonds can be considered as zero-beta assets, and the structure of a bond's cash flow depends on the occurrence of the predefined catastrophe and the attachment point of trigger indicators (e.g., Litzenberger et al., 1996; Cummins and Weiss, 2009; Tao, 2011). Therefore, the accurate assessment of catastrophe risk is a crucial step. Catastrophe risks have obvious heavy tail characteristics, and the pricing of CAT bonds is mainly related to the extremes of the tail. Extreme value theory (EVT) is an obbligato tool in the study of catastrophe; it can effectively describe the tail characteristics of the trigger indicators' distribution. It is applied widely in various fields, such as earth science, meteorology and hydrology (Ma et al., 2021). By using the PCS loss index data from 2001 to 2010, Ma et al. (2017) adopted the peaks-over-threshold (POT) model in EVT to characterize the tail characteristics of catastrophic loss distribution when evaluating zero-coupon CAT bonds. Chao and Zou (2018) evaluated the distributions of flood catastrophe losses and deaths according to EVT by using the global flood data since 1985. Based on the data on global drought catastrophe losses from 1900 to 2018, Deng et al. (2020) studied the pricing of drought CAT bonds by using EVT and high quantile estimation.
When multiple variables related to a catastrophic event are used as trigger indicators, the correlation between them should be considered. The copula modeling is a more general method that can be used to describe the dependence of a set of random variables (Chebbi and Hedhli, 2020). Frees and Valdez (1998) were among the first to use it in insurance applications. Then, the literature on the applications of copulas in catastrophe insurance has been developing rapidly in recent years. Bokusheva (2014) presented a copula-based approach for rating weather index insurance designed to provide coverage for yield losses due to extreme weather events; their results indicated that the application of the copula approach might improve the performance of weather index insurance. Xu et al. (2013) built a drought disaster evaluation model based on the copula-EVT model; it reflected the extreme process about drought duration and drought intensity. Liu et al. (2019) modeled the direct economical losses and tolls of earthquakes by using a copula-mixed distribution model. Chao (2021) employed the copula function to describe the dependence of damage areas and deaths due to flood catastrophe; they developed a framework for valuing multirisk catastrophe reinsurance contracts.
Earthquake disaster is the most severe of all disasters, as it has caused considerable losses to human beings. In this paper, we take the earthquake disaster as the insurance subject and construct the pricing model of multi-year hybrid-triggered CAT bonds under stochastic interest rates. According to the characteristics of a earthquake disaster, we choose the economic loss and magnitude of a single earthquake disaster as trigger indicators, and their marginal distributions are depicted according to EVT. By introducing an Archimedean copula function and its survival function, the valuation formula for CAT bonds' price is derived. Furthermore, the impacts of the maturity, trigger level, dependence of trigger indicators, market interest rate and coupon rate on CAT bond pricing are analyzed.
Compared with previous research, the main contributions of this paper to the field of CAT bond pricing are twofold. First, a framework is developed for valuing CAT bonds based on the hybrid trigger mechanism; it combines the indemnity trigger and index trigger. The application of the hybrid trigger mechanism can effectively enhance the investment attractiveness and market competitiveness of CAT bonds. However, compared with indemnity-triggered or index-triggered CAT bonds, the payment structure of hybrid-triggered CAT bonds is more complicated, and there are relatively few related studies. Hence, this study aims to develop a model to value hybrid-triggered CAT bonds. Second, a multi-year hybrid-triggered CAT bond pricing formula is derived using the copula-EVT model. The general statistical distributions are not suitable for fitting the heavy-tailed data since the catastrophe risk data contains some extreme values. And, there may be a certain degree of correlation between the related variables of the same catastrophic event. Therefore, the goal was to combine the EVT and copula function to analyze the distribution characteristics and the dependence of the catastrophe variables.
The rest of the paper is organized as follows. The hybrid-triggered CAT bond pricing formula is introduced and proved in Section 2. Section 3 shows the parameters estimation of the pricing model. Section 4 presents the numerical analysis. Finally, we end our work with a conclusion in Section 5.
We consider a coupon-paying CAT bond with maturity T years, and choose the economy loss and magnitude as trigger indicators. The principal is F and coupon payments Ct are paid annually at the times t=1,2,...,T. The coupon rate is denoted by R. {Let (Ω,F,P) define a probability space, where Ω is the sample space, F is a σ-algebra of subsets of Ω and P is the probability measure.} Let {Nt,t=1,2,...,T} denote the number of earthquake disasters that occurred in the interval [t−1,t], which is a sequence of independent and identically distributed (i.i.d.) random variables that obeys a Poisson distribution with an intensity λ>0. The occurrence time of the i-th earthquake disaster is denoted as di, i=1,2,...,m (m=∑Tt=1Nt). Let {Xd,d=d1,d2,...,dm} denote the economy loss caused by the earthquake disaster that occurred at the time d, and let {Yd,d=d1,d2,...,dm} denote the magnitude of the earthquake disaster that occurred at the time d. {Xd,d=d1,d2,...,dm} and {Yd,d=d1,d2,...,dm} are two sequences of i.i.d. positive random variables with cumulative distribution functions FX and FY, respectively. Both of them are independent of {Nt,t=1,2,...,T}.
In the event of a catastrophe, the coupons and principal are at risk. When one of the indicators is triggered, the current and future coupons will no longer be paid. When both indicators are triggered simultaneously, the principal at maturity will no longer be returned to investors.
We define the stopping time τC as the first time when the economy loss exceeds the trigger level TrX or the magnitude exceeds the trigger level TrY:
τC=min{d=d1,d2,...,dm|Xd>TrX or Yd>TrY}, | (1) |
Then, the payoff structure of coupons can be described as
Ct=FRI{τC>t},t=1,2,...,T, | (2) |
where I{⋅} is an indicator function.
Next, we define the stopping time τF as the first time when the economy loss and magnitude exceed their trigger levels simultaneously:
τF=min{d=d1,d2,...,dm|Xd>TrX and Yd>TrY}. | (3) |
Then, the payoff structure of the principal is given by the following:
FT=FI{τF>T}. | (4) |
Combining Equations (2) and (4), the cash flows to the bondholder can be described as
CF(t)={FRI{τC>t},t=1,2,...,T−1,FRI{τC>t}+FI{τF>t},t=T. | (5) |
There are two principal kinds of models for extreme values, namely, the block maxima (BM) model and POT model. The focus of the former is on modeling the BM of the variables with the generalized extreme value distribution (Acero et al., 2018). However, it is only interested in the behavior of the sample maximum, which could result in cause vast amounts of missing valid data (Chao and Zou, 2018). The POT model is based on the asymptotic convergence of the exceedances of a high threshold u to a generalized Pareto distribution (GPD) when u tends to infinity. It is generally considered the most useful for practical applications due to its more efficient use of data with extreme values. Hence, the GPD is chosen to analyze the marginal distributions of trigger indicators.
The GPD is a two-parameter distribution with distribution function:
Gξ,β(y)=1−(1+ξy/β)−1/ξ, | (6) |
where β>0, and the support is y≥0 when ξ≥0 and 0≤y≤−β/ξ when ξ<0.
This distribution is generalized in the sense that it subsumes certain other distributions under a common parametric form (McNeil and Frey, 2000). If ξ>0, then Gξ,β is a reparametrized version of the ordinary Pareto distribution; ξ<0 is known as a Pareto type Ⅱ distribution; ξ=0 corresponds to the exponential distribution, that is, G0,β=1−exp(−y/β).
Consider X1,X1,...,Xn to be a sequence of i.i.d. random variables with a cumulative continuous distribution function F(x). We choose a sufficiently large value threshold u, and the number of observations exceeding u are set as Nu. Yi=Xi−u denotes the excesses of Xi over the threshold u. Then, the distribution of Yi, known as the "conditional excess distribution function", is given by
Fu(y)=P(X−u≤y|X>u)=F(y+u)−F(u)1−F(u)=F(x)−F(u)1−F(u). | (7) |
When the threshold u is large enough, Fu(y) can be approximated by the GPD (see Balkema and Haan, 1974; Pickands, 1975). Substituting Equation (6) for Fu(y) into Equation (7), we obtain
F(x)=(1−F(u))(1−(1+ξ(x−u)/β)−1/ξ)+F(u),x>u, | (8) |
There is an empirical estimator of F(u), that is, (n−Nu)/n, which can be substituted into Equation (8) to obtain the following expression (see McNeil and Frey, 2000):
ˆF(x)=1−Nun(1+ξ(x−u)/β)−1/ξ,x>u. | (9) |
A copula is a function that factors a joint function into a marginal distribution and a dependence function C. The concept of copula goes back to Sklar (see Sklar, 1959). For two random variables X and Y, the joint distribution function is defined as
FXY(x,y)=P(X≤x,Y≤y). | (10) |
The corresponding survival function is defined as
¯FXY(x,y)=P(X>x,Y>y). | (11) |
The dependence function C is defined as
FXY(x,y)=C(w,v), | (12) |
where w=FX(x) and v=FY(y), and FX(x) and FY(y) are the marginal distribution functions of X and Y, respectively. If FX(x) and FY(y) are continuous, then the copula function C is unique. Conversely, when w, v are marginal distribution functions and C is a copula, the function F defined by Equation (12) is a multivariate distribution with marginal distributions FX(x) and FY(y) (Yao et al., 2017).
Another notation that is used in this paper is a survival copula of C, denoted by ¯C:
¯FXY(x,y)=¯C(¯FX(x),¯FY(y)). | (13) |
where ¯FX(x) and ¯FY(y) are the marginal survival functions of X and Y, respectively. The survival copula ¯C is related to the copula C by
¯C(w,v)=C(1−w,1−v)+w+v−1. | (14) |
Combining Equation (13) with Equation (14), we obtain
¯FXY(x,y)=C(FX(x),FY(y))−FX(x)−FY(y)+1. | (15) |
Our aim in this section is to prove the CAT bond pricing formula. We follow Merton (1976) to assume that the overall economy is only marginally influenced by localized catastrophic events and the catastrophe losses pertain to idiosyncratic shocks to the capital markets. That is, the catastrophic shocks will represent as "nonsystematic" and have a zero-risk premium. Under the risk-neutral probability measure Q, those events that depend only on financial variables are independent of those events that depend on catastrophic risk variables. As a result, the stochastic processes and distributions under Q retain the same characteristics as those under P (for a detailed discussion of this point, please refer to Cox and Pedersen, 2000; Lee and Yu, 2002; Braun, 2011; Ma and Ma, 2013). Assuming that the cash flows CF(t) on CAT bonds depend only on the catastrophic risk variables, then the price of CAT bonds can be calculated by the following formula:
P0(CF)=EQ[T∑t=1e−∫t0rsdsCF(t)]=p(0,t)EQ[T∑t=1CF(t)], | (16) |
where rs denotes the instantaneous interest rate and p(0,t)=EQ[e−∫t0rsds] stands for the price of a non-defaultable zero-coupon bond with a face value of 1 maturing at a time t.
In this paper, the instantaneous interest rate is set to be governed by a Cox-Ingersoll-Ross (CIR) model (see Cox et al., 1985). Then, the dynamics of the instantaneous interest rate under the risk-neutral probability measure Q can be described as
drt=κ(θ−rt)dt+ε√rtdWt. | (17) |
where κ>0 is the mean-reverting force measurement, ε>0 is the volatility parameter, θ>0 is the long-run mean of the interest rate, and Wt denotes a standard Wiener process under Q. The initial interest rate r0 is labeled as r; then, p(0,t) can be expressed as
p(0,t)=A(0,t)e−B(0,t)r, | (18) |
where
A(0,t)=[2ηe(κ+η)t/2(κ+η)(eηt−1)+2η]2κθ/ε2,B(0,t)=2(eηt−1)(κ+η)(eηt−1)+2η,η=√κ2+2ε2. |
Combining Equation (5) and Equation (16), the price of the bond at time t=0 is given by
P0(CF)=T∑t=1EQ[Ct]p(0,t)+EQ[FT]p(0,T)=T∑t=1FREQ[I{τC>t}]p(0,t)+FEQ[I{τF>T}]p(0,T)=F(RT∑t=1P(τC>t)p(0,t)+P(τF>T)p(0,T)). | (19) |
The probabilities that a single earthquake disaster satisfies the trigger conditions for a coupon and principal are marked as γC and γF, respectively. According to Equation (12) and Equation (15), we obtain
γC=1−C(FX(Trx),FY(Try)), | (20) |
and
γF=C(FX(Trx),FY(Try))−FX(Trx)−FY(Try)+1. | (21) |
The probability that the trigger condition for the coupon is satisfied in the interval [t−1,t] is marked as υC; then, it can be expressed as follows:
υC=1−EQ[Nt∏i=1(1−γC)]=1−∞∑n=0EQ[Nt∏i=1(1−γC)|Nt=n]P(Nt=n)=1−∞∑n=0EQ[n∏i=1(1−γC)]P(Nt=n)=1−∞∑n=0(1−γC)nλne−λn!=1−e−λ∞∑n=0[λ(1−γC)]nn!=1−e−λeλ(1−γC)=1−e−λγC. | (22) |
Similarly, the probability that the trigger condition for the principal is satisfied in the interval [0,T] is marked as υF. Since {Nt,t=1,2,...,T} is a sequence of i.i.d. random variables that obeys a Poisson distribution with an intensity λ, then ∑Tt=1Nt obeys a Poisson distribution with the intensity λT. Then, we can get
υF=1−EQ[T∑t=1Nt∏i=1(1−γF)]=1−∞∑n=0EQ[T∑t=1Nt∏i=1(1−γF)|T∑t=1Nt=n]P(T∑t=1Nt=n)=1−∞∑n=0EQ[n∏i=1(1−γF)]P(T∑t=1Nt=n)=1−∞∑n=0(1−γF)n(λT)ne−λTn!=1−e−λT∞∑n=0[λT(1−γF)]nn!=1−e−λTeλT(1−γF)=1−e−λTγF. | (23) |
Combining Equations (19)–(23), we can derive the CAT bond pricing formula as follows:
P0(CF)=F(RT∑t=1(1−υC)tp(0,t)+(1−υF)p(0,T))=F(RT∑t=1e−λtγCp(0,t)+e−λTγFp(0,T))=F(RT∑t=1e−λt(1−C(FX(Trx),FY(Try)))p(0,t)+e−λT(C(FX(Trx),FY(Try))−FX(Trx)−FY(Try)+1)p(0,T)). | (24) |
The data cover economy losses and magnitudes resulting from earthquake events on the Chinese mainland that occurred between 1990 and 2020. The data comes from the "Review of earthquake damage losses in the mainland of China (1990–2020)" compiled by the China Earthquake Networks Center; a total of 344 pairs of observations for economy losses and magnitudes were picked out. To eliminate the impact of time, the GDP for 2020 in China was chosen as the fixed base index to deal with the economy loss data.
In this section, we construct the marginal distributions of economy losses and magnitudes based on the POT model. In the field of EVT, threshold selection is critical, especially when the threshold-based methods are considered (Chukwudum et al., 2020). Hence, threshold selecting has to be done first before we can fit the GPD for the given data. We combine the mean excess function (MEF) plot and Hill plot to select the optimal threshold.
Figure 1 shows the MEF plot of economic losses, which is approximately linearly inclined upward, indicating that this group of data is a heavy tailed-distribution (ξ>0). Figure 2 is the Hill plot of economic losses, and the image tends to be stable roughly from the 50th point. To select the optimal threshold, the sequence statistics corresponding to each point are set as the threshold starting from the 45th point, and the maximum likelihood estimation is adopted for fitting. According to the results of Kolmogorov-Smirnov (K-S) testing, it can be known that the model fitting is optimal at the 48th point, that is, the threshold is uX=336975.39. Figure 3 shows the MEF plot of magnitudes, which is approximately linearly inclined downward, indicating that the data is light-tailed. Figure 4 is the Hill plot of magnitudes. Using a similar method to that used to determine the threshold of economic losses, the threshold of magnitudes can be determined as uY=6.6. Other parameter estimation results are shown in Table 1.
ξ | β | u | Nu | K-S statistics | p value | |
Economic losses | 1.1266 | 280471.11 | 336975.39 | 47 | 0.0659 | 0.9765 |
Magnitudes | −0.4789 | 0.8650 | 6.6 | 24 | 0.1926 | 0.2961 |
Substituting the estimated values of parameters into equation (9), the marginal distribution functions of economic losses and magnitudes can be obtained as follows:
FX(x)=1−47344(1+1.1266(x−336975.39)/280471.11)−1/1.1266, | (25) |
and
FY(y)=1−24344(1−0.4789(y−6.6)/0.8650)1/0.4789. | (26) |
Archimedean copula functions have the advantage of better describing the asymmetric correlation between variables and tail features. In this section, the maximum likelihood method is used to estimate the parameters of three commonly used binary Archimedean copula functions. The results are listed in Table 2. The empirical copula function is introduced to calculate the squared Euclidean distance between the estimation results for the three copula functions and the empirical copula function of the sample. The smaller the squared Euclidean distance, the better the fitness of the copula function (Shen et al., 2018; Gu et al., 2019; Cai et al., 2019). The empirical copula function is defined as
ˆC(w,v)=1nn∑i=1I{FX(xi)⩽w}I{FY(yi)⩽v}. | (27) |
Copula family | Function expression | Parameter α |
Frank copula | C(w,v)=−1αln{1+[exp(−αw)−1][exp(−αv)−1]exp(−α)−1}, α≠0 | 4.2634 |
Clayton copula | C(w,v)=(w−α+v−α−1)−1/α, α∈[−1,∞)\{0} | 1.0696 |
Gumbel copula | C(w,v)=exp{−[(−lnw)α+(−lnv)α]1/α}, α∈[1,∞) | 1.6176 |
Then, the squared Euclidean distance can be given by
d2=n∑i=1|ˆC(w,v)−C(w,v)|2. | (28) |
As shown in Table 3, the square Euclidean distance d2 of the Gumbel copula function is the smallest. Thus, the Gumbel copula function is the optimal joint distribution function; its function form is as follows:
C(w,v)=exp{−[(−lnw)1.6176+(−lnv)1.6176]1/1.6176}. | (29) |
Copula family | Squared Euclidean distance d2 |
Frank copula | 0.0773 |
Clayton copula | 0.2035 |
Gumbel copula | 0.0266 |
The correlation coefficient between economic losses and magnitudes can be measured by the Kendall rank correlation coefficient, which can be calculated as ρ=1−1/α=0.3818.
For a reference point of numerical analysis, we establish the set of parameters summarized in Table 4. The attachment points of the two trigger indicators are measured by the risk value VaRp, which refers to the maximum loss that a risk product may incur in a specified future period of time under the confidence level p. The risk value VaRP can be described as
VaRp(X)=F−1(p)=inf{x∈R:F(X)≥p},0<p<1. | (30) |
Parameter | Definitions | Base values |
T | Maturity (year) | 1–5 |
F | Principal (RMB) | 100 |
R | Coupon rate | 0.06 |
r | Initial interest rate value | 0.04 |
κ | Magnitude of mean-reverting force | 0.2 |
θ | Long-run mean of the interest rate | 0.05 |
ε | Volatility of the interest rate | 0.1 |
λ | Poisson intensity | 11 |
ρ | Kendall rank correlation coefficient | 0.3818 |
TrX | Attachment point of economic losses (million RMB) | 4824080.67 |
TrY | Attachment point of magnitudes | 7.7 |
Combining Equation (9) and Equation (30), the estimate of the risk value VaRp based on the GPD can be obtained as
VaRp(X)=u+βξ{[nNu(1−p)]−ξ−1}. | (31) |
In our case, the confidence level p is set equal to 99%. Then, the attachment points of the two trigger indicators can be calculated as shown in Table 4. The expectation and variance of the Poisson distribution are equal. Based on the data for the frequency of the earthquake disaster, we calculated the mean to be approximately equal to 11 for the estimation. Then, K-S testing was done; the result shows that the distribution of the frequency of earthquake disasters in China obeys the Poisson distribution. The interest rate parameters for the CIR model were set to r=0.04, κ=0.2, θ=0.05 and ε=0.1. These interest rate parameter values are all within the range typically used in previous studies (e.g., Lee and Yu, 2007; Lo, et al., 2013; Chao, 2021).
The prices of CAT bonds with maturities of 1–5 years are shown in Table 5. The CAT bond prices tend to decrease with the maturity period. The main reason is that, as the maturity period increases, there will be more catastrophe events that meet the trigger conditions, and the probability of losing coupons and the principal will be higher. The decline of the CAT bond price is slower than that of the principal; it is due to the coupon payment structure, that is, with the extension of the maturity period, the accumulated coupon amount also increases, making up for part of the principal loss.
Price | T =1 | T =2 | T =3 | T =4 | T =5 |
Cumulative discounted value of coupon | 4.87 | 8.06 | 10.00 | 11.03 | 11.41 |
Discounted value of principal | 91.17 | 83.00 | 75.50 | 68.63 | 62.37 |
CAT bond price | 96.04 | 91.06 | 85.50 | 79.66 | 73.78 |
The correlation between the two trigger indicators and the trigger levels are both important factors affecting the pricing of CAT bonds, which are shown in Figure 5. When the trigger levels of the two trigger indicators are increased simultaneously, that is, when the confidence level p of VaR is increased, the price is significantly increased. As can be seen in Figure 6, CAT bond prices are positively correlated with TrX and TrY. Raising the trigger levels means that the probability that the trigger conditions are met becomes smaller, and the possibility of the loss of coupons and principal is also reduced. Furthermore, we pay attention to the respective sensitivities of CAT bond prices to TrX and TrY.
Tables 6 and 7 show the changes in 1-year CAT bond prices after the trigger levels of the two trigger indicators were increased or decreased by 2.5% and 5%, respectively. It can be seen that, under the adjustment ratio, TrY has a greater impact on the CAT bond price.
Adjustment ratio for TrX | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+5% | +0.13% | +0.27% | +0.43% | +0.60% | +0.77% |
+2.5% | +0.06% | +0.14% | +0.22% | +0.30% | +0.39% |
-2.5% | -0.06% | -0.14% | -0.23% | -0.31% | -0.40% |
-5% | -0.13% | -0.29% | -0.46% | -0.64% | -0.82% |
Adjustment ratio for TrY | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+5% | +3.63% | +7.55% | +11.72% | +16.13% | +20.76% |
+2.5% | +1.77% | +3.73% | +5.84% | +8.06% | +10.35% |
-2.5% | -1.49% | -3.26% | -5.18% | -7.15% | -9.12% |
-5% | -2.72% | -6.06% | -9.61% | -13.15% | -16.53% |
It can also be seen in Figure 5 that the CAT bond price is negatively correlated with ρ. Table 8 shows the price of the 5-year CAT bonds at different values of the Kendall rank correlation coefficient ρ for the two trigger indicators. As ρ increases, the coupon discount changes less, while the principal discount decreases significantly. The stronger correlation between losses and magnitudes implies a greater probability of triggering two trigger indicators at the same time during a catastrophe event, and a greater possibility of principal damage.
Price | ρ=0.2818 | ρ=0.3818 | ρ=0.4818 | ρ=0.5818 | ρ=0.6818 |
Cumulative discounted value of coupon | 10.7437 | 11.4073 | 12.0641 | 12.7116 | 13.3477 |
Discounted value of principal | 66.2205 | 62.3683 | 58.9730 | 55.9689 | 53.3018 |
CAT bond price | 76.9643 | 73.7757 | 71.0370 | 68.6805 | 66.6495 |
Both the market interest rate and the coupon rate are important factors that directly affect the pricing of CAT bonds. According to the structure of the CIR model, the market interest rate is further divided into the initial interest rate and the long-run mean of the interest rate. As shown in Tables 9 and 10, CAT bond prices are negatively correlated with both r and θ, which is consistent with bond pricing theory. The adjustment of r has a greater impact on the pricing of CAT bonds, while the change of θ has little impact in the short term. It can be ascertained from Table 11 that there is a positive correlation between CAT bond price and R. Raising the coupon rate can increase the cash flow, and the CAT bond price will increase accordingly. Similarly, the adjustment of R has a greater impact on the price of multi-year CAT bonds, which is related to the pricing structure of CAT bonds. According to Table 5, the longer the maturity period, the larger the proportion of the cumulative discounted value of the coupon in the bond price, so the change of R will have a greater impact on the bond price.
Adjustment ratio for r | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+20% | -0.72% | -1.28% | -1.71% | -2.05% | -2.31% |
+10% | -0.36% | -0.64% | -0.86% | -1.03% | -1.16% |
-10% | +0.36% | +0.64% | +0.87% | +1.04% | +1.18% |
-20% | +0.73% | +1.29% | +1.74% | +2.09% | +2.37% |
Adjustment ratio for θ | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+20% | -0.09% | -0.34% | -0.70% | -1.14% | -1.65% |
+10% | -0.05% | -0.17% | -0.35% | -0.57% | -0.83% |
-10% | +0.06% | +0.17% | +0.35% | +0.57% | +0.83% |
-20% | +0.09% | +0.34% | +0.70% | +1.15% | +1.67% |
Adjustment ratio for R | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+20% | +1.01% | +1.77% | +2.34% | +2.77% | +3.09% |
+10% | +0.51% | +0.88% | +1.17% | +1.38% | +1.55% |
-10% | -0.51% | -0.88% | -1.17% | -1.38% | -1.55% |
-20% | -1.01% | -1.77% | -2.34% | -2.77% | -3.09% |
We developed a hybrid-triggered CAT bond pricing model and derived the pricing formula based on the copula-EVT model. We took earthquake CAT bonds as an example for model construction and numerical analysis. Taking economic loss and magnitude as dual trigger indicators, the advantages of the indemnity trigger and index trigger have been combined to increase the interest of investors with low-risk preference. The marginal distribution of the two trigger indicators were fitted to a POT model, and the joint distribution was established by using a copula function. Furthermore, we illustrated an application to the earthquake data for China from 1990 to 2020 and analyzed the sensitivity of CAT bond prices to model parameters.
The numerical experiments showed that the maturity, trigger level, dependence of trigger indicators, market interest rate and coupon rate have important implications for the CAT bond pricing model. First, the extension of the maturity period increases the probability of catastrophic events that meet the trigger conditions, and the CAT bond price decreases. In the development of multi-year CAT bond products, we can increase the yield of investors by raising the coupon rate or strengthening the trigger conditions so that CAT bonds will be more attractive in the investment market. Second, the stronger the correlation between the two trigger indicators, the more likely they will be triggered at the same time during a catastrophe event, which will weaken the advantages of the hybrid trigger mechanism. The Kendall rank correlation coefficient between economic loss and magnitude was 0.3818. It is appropriate to choose them as trigger indicators. Third, adjusting the trigger values of the two trigger indicators will change the price of CAT bonds in the opposite direction. Under the same adjustment ratio, a change in the magnitude trigger level has a greater impact on the price of CAT bonds. It is essential to select a reasonable attachment point of magnitude for the pricing of CAT bonds. Furthermore, the CAT bond price is positively correlated with the coupon rate and negatively correlated with the market interest rate. In the pricing of multi-year CAT bonds, the impact of a forward rate of interest on pricing should be fully considered.
There are several possible extensions and potential improvements in this direction. Firstly, the hybrid trigger indicator applied in our model is composed of an indemnity trigger indicator and a parametric trigger indicator. In the follow-up study, one can consider the application of other index triggers, such as the industry loss index and modeled loss index, in the construction of a hybrid-triggered CAT bond pricing model. Secondly, the periodicity and trend of the frequency of catastrophic disasters can be taken into consideration in our research. It is more realistic to apply the stochastic process with uncertain intensity to describe the frequency of catastrophic events. Thirdly, in addition to earthquakes, other disasters such as floods, typhoons and pandemics can also cause substantial economic losses. The characteristics of various disasters are different, and the CAT bonds with different catastrophe risks as the guarantee object have differences in the selection of trigger indicators and trigger thresholds, which will be the focus of future research.
This work was supported by the National Social Science Foundation of China (Grant No. 17BJY204), and the Ministry of Education of Humanities and Social Science Project (Grant No. 20YJA910001).
All authors declare no conflicts of interest in this paper.
[1] |
A. Adamu, A. Adam, Approximation of solutions of split equality fixed point problems with applications, Carpathian J. Math., 37 (2021), 381–392. https://doi.org/10.37193/CJM.2021.03.02 doi: 10.37193/CJM.2021.03.02
![]() |
[2] | R. Agarwal, D. O'Regan, D. Sahu, Iterative construction of fixed points of nearly asymptotically nonexpansive mappings, J. Nonlinear Convex Anal., 8 (2007), 61–79. |
[3] |
K. Aremu, L. Jolaoso, C. Izuchukwu, O. Mewomo, Approximation of common solution of finite family of monotone inclusion and fixed point problems for demicontractive multivalued mappings in CAT(0) spaces, Ricerche Mat., 69 (2020), 13–34. https://doi.org/10.1007/s11587-019-00446-y doi: 10.1007/s11587-019-00446-y
![]() |
[4] |
T. Bantaojai, C. Garodia, I. Uddin, N. Pakkaranang, P. Yimmuang, A novel iterative approach for solving common fixed point problems in geodesic spaces with convergence analysis, Carpathian J. Math., 37 (2021), 145–160. https://doi.org/10.37193/CJM.2021.02.01 doi: 10.37193/CJM.2021.02.01
![]() |
[5] |
I. Berg, I. Nikolaev, Quasilinearization and curvature of Aleksandrov spaces, Geom. Dedicata, 133 (2008), 195–218. https://doi.org/10.1007/s10711-008-9243-3 doi: 10.1007/s10711-008-9243-3
![]() |
[6] |
V. Berinde, Approximating fixed points of enriched nonexpansive mappings by Krasnoselskij iteration in Hilbert spaces, Carpathian J. Math., 35 (2019), 293–304. https://doi.org/10.37193/cjm.2019.03.04 doi: 10.37193/cjm.2019.03.04
![]() |
[7] |
V. Berinde, M. Păcurar, Approximating fixed points of enriched contractions in Banach spaces, J. Fixed Point Theory Appl., 22 (2020), 38. https://doi.org/10.1007/s11784-020-0769-9 doi: 10.1007/s11784-020-0769-9
![]() |
[8] |
V. Berinde, Approximating fixed points of demicontractive mappings via the quasi-nonexpansive case, Carpathian J. Math., 39 (2023), 73–84. https://doi.org/10.37193/cjm.2023.01.04 doi: 10.37193/cjm.2023.01.04
![]() |
[9] |
V. Berinde, M. Pǎcurar, Fixed points theorems for unsaturated and saturated classes of contractive mappings in Banach spaces, Symmetry, 13 (2021), 713. https://doi.org/10.3390/sym13040713 doi: 10.3390/sym13040713
![]() |
[10] | M. Bridson, A. Haefliger, Metric spaces of non-positive curvature, Berlin: Springer, 1999. https://doi.org/10.1007/978-3-662-12494-9 |
[11] |
F. Browder, Semicontractive and semiaccretive nonlinear mappings in Banach spaces, Bull. Amer. Math. Soc., 74 (1968), 660–665. https://doi.org/10.1090/s0002-9904-1968-11983-4 doi: 10.1090/s0002-9904-1968-11983-4
![]() |
[12] |
F. Browder, W. Petryshyn, The solution by iteration of nonlinear functional equations in Banach spaces, Bull. Amer. Math. Soc., 72 (1966), 571–575. https://doi.org/10.1090/s0002-9904-1966-11544-6 doi: 10.1090/s0002-9904-1966-11544-6
![]() |
[13] |
F. Browder, W. Petryshyn, Construction of fixed points of nonlinear mappings in Hilbert space, J. Math. Anal. Appl., 20 (1967), 197–228. https://doi.org/10.1016/0022-247x(67)90085-6 doi: 10.1016/0022-247x(67)90085-6
![]() |
[14] |
F. Bruhat, J. Tits, Groupes réductifs sur un corps local, Publications Mathématiques de l'Institut des Hautes Études Scientifiques, 41 (1972), 5–251. https://doi.org/10.1007/bf02715544 doi: 10.1007/bf02715544
![]() |
[15] |
K. Calderón, M. Khamsi, J. Martínez-Moreno, Perturbed approximations of fixed points of nonexpansive mappings in CATp(0) spaces, Carpathian J. Math., 37 (2021), 65–79. https://doi.org/10.37193/CJM.2021.01.07 doi: 10.37193/CJM.2021.01.07
![]() |
[16] |
M. Choban, About convex structures on metric spaces, Carpathian J. Math., 38 (2022), 391–404. https://doi.org/10.37193/cjm.2022.02.10 doi: 10.37193/cjm.2022.02.10
![]() |
[17] |
S. Dhompongsa, W. Kirk, B. Sims, Fixed points of uniformly Lipschitzian mappings, Nonlinear Anal., 65 (2006), 762–772. https://doi.org/10.1016/j.na.2005.09.044 doi: 10.1016/j.na.2005.09.044
![]() |
[18] |
S. Dhompongsa, B. Panyanak, On Δ-convergence theorems in CAT (0) spaces, Comput. Math. Appl., 56 (2008), 2572–2579. https://doi.org/10.1016/j.camwa.2008.05.036 doi: 10.1016/j.camwa.2008.05.036
![]() |
[19] |
W. Dotson, Fixed points of quasi-nonexpansive mappings, J. Aust. Math. Soc., 13 (1972), 167–170. https://doi.org/10.1017/s144678870001123x doi: 10.1017/s144678870001123x
![]() |
[20] |
M. Eslamian, General algorithms for split common fixed point problem of demicontractive mappings, Optimization, 65 (2016), 443–465. https://doi.org/10.1080/02331934.2015.1053883 doi: 10.1080/02331934.2015.1053883
![]() |
[21] |
Q. Fan, J. Peng, H. He, Weak and strong convergence theorems for the split common fixed point problem with demicontractive operators, Optimization, 70 (2021), 1409–1423. https://doi.org/10.1080/02331934.2021.1890074 doi: 10.1080/02331934.2021.1890074
![]() |
[22] |
A. Hanjing, S. Suantai, The split fixed point problem for demicontractive mappings and applications, Fixed Point Theory, 21 (2020), 507–524. https://doi.org/10.24193/fpt-ro.2020.2.37 doi: 10.24193/fpt-ro.2020.2.37
![]() |
[23] |
A. Hanjing, S. Suantai, Y. Cho, Hybrid inertial accelerated extragradient algorithms for split pseudomonotone equilibrium problems and fixed point problems of demicontractive mappings, Filomat, 37 (2023), 1607–1623. https://doi.org/10.2298/FIL2305607H doi: 10.2298/FIL2305607H
![]() |
[24] |
T. Hicks, J. Kubicek, On the Mann iteration process in a Hilbert space, J. Math. Anal. Appl., 59 (1977), 498–504. https://doi.org/10.1016/0022-247x(77)90076-2 doi: 10.1016/0022-247x(77)90076-2
![]() |
[25] | A. Inuwa, P. Chaipunya, P. Kumam, S. Salisu, Equilibrium problems and proximal algorithm using tangent space products, Carpathian J. Math., in press. |
[26] |
S. Ishikawa, Fixed points by a new iteration method, Proc. Amer. Math. Soc., 44 (1974), 147–150. https://doi.org/10.1090/s0002-9939-1974-0336469-5 doi: 10.1090/s0002-9939-1974-0336469-5
![]() |
[27] |
S. Ishikawa, Fixed points and iteration of a nonexpansive mapping in a Banach space, Proc. Amer. Math. Soc., 59 (1976), 65–71. https://doi.org/10.1090/s0002-9939-1976-0412909-x doi: 10.1090/s0002-9939-1976-0412909-x
![]() |
[28] |
P. Jailoka, S. Suantai, Split null point problems and fixed point problems for demicontractive multivalued mappings, Mediterr. J. Math., 15 (2018), 204. https://doi.org/10.1007/s00009-018-1251-4 doi: 10.1007/s00009-018-1251-4
![]() |
[29] |
P. Jailoka, S. Suantai, Split common fixed point and null point problems for demicontractive operators in Hilbert spaces, Optim. Method. Softw., 34 (2019), 248–263. https://doi.org/10.1080/10556788.2017.1359265 doi: 10.1080/10556788.2017.1359265
![]() |
[30] |
P. Jailoka, S. Suantai, The split common fixed point problem for multivalued demicontractive mappings and its applications, RACSAM, 113 (2019), 689–706. https://doi.org/10.1007/s13398-018-0496-x doi: 10.1007/s13398-018-0496-x
![]() |
[31] |
K. Juagwon, W. Phuengrattana, Iterative approaches for solving equilibrium problems, zero point problems and fixed point problems in Hadamard spaces, Comp. Appl. Math., 42 (2023), 75. https://doi.org/10.1007/s40314-023-02209-w doi: 10.1007/s40314-023-02209-w
![]() |
[32] |
Y. Kimura, Resolvents of equilibrium problems on a complete geodesic space with curvature bounded above, Carpathian J. Math., 37 (2021), 463–476. https://doi.org/10.37193/CJM.2021.03.09 doi: 10.37193/CJM.2021.03.09
![]() |
[33] | W. Kirk, N. Shahzad, Fixed point theory in distance spaces, Cham: Springer, 2014. https://doi.org/10.1007/978-3-319-10927-5 |
[34] | M. Krasnosel'skii, Two remarks on the method of successive approximations, Uspekhi Matematicheskikh Nauk, 10 (1955), 123–127. |
[35] |
W. Mann, Mean value methods in iteration, Proc. Amer. Math. Soc., 4 (1953), 506–510. https://doi.org/10.1090/s0002-9939-1953-0054846-3 doi: 10.1090/s0002-9939-1953-0054846-3
![]() |
[36] |
M. Minjibir, S. Salisu, Strong and Δ-convergence theorems for a countable family of multivalued demicontractive maps in Hadamard spaces, Nonlinear Functional Analysis and Applications, 27 (2022), 45–58. https://doi.org/10.22771/nfaa.2022.27.01.03 doi: 10.22771/nfaa.2022.27.01.03
![]() |
[37] | C. Moore, Iterative approximation fixed points of demicontractive maps, Proceedings of The Abdus Salam International Centre for Theoretical Physics, 1998, 1–12. |
[38] |
Ş. Mǎruşter, Sur le calcul des zéros d'un opérateur discontinu par itération, Can. Math. Bull., 16 (1973), 541–544. https://doi.org/10.4153/cmb-1973-088-7 doi: 10.4153/cmb-1973-088-7
![]() |
[39] |
Ş. Mǎruşter, The solution by iteration of nonlinear equations in Hilbert spaces, Proc. Amer. Math. Soc., 63 (1977), 69–73. https://doi.org/10.1090/s0002-9939-1977-0636944-2 doi: 10.1090/s0002-9939-1977-0636944-2
![]() |
[40] |
M. Nnakwe, J. Ezeora, Strong convergence theorems for variational inequalities and fixed point problems in Banach spaces, Carpathian J. Math., 37 (2021), 477–487. https://doi.org/10.37193/cjm.2021.03.10 doi: 10.37193/cjm.2021.03.10
![]() |
[41] |
M. Noor, New approximation schemes for general variational inequalities, J. Math. Anal. Appl., 251 (2000), 217–229. https://doi.org/10.1006/jmaa.2000.7042 doi: 10.1006/jmaa.2000.7042
![]() |
[42] |
F. Ogbuisi, F. Isiogugu, A new iterative algorithm for pseudomonotone equilibrium problem and a finite family of demicontractive mappings, Abstr. Appl. Anal., 2020 (2020), 3183529. https://doi.org/10.1155/2020/3183529 doi: 10.1155/2020/3183529
![]() |
[43] |
A. Padcharoen, P. Kumam, Y. Cho, Split common fixed point problems for demicontractive operators, Numer. Algor., 82 (2019), 297–320. https://doi.org/10.1007/s11075-018-0605-0 doi: 10.1007/s11075-018-0605-0
![]() |
[44] |
S. Panja, M. Saha, R. Bisht, Existence of common fixed points of non-linear semigroups of enriched Kannan contractive mappings, Carpathian J. Math., 38 (2022), 169–178. https://doi.org/10.37193/cjm.2022.01.14 doi: 10.37193/cjm.2022.01.14
![]() |
[45] |
S. Salisu, M. Minjibir, P. Kumam, S. Sriwongsa, Convergence theorems for fixed points in CATp(0) spaces, J. Appl. Math. Comput., 69 (2023), 631–650. https://doi.org/10.1007/s12190-022-01763-6 doi: 10.1007/s12190-022-01763-6
![]() |
[46] |
S. Salisu, P. Kumam, S. Sriwongsa, On fixed points of enriched contractions and enriched nonexpansive mappings, Carpathian J. Math., 39 (2022), 237–254. https://doi.org/10.37193/cjm.2023.01.16 doi: 10.37193/cjm.2023.01.16
![]() |
[47] |
S. Salisu, P. Kumam, S. Sriwongsa, Strong convergence theorems for fixed point of multi-valued mappings in Hadamard spaces, J. Inequal. Appl., 2022 (2022), 143. https://doi.org/10.1186/s13660-022-02870-5 doi: 10.1186/s13660-022-02870-5
![]() |
[48] |
S. Salisu, P. Kumam, S. Sriwongsa, Convergence theorems for monotone vector field inclusions and minimization problems in Hadamard spaces, Anal. Geom. Metr. Space., 11 (2023), 20220150. https://doi.org/10.1515/agms-2022-0150 doi: 10.1515/agms-2022-0150
![]() |
[49] |
S. Suantai, P. Jailoka, A self-adaptive algorithm for split null point problems and fixed point problems for demicontractive multivalued mappings, Acta Appl. Math., 170 (2020), 883–901. https://doi.org/10.1007/s10440-020-00362-6 doi: 10.1007/s10440-020-00362-6
![]() |
[50] |
R. Suparatulatorn, P. Cholamjiak, S. Suantai, Self-adaptive algorithms with inertial effects for solving the split problem of the demicontractive operators, RACSAM, 114 (2020), 40. https://doi.org/10.1007/s13398-019-00737-x doi: 10.1007/s13398-019-00737-x
![]() |
[51] | A. Wang, J. Zhao, Self-adaptive iterative algorithms for the split common fixed point problem with demicontractive operators, J. Nonlinear Var. Anal., 5 (2021), 573–587. |
[52] |
Y. Wang, X. Fang, J. Guan, T. Kim, On split null point and common fixed point problems for multivalued demicontractive mappings, Optimization, 70 (2021), 1121–1140. https://doi.org/10.1080/02331934.2020.1764952 doi: 10.1080/02331934.2020.1764952
![]() |
[53] |
F. Wang, The split feasibility problem with multiple output sets for demicontractive mappings, J. Optim. Theory Appl., 195 (2022), 837–853. https://doi.org/10.1007/s10957-022-02096-x doi: 10.1007/s10957-022-02096-x
![]() |
[54] |
J. Xiao, Y. Wang, A viscosity method with inertial effects for split common fixed point problems of demicontractive mappings, J. Nonlinear Funct. Anal., 2022 (2022), 17. https://doi.org/10.23952/jnfa.2022.17 doi: 10.23952/jnfa.2022.17
![]() |
[55] |
Y. Yao, Y. Liou, M. Postolache, Self-adaptive algorithms for the split problem of the demicontractive operators, Optimization, 67 (2018), 1309–1319. https://doi.org/10.1080/02331934.2017.1390747 doi: 10.1080/02331934.2017.1390747
![]() |
[56] |
Y. Yao, J. Yao, Y. Liou, M. Postolache, Iterative algorithms for split common fixed points of demicontractive operators without priori knowledge of operator norms, Carpathian J. Math., 34 (2018), 459–466. https://doi.org/10.37193/cjm.2018.03.23 doi: 10.37193/cjm.2018.03.23
![]() |
[57] | C. Zhang, Y. Li, Y. Wang, On solving split generalized equilibrium problems with trifunctions and fixed point problems of demicontractive multi-valued mappings, J. Nonlinear Convex Anal., 21 (2020), 2027–2042. |
1. | Wulan Anggraeni, Sudradjat Supian, Nurfadhlina Binti Abdul Halim, Earthquake Catastrophe Bond Pricing Using Extreme Value Theory: A Mini-Review Approach, 2022, 10, 2227-7390, 4196, 10.3390/math10224196 | |
2. | Riza Andrian Ibrahim, Moch Panji Agung Saputra, Yuyun Hidayat, Hafizan Juahir, Igif Gimin Prihanto, Nurfadhlina Binti Abdul Halim, Modeling Multiple-Event Catastrophe Bond Prices Involving the Trigger Event Correlation, Interest, and Inflation Rates, 2022, 10, 2227-7390, 4685, 10.3390/math10244685 | |
3. | Jiayi Li, Zhiyan Cai, Yixuan Liu, Chengxiu Ling, Extremal Analysis of Flooding Risk and Its Catastrophe Bond Pricing, 2022, 11, 2227-7390, 114, 10.3390/math11010114 | |
4. | Wulan Anggraeni, Sudradjat Supian, Nurfadhlina Abdul Halim, Single Earthquake Bond Pricing Framework with Double Trigger Parameters Based on Multi Regional Seismic Information, 2023, 11, 2227-7390, 689, 10.3390/math11030689 | |
5. | Yifan Tang, Conghua Wen, Chengxiu Ling, Yuqing Zhang, Pricing Multi-Event-Triggered Catastrophe Bonds Based on a Copula–POT Model, 2023, 11, 2227-9091, 151, 10.3390/risks11080151 | |
6. | Wulan Anggraeni, Sudradjat Supian, Nurfadhlina Abdul Halim, Catastrophe Bond Diversification Strategy Using Probabilistic–Possibilistic Bijective Transformation and Credibility Measures in Fuzzy Environment, 2023, 11, 2227-7390, 3513, 10.3390/math11163513 | |
7. | Riza Andrian Ibrahim, Herlina Napitupulu, Rose Irnawaty Ibrahim, How to Price Catastrophe Bonds for Sustainable Earthquake Funding? A Systematic Review of the Pricing Framework, 2023, 15, 2071-1050, 7705, 10.3390/su15097705 | |
8. | Riza Andrian Ibrahim, Herlina Napitupulu, Rose Irnawaty Ibrahim, Earthquake Bond Pricing Model Involving the Inconstant Event Intensity and Maximum Strength, 2024, 12, 2227-7390, 786, 10.3390/math12060786 | |
9. | Riza Andrian Ibrahim, Herlina Napitupulu, Rose Irnawaty Ibrahim, Modeling earthquake bond prices with correlated dual trigger indices and the approximate solution using the Monte Carlo algorithm, 2025, 10, 2473-6988, 2223, 10.3934/math.2025103 |
ξ | β | u | Nu | K-S statistics | p value | |
Economic losses | 1.1266 | 280471.11 | 336975.39 | 47 | 0.0659 | 0.9765 |
Magnitudes | −0.4789 | 0.8650 | 6.6 | 24 | 0.1926 | 0.2961 |
Copula family | Function expression | Parameter α |
Frank copula | C(w,v)=−1αln{1+[exp(−αw)−1][exp(−αv)−1]exp(−α)−1}, α≠0 | 4.2634 |
Clayton copula | C(w,v)=(w−α+v−α−1)−1/α, α∈[−1,∞)\{0} | 1.0696 |
Gumbel copula | C(w,v)=exp{−[(−lnw)α+(−lnv)α]1/α}, α∈[1,∞) | 1.6176 |
Copula family | Squared Euclidean distance d2 |
Frank copula | 0.0773 |
Clayton copula | 0.2035 |
Gumbel copula | 0.0266 |
Parameter | Definitions | Base values |
T | Maturity (year) | 1–5 |
F | Principal (RMB) | 100 |
R | Coupon rate | 0.06 |
r | Initial interest rate value | 0.04 |
κ | Magnitude of mean-reverting force | 0.2 |
θ | Long-run mean of the interest rate | 0.05 |
ε | Volatility of the interest rate | 0.1 |
λ | Poisson intensity | 11 |
ρ | Kendall rank correlation coefficient | 0.3818 |
TrX | Attachment point of economic losses (million RMB) | 4824080.67 |
TrY | Attachment point of magnitudes | 7.7 |
Price | T =1 | T =2 | T =3 | T =4 | T =5 |
Cumulative discounted value of coupon | 4.87 | 8.06 | 10.00 | 11.03 | 11.41 |
Discounted value of principal | 91.17 | 83.00 | 75.50 | 68.63 | 62.37 |
CAT bond price | 96.04 | 91.06 | 85.50 | 79.66 | 73.78 |
Adjustment ratio for TrX | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+5% | +0.13% | +0.27% | +0.43% | +0.60% | +0.77% |
+2.5% | +0.06% | +0.14% | +0.22% | +0.30% | +0.39% |
-2.5% | -0.06% | -0.14% | -0.23% | -0.31% | -0.40% |
-5% | -0.13% | -0.29% | -0.46% | -0.64% | -0.82% |
Adjustment ratio for TrY | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+5% | +3.63% | +7.55% | +11.72% | +16.13% | +20.76% |
+2.5% | +1.77% | +3.73% | +5.84% | +8.06% | +10.35% |
-2.5% | -1.49% | -3.26% | -5.18% | -7.15% | -9.12% |
-5% | -2.72% | -6.06% | -9.61% | -13.15% | -16.53% |
Price | ρ=0.2818 | ρ=0.3818 | ρ=0.4818 | ρ=0.5818 | ρ=0.6818 |
Cumulative discounted value of coupon | 10.7437 | 11.4073 | 12.0641 | 12.7116 | 13.3477 |
Discounted value of principal | 66.2205 | 62.3683 | 58.9730 | 55.9689 | 53.3018 |
CAT bond price | 76.9643 | 73.7757 | 71.0370 | 68.6805 | 66.6495 |
Adjustment ratio for r | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+20% | -0.72% | -1.28% | -1.71% | -2.05% | -2.31% |
+10% | -0.36% | -0.64% | -0.86% | -1.03% | -1.16% |
-10% | +0.36% | +0.64% | +0.87% | +1.04% | +1.18% |
-20% | +0.73% | +1.29% | +1.74% | +2.09% | +2.37% |
Adjustment ratio for θ | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+20% | -0.09% | -0.34% | -0.70% | -1.14% | -1.65% |
+10% | -0.05% | -0.17% | -0.35% | -0.57% | -0.83% |
-10% | +0.06% | +0.17% | +0.35% | +0.57% | +0.83% |
-20% | +0.09% | +0.34% | +0.70% | +1.15% | +1.67% |
Adjustment ratio for R | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+20% | +1.01% | +1.77% | +2.34% | +2.77% | +3.09% |
+10% | +0.51% | +0.88% | +1.17% | +1.38% | +1.55% |
-10% | -0.51% | -0.88% | -1.17% | -1.38% | -1.55% |
-20% | -1.01% | -1.77% | -2.34% | -2.77% | -3.09% |
ξ | β | u | Nu | K-S statistics | p value | |
Economic losses | 1.1266 | 280471.11 | 336975.39 | 47 | 0.0659 | 0.9765 |
Magnitudes | −0.4789 | 0.8650 | 6.6 | 24 | 0.1926 | 0.2961 |
Copula family | Function expression | Parameter α |
Frank copula | C(w,v)=−1αln{1+[exp(−αw)−1][exp(−αv)−1]exp(−α)−1}, α≠0 | 4.2634 |
Clayton copula | C(w,v)=(w−α+v−α−1)−1/α, α∈[−1,∞)\{0} | 1.0696 |
Gumbel copula | C(w,v)=exp{−[(−lnw)α+(−lnv)α]1/α}, α∈[1,∞) | 1.6176 |
Copula family | Squared Euclidean distance d2 |
Frank copula | 0.0773 |
Clayton copula | 0.2035 |
Gumbel copula | 0.0266 |
Parameter | Definitions | Base values |
T | Maturity (year) | 1–5 |
F | Principal (RMB) | 100 |
R | Coupon rate | 0.06 |
r | Initial interest rate value | 0.04 |
κ | Magnitude of mean-reverting force | 0.2 |
θ | Long-run mean of the interest rate | 0.05 |
ε | Volatility of the interest rate | 0.1 |
λ | Poisson intensity | 11 |
ρ | Kendall rank correlation coefficient | 0.3818 |
TrX | Attachment point of economic losses (million RMB) | 4824080.67 |
TrY | Attachment point of magnitudes | 7.7 |
Price | T =1 | T =2 | T =3 | T =4 | T =5 |
Cumulative discounted value of coupon | 4.87 | 8.06 | 10.00 | 11.03 | 11.41 |
Discounted value of principal | 91.17 | 83.00 | 75.50 | 68.63 | 62.37 |
CAT bond price | 96.04 | 91.06 | 85.50 | 79.66 | 73.78 |
Adjustment ratio for TrX | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+5% | +0.13% | +0.27% | +0.43% | +0.60% | +0.77% |
+2.5% | +0.06% | +0.14% | +0.22% | +0.30% | +0.39% |
-2.5% | -0.06% | -0.14% | -0.23% | -0.31% | -0.40% |
-5% | -0.13% | -0.29% | -0.46% | -0.64% | -0.82% |
Adjustment ratio for TrY | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+5% | +3.63% | +7.55% | +11.72% | +16.13% | +20.76% |
+2.5% | +1.77% | +3.73% | +5.84% | +8.06% | +10.35% |
-2.5% | -1.49% | -3.26% | -5.18% | -7.15% | -9.12% |
-5% | -2.72% | -6.06% | -9.61% | -13.15% | -16.53% |
Price | ρ=0.2818 | ρ=0.3818 | ρ=0.4818 | ρ=0.5818 | ρ=0.6818 |
Cumulative discounted value of coupon | 10.7437 | 11.4073 | 12.0641 | 12.7116 | 13.3477 |
Discounted value of principal | 66.2205 | 62.3683 | 58.9730 | 55.9689 | 53.3018 |
CAT bond price | 76.9643 | 73.7757 | 71.0370 | 68.6805 | 66.6495 |
Adjustment ratio for r | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+20% | -0.72% | -1.28% | -1.71% | -2.05% | -2.31% |
+10% | -0.36% | -0.64% | -0.86% | -1.03% | -1.16% |
-10% | +0.36% | +0.64% | +0.87% | +1.04% | +1.18% |
-20% | +0.73% | +1.29% | +1.74% | +2.09% | +2.37% |
Adjustment ratio for θ | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+20% | -0.09% | -0.34% | -0.70% | -1.14% | -1.65% |
+10% | -0.05% | -0.17% | -0.35% | -0.57% | -0.83% |
-10% | +0.06% | +0.17% | +0.35% | +0.57% | +0.83% |
-20% | +0.09% | +0.34% | +0.70% | +1.15% | +1.67% |
Adjustment ratio for R | T = 1 | T = 2 | T = 3 | T = 4 | T = 5 |
+20% | +1.01% | +1.77% | +2.34% | +2.77% | +3.09% |
+10% | +0.51% | +0.88% | +1.17% | +1.38% | +1.55% |
-10% | -0.51% | -0.88% | -1.17% | -1.38% | -1.55% |
-20% | -1.01% | -1.77% | -2.34% | -2.77% | -3.09% |