Next Article in Journal
Constraint Qualifications for Vector Optimization Problems in Real Topological Spaces
Next Article in Special Issue
On the Study of Pseudo 𝒮-Asymptotically Periodic Mild Solutions for a Class of Neutral Fractional Delayed Evolution Equations
Previous Article in Journal
Moderate Deviation Principle for Linear Processes Generated by Dependent Sequences under Sub-Linear Expectation
Previous Article in Special Issue
Determinantal Expressions, Identities, Concavity, Maclaurin Power Series Expansions for van der Pol Numbers, Bernoulli Numbers, and Cotangent
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pricing of Credit Risk Derivatives with Stochastic Interest Rate

College of Science, Donghua University, Shanghai 201620, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(8), 782; https://doi.org/10.3390/axioms12080782
Submission received: 15 June 2023 / Revised: 27 July 2023 / Accepted: 9 August 2023 / Published: 12 August 2023

Abstract

:
This paper deals with a credit derivative pricing problem using the martingale approach. We generalize the conventional reduced-form credit risk model for a credit default swap market, assuming that the firms’ default intensities depend on the default states of counterparty firms and that the stochastic interest rate follows a jump-diffusion Cox–Ingersoll–Ross process. First, we derive the joint Laplace transform of the distribution of the vector process ( r t , R t ) by applying piecewise deterministic Markov process theory and martingale theory. Then, using the joint Laplace transform, we obtain the explicit pricing of defaultable bonds and a credit default swap. Lastly, numerical examples are presented to illustrate the dynamic relationships between defaultable securities (defaultable bonds, credit default swap) and the maturity date.

1. Introduction

Efficient tools for the management of credit risk in portfolios are always in high demand in financial institutions. This leads to a specialized market for credit derivatives. Due to the initiation of the credit derivatives market in the early 1990s, credit derivatives have maintained robust growth in the past decades. Being an enormous market that already surpasses the debt one, it has been attracting increasing research interest. For the modeling of credit risk in credit-risky securities markets, depending on the formulation, there are two well-established fundamental models available at present, i.e. the so-called reduced-form model [1,2,3] and the structural model [4,5]. In the reduced-form model, the influences on the default are assumed to be dominated by exogenous factors, for example, policies from the government; meanwhile, structural models mainly focus on endogenous factors—more specifically, the asset allocation of firms. Compared with the structural one, the reduced-form model is believed to be easier for constructing a tractable formula explaining the price of credit-risky securities in terms of economic covariates, which facilitates the estimation. In this work, we mainly focus on the reduced-form model.
There is a long history of insights into credit risk in the credit derivatives market with the reduced-form model, going back to the 1990s. When describing dependency default risks, bottom-up models as well as top-down models are commonly used. Bottom-up models focus on modeling default intensities of individual reference entities to characterize a portfolio default intensity. For more developments of bottom-up models, references [6,7,8] are highly recommended. The top-down models mainly focus on modeling default at the portfolio level; in other words, the default intensity of the entire investment portfolio is modeled without reference to the constituent names. Some works on top-down models are included in [9,10,11]. In our paper, we focus on bottom-up models.
In bottom-up models, the default intensity is generally assumed to be a stochastic process, while the relevant randomness is built up from state variables. Moreover, a further consideration of credit contagion is then pioneered by [12] (the so-called DL model) to incorporate the concentration risk in credit portfolios. In the meantime, the emergence of certain special phenomena—for instance, the banking crisis—reveals one of the defects of the conventional reduced-form models, which did not possess a full consideration of credit risk sources. Accordingly, a model with the concept of counterparty risk was brought forward by [7]. In recent years many scholars have considered this factor, for example in [13,14,15,16,17,18]. In conjunction with the credit contagion, this model emphasizes the so-called primary–secondary framework, which simplifies the payoff structure. The default intensity in such a framework is assumed to be mainly decided by the default state of the counterparty and the economy-wide state variables that follow the Vasicek process (see e.g., [19]). In addition, [13,15,16,20,21] also involved the effects of various other factors (such as the settlement risk, the replacement cost, credit rating, Markov Regime-Switching) in the pricing of CDSs. However, further effects are still desired for improving these models to better characterize the credit derivatives market.
One of the essential aspects of a fine modeling of credit risk is the proper addressing of the dynamic changes of economy-wide state variables. Among all economy-wide state variables, we only consider the risk-free interest rate in this work. As the mean reversion may come from the underlying macro-economic currents, or the corrections on account of bond market overreaction, the Cox–Ingersoll–Ross (CIR) process (for more details, see [22]) is applied to characterize the interest rate fluctuation [23,24,25,26,27]. This process covers the defect of the Vasicek model that makes it possible to get negative interest rates. In the CIR process, the volatility term θ r ( t ) will approach zero when the interest rate gets close to zero, which cancels the effects of the randomness. Thus, the interest rate in such a case remains positive. Meanwhile, the volatility takes high value when the interest rate is high, which is a desired property for modeling. Moreover, in practice, when a primary event takes place, e.g., a change in monetary policies, there will be a positive jump in the interest rate process. Such a jump will be suppressed by the government over time. Then, the process continues until it is taken over by another jump caused by a new event. Therefore, to depict such a phenomenon, we consider a jump-diffusion Cox–Ingersoll–Ross (JCIR) process.
In this paper, we propose a revised reduced-form model [28]. It incorporates the assumption that the firms’ default intensities are related to the default states of counterparty firms and the stochastic interest rate driven by the jump-diffusion Cox–Ingersoll–Ross (JCIR) process. This model, to some extent, is more in accordance with the actual CDS market than that of classical work [7]. The primary intention of this article is to derive an explicit expression for the joint Laplace transform of the distribution of the risk-free interest rate-related vector process ( r t , R t ) by applying the piecewise deterministic Markov process theory. The piecewise deterministic Markov process theory was developed by [29] and has been proven to be a very powerful mathematical tool in the field of non-diffusion models. Using the joint Laplace transform, we obtain the explicit pricing of defaultable bonds and default swap rates. Finally, according to this explicit expression, we present the dynamic relations between defaultable securities and the maturity date.
The rest of this paper is organized as follows. In Section 2, we present the setup of the counterparty structure. In Section 3, we derive the joint Laplace transform of the distribution of the vector process ( r t , R t ) , and the pricing formulae for defaultable bonds are derived. Meanwhile, with several examples, we also explore the impact of counterparty risk and spot interest rate on bond pricing. In Section 4, we extend our model to the pricing of credit default swaps and the importance of default correlation and spot interest rate models in pricing a default swap. In Section 5, conclusive remarks and further academic challenges are presented.

2. Reduced-Form Model

In this section, to construct the model, a collection of doubly stochastic Poisson processes are assumed to represent the default. Then, the primary–secondary framework is applied to further simplify the subsequent exposition.

2.1. Model Setup

To set up the model, the filtered probability space ( Ω , F , { F } t = 0 T * , P ) is used for the uncertainty in the economy, while X t , an R d -valued process, and point processes N i ( i = 1 , , I ) are assumed to indicate the d economy-wide state variables and the default processes of firm i, respectively. N i is 0 initially and changes to 1 once the default takes place. Denoting τ i as the relevant default time, the default process of firm i could be presented as N t i = 1 { τ i t } . Then, the filtration could be generated according to X t and N t i as
F t = F t X F t 1 F t I ,
with
F t X = σ ( X s , 0 s t ) and F t i = σ ( N s i , 0 s t )
being the filtrations obtained from X t and N t i , respectively.
In addition, we indicate the filtration generated by the default process of all firms but the ith one as
F t i = F t 1 F t i 1 F t i + 1 F t I .
Thus, F T * X F T * i embraces the full information up to time T * on the state variables and the default processes of all firms other than that of the ith. We could also choose an intensity process λ t i , for t [ 0 , T * ] , being non-negative and F T * X F T * i -measurable, with
0 t λ s i d s < , P a . s . ,
so that a Poisson process N i can be defined.

2.2. Primary–Secondary Framework

As a typical case, in this work, we discuss the default probability of a commercial bank. We assume that bank A owns a large amount of debt from firm B. However, it is normally not vice versa, and thus not possible to affect B’s default probability. Thus, we implement the constraints as follows. All I firms considered could be separated into two incompatible types, i.e., S 1 primary firms and the rest I S 1 secondary firms. The default processes of the primary firms only depend on state variables, while those of secondary firms depend on both the state variables and the default status of the primary firms.
To start with, under probability measure P , for the S 1 primary firms, there is a sequence of independent unit exponential random variables, { e i , 1 i S 1 } , being also independent of X t . Thus, the relevant default times can be defined as
τ i = inf { t : 0 t λ s i d s e i } , 1 i S 1 ,
with λ t i adapted to F t X , while the conditional and unconditional survival probability distributions of primary firm i are given by
P ( τ i > t | F T * X ) = exp ( 0 t λ s i d s ) , t [ 0 , T * ] ,
P ( τ i > t ) = E [ exp ( 0 t λ s i d s ) ] , t [ 0 , T * ] .
Since the primary firms’ default processes only depend on X t , we denote their default intensities by
λ t i = a 0 , t i , 1 i S 1 .
Then, we proceed to the secondary default processes. Similarly, for the I S 1 secondary firms, there is also a sequence of relevant independent variables, { e j , S 1 + 1 j I } , being independent of X t and { τ i , 1 i S 1 } as well. Hence, the default times can be defined as
τ j = inf { t : 0 t λ s j d s e j } , S 1 + 1 j I ,
with λ t j adapted to F t X F t 1 F t S 1 , while the conditional and unconditional survival probability distributions of secondary firm j are expressed as
P ( τ j > t | F T * X F T * 1 F T * S 1 ) = exp ( 0 t λ s j d s ) , t [ 0 , T * ] ,
P ( τ j > t ) = E [ exp ( 0 t λ s j d s ) ] , t [ 0 , T * ] .
Since the default processes of secondary firms depend on X t and the default processes of primary firms, we denote the default intensities by
λ t j = a 0 , t j + k = 1 S 1 a k , t j 1 { τ k t } , S 1 + 1 j I ,
with a k , t j adapted to F t X . The default of firm k, according to the sign of a k , t j , can either increase ( a k , t j > 0 ) or decrease ( a k , t j < 0 ) the default probability of firm j. In addition, a more general default intensity can be achieved by involving the interactions of the indicator functions.
According to the assumptions above, we now derive the explicit pricing formulae of defaultable bonds and default swap rates in Section 3 and Section 4, respectively.

3. Defaultable Bonds

In this section, we derive the pricing formulae for defaultable bonds. First, we present the general formulae for bonds in the primary–secondary framework that fulfills the assumptions in the previous sections. Under the equivalent martingale measure P , the prices of the default-free and defaultable bonds are, respectively, given by
p ( t , T ) = E [ exp ( t T r s d s ) | F t ] ,
V i ( t , T ) = β i p ( t , T ) + 1 { τ i > t } ( 1 β i ) E [ exp ( t T ( r s + λ s i ) d s ) | F t ] ,
where β i is the recovery rate of defaultable bond i, and t T < T * . Here, to reduce derivation complexity, we assume the recovery rate to be a constant for obtaining a simple martingale pricing formula. Obviously, due to (10), the prices of risky bonds consist of two parts, i.e., the recovery rate β i and the rest factor 1 β i , both discounted to t. The first part could be calculated explicitly, while the second one could be obtained in the absence of default. Moreover, the second part, being discounted at the adjusted spot rate, presents the default risk in the model. Therefore, the pricing under P depends on the expectation in (10).
Then, in the following, we assume firm A and B to be the primary and secondary firm, respectively. The default of firm A is independent of the default risk of others, while that of firm B depends on the default risk of firm A. Both firm A and B’s default risk depends on the interest rate, which could be owing to the firms’ big cash demands for debt repayments. In this work, we only consider the simple case, with the default intensities following the linear relations:
λ t A = b 0 A + b 1 A r t ,
λ t B = b 0 B + b 1 B r t + b 1 { τ A t } ,
where b 0 A , b 1 A , b 0 B , b 1 B , and b are positive constants. The only economy-wide state variable we account for here is the risk-free interest rate, which, more specifically, follows an extended CIR process with positive jumps:
d r t = α ( η r t ) d t + θ r t d W t + d J t ,
where we assume that α , η , θ 0 , J t = i = 1 M t Y i . M t is a Poisson process with frequency ρ , which indicates the total number of jumps up to time t. { Y i , i 1 } are the sizes of the jumps and are assumed to be independent and identically distributed random variables following a distribution function F ( x ) ( x > 0 ) . When the value of the CIR process approaches 0, the behavior of such a process relies on its parameters. There are several conditions that assure that the jump-CIR process (13) remains strictly positive. We list those conditions as follows:
  • θ 2 2 α η . This condition is needed to guarantee r t > 0 with probability 1 if r 0 > 0 , 0 is not accessible, and the process never blows up to . Since r t follows a classical CIR process, the Feller condition θ 2 2 α η assures that r t > 0 . Specifically, the probability of no jumps on a given time interval is positive, which also strengthens the importance of θ 2 2 α η . If θ 2 > 2 α η , the process r t will repeatedly cross zero and reflect at zero.
  • F ( Y i ) F t F ( Y i ) is the distribution of Y i , which ensures the Markov property of the r t . That is, if one can predict the exact jump timing and jump-size at an instant before a jump occurs, one can make an arbitrarily large profit with certainty.
  • E Y i [ 0 , ) , r t < Y i < at M t = 1 , which preserves non-negative jump intensity and a positive instantaneous volatility. For a more specific description about the conditions where the jump-CIR process remains strictly positive, one can refer to [30].
Comparing with the Vasicek model, the CIR process possesses the properties of mean-reverting and non-negativity, which is of better modeling for interest rates. When η = θ = 0 , then JCIR processes reduce to shot noise processes, which are widely used in various fields, e.g., actuarial science, mathematical finance and electronics. Accordingly, investigating JCIR processes has important academic value and application value.
In the following, a critical lemma is brought out to price credit securities. We can refer to [31] for specific proof methods via martingale theory as well as piecewise deterministic Markov process theory.
Lemma 1.
Assume that μ 0 , k 0 , then the joint Laplace transform of the distribution of ( r t , R t ) (let R t = 0 t r u d u ) is
g ( μ , k , t , T ) = E { e μ r T e k ( R T R t ) | F t r } = exp ( B μ , k ( t , T ) r t ) exp α 2 η ( T t ) θ 2 ( C μ , k ( t , T ) ) 2 α η θ 2 × exp ( ρ 0 T t [ 1 h ( B μ , k ( 0 , u ) ) ] d u ) ,
where
h ( B μ , k ( 0 , u ) ) = 0 e B μ , k ( 0 , u ) x d F ( x ) ,
B μ , k ( t , T ) = ( 2 k α μ ) + μ α 2 + 2 k θ 2 c o t h ( α 2 + 2 k θ 2 ( T t ) / 2 ) ( θ 2 μ + α ) + α 2 + 2 k θ 2 c o t h ( α 2 + 2 k θ 2 ( T t ) / 2 ) ,
C μ , k ( t , T ) = c o s h ( α 2 + 2 k θ 2 ( T t ) / 2 ) + ( θ 2 μ + α ) ( α 2 + 2 k θ 2 ) 1 2 s i n h ( α 2 + 2 k θ 2 ( T t ) / 2 ) .
Proof. 
1 . (Construct a martingale) Setting the infinitesimal generator operator A to 0, we construct a martingale. Let f ( R , r , t ) = exp { G ( t ) r k R + H ( t ) } . By [32], letting A f ( R , r , t ) = 0 , we have
G ( t ) r + H ( t ) k r α ( η r ) G ( t ) + 1 2 θ 2 r G 2 ( t ) + ρ [ y ( G ( t ) ) 1 ] = 0 .
Solving Equation (18), we get
G ( t ) k + α G ( t ) + 1 2 θ 2 G 2 ( t ) = 0 , H ( t ) α η G ( t ) + ρ [ y ( G ( t ) ) 1 ] = 0 .
By solving the above coupled equations, we can get the explicit expression of G ( t ) and H ( t ) . Then for l > 0 and k 0 , when 0 t < l α 2 + 2 k θ 2 , exp { G ( t ) r t k R t + ρ 0 t [ 1 y ( G ( u ) ) ] d u + α η 0 t G ( u ) d u } is a martingale, where
y ( ϕ ) = 0 e ϕ x d F ( x ) ,
and
G ( t ) = α θ 2 α 2 + + 2 k θ 2 θ 2 coth α 2 + 2 k θ 2 t l 2 .
2 . (Laplace transforms of the distributions of r t and R t ) For an arbitrary fixed time t ( 0 s t < l α 2 + 2 k θ 2 ) , by the properties of martingale, we have
E { exp { G ( t ) r t k R t + ρ 0 t [ 1 y ( G ( u ) ) ] d u + α η 0 t G ( u ) d u } | F s r } = exp { G ( s ) r s k R s + ρ 0 s [ 1 y ( G ( u ) ) ] d u + α η 0 s G ( u ) d u } .
Then,
E { exp { G ( t ) r t k ( R t R s ) } | F s r } = exp { G ( s ) r s } exp { ρ s t [ 1 y ( G ( u ) ) ] d u } exp { α η s t G ( u ) d u } .
Set G ( t ) = μ 0 . By (21), we have
l = α 2 + 2 k θ 2 t ln μ θ 2 + α α 2 + 2 k θ 2 μ θ 2 + α + α 2 + 2 k θ 2 .
From (24), we easily get that l > α 2 + 2 k θ 2 t , i.e., t < l α 2 + 2 k θ 2 . Substituting (24) in to G ( s ) , we get
G ( s ) = B μ , k ( s , t ) = ( 2 k α μ ) + μ α 2 + 2 k θ 2 c o t h ( α 2 + 2 k θ 2 ( t s ) / 2 ) ( θ 2 μ + α ) + α 2 + 2 k θ 2 c o t h ( α 2 + 2 k θ 2 ( t s ) / 2 ) ,
then, we have
E { exp { G ( t ) r t k ( R t R s ) } | F s r } = exp { B μ , k ( s , t ) r s } exp { ρ s t [ 1 y ( B μ , k ( u , t ) ) ] d u } exp { α η s t B μ , k ( u , t ) d u } .
For the right side of (26), letting v = t u and by integral calculation, we get the results of Lemma 1. □
Until now, the closed-form solution of the Laplace transforms of r t and R t are derived. Setting μ = 0 and k = 0 in (14) gives the following corollaries.
Corollary 1.
Assume that μ 0 , k 0 ; then, the Laplace transforms of the distributions of r t and R t are, respectively,
E { e μ r T | F t r } = exp ( B μ , 0 ( t , T ) r t ) exp α 2 η ( T t ) θ 2 ( C μ , 0 ( t , T ) ) 2 α η θ 2 × exp ( ρ 0 T t [ 1 h ( B μ , 0 ( 0 , u ) ) ] d u ) ,
and
E { e k ( R T R t ) | F t r } = exp ( B 0 , k ( t , T ) r t ) exp α 2 η ( T t ) θ 2 ( C 0 , k ( t , T ) ) 2 α η θ 2 × exp ( ρ 0 T t [ 1 h ( B 0 , k ( 0 , u ) ) ] d u ) .
For mathematical simplicity, the jump size { Y i , i 1 } is assumed to be exponentially distributed; in other words,
F ( x ) = 1 e ω x ( ω > 0 , x > 0 ) .
As a special case of Corollary 1, the following corollary is presented.
Corollary 2.
Assume that μ 0 , k 0 , and the distribution of jump sizes is F ( x ) = 1 e ω x ( ω > 0 , x > 0 ) ; then, the Laplace transforms of the distributions of r t and R t are
E { e μ r T | F t r } = exp ( B μ , 0 ( t , T ) r t ) exp α 2 η ( T t ) θ 2 ( C μ , 0 ( t , T ) ) 2 α η θ 2 × 2 α ( ω + μ ) exp ( α ( T t ) ) ω ( θ 2 μ + 2 α ) exp ( α ( T t ) ) + ( 2 α μ ω θ 2 μ ) 2 ρ 2 α ω θ 2 ,
and
E { e k ( R T R t ) | F t r } = exp ( B 0 , k ( t , T ) r t ) exp { H 1 ( k ) ( T t ) + H 2 ( k ) ln D k ( t , T ) H 3 ln C 0 , k ( t , T ) } .
where
D k ( t , T ) = cosh ( α 2 + 2 k θ 2 ( T t ) / 2 ) + ω α + 2 k α α 2 + 2 k θ 2 sinh ( α 2 + 2 k θ 2 ( T t ) / 2 ) ,
H 1 ( k ) = α 2 η ( T t ) θ 2 2 k ρ ω ( α 2 + 2 k θ 2 + α ) + 2 k ω ρ α 2 + 2 k θ 2 2 k + 2 ω α ω 2 θ 2 ,
H 2 ( k ) = 2 ω ρ 2 k + 2 ω α ω 2 θ 2 ,
H 3 = 2 α η θ 2 .
Thus, from (10) and Lemma 1, we can derive the pricing formulae of defaultable bonds. The time-t prices of the A- and B-bond with maturity date T will be shown in the following theorem.
Theorem 1.
In the primary–secondary framework, assume that firm A and firm B have the same maturity date with 0-valued recovery rate ( β A = 0 , β B = 0 ) . Suppose that intensity process λ t A satisfies (11) and process λ t B satisfies (12); then, the time-t prices of the bonds issued by primary firm A and secondary firm B, respectively, are
V A ( t , T ) = E [ e ( 1 + b 1 A ) ( R T R t ) | F t r ] e b 0 A ( T t ) = g ( 0 , 1 + b 1 A , t , T ) e b 0 A ( T t ) ,
and
V B ( t , T ) = e ( b 0 B + b ) ( T t ) g ( 0 , 1 + b 1 B , t , T ) + b e ( b 0 B + b ) T + ( b 0 B + b 0 A ) t × t T e ( b 0 A b ) s · exp α 2 η ( T s ) θ 2 ( C 0 , 1 + b 1 B ( s , T ) ) 2 α η θ 2 × exp ( ρ 0 T s [ 1 h ( B 0 , 1 + b 1 B ( 0 , u ) ) ] d u ) × g ( B 0 , 1 + b 1 B ( s , T ) , 1 + b 1 B + b 1 A , t , s ) d s .
Proof. 
Firstly, from (11) and Lemma 1, we can easily illustrate that (31) holds. Secondly, according to (10), (12), and the property of conditional expectation, we obtain the time-t price of the bond with maturity date T issued by firm B
V B ( t , T ) = E [ exp ( t T ( r s + λ s B ) d s ) | F t ] = E [ exp ( b 0 B ( T t ) ( 1 + b 1 B ) ( R T R t ) b ( T τ A ) 1 { τ A T } ) | F t ] = E [ exp ( b 0 B ( T t ) ( 1 + b 1 B ) ( R T R t ) ) × E [ exp ( b ( T τ A ) 1 { τ A T } ) | F t F T * r ] | F t ] .
To calculate the time-t price of the B-bond with maturity T, we need to evaluate the conditional expectation embedded in (33). According to (2) and the law of integration by parts, we obtain that
E [ exp ( b ( T τ A ) 1 { τ A T } ) | F t F T * r ] = ( t T + T ) e b ( T s ) 1 { s T } d ( 1 e b 0 A ( s t ) b 1 A ( R s R t ) ) = e b ( T t ) ( 1 + b t T e ( b 0 A b ) ( s t ) b 1 A ( R s R t ) d s ) .
Hence, according to the expectation obtained, we get
V B ( t , T ) = e ( b 0 B + b ) ( T t ) E [ e ( 1 + b 1 B ) ( R T R t ) | F t ] + b e ( b 0 B + b ) T + ( b 0 B + b 0 A ) t × t T e ( b 0 A b ) s E [ e ( 1 + b 1 B + b 1 A ) ( R s R t ) ( 1 + b 1 B ) ( R T R s ) | F t ] d s .
For the first expectation term, we easily obtain
E [ e ( 1 + b 1 B ) ( R T R t ) | F t ] = g ( 0 , 1 + b 1 B , t , T ) .
Furthermore, using the properties of iterated conditional expectation, we get
E [ e ( 1 + b 1 B + b 1 A ) ( R s R t ) ( 1 + b 1 B ) ( R T R s ) | F t ] = E [ e ( 1 + b 1 B + b 1 A ) ( R s R t ) · E [ e ( 1 + b 1 B ) ( R T R s ) | F s ] | F t ] = exp α 2 η ( T s ) θ 2 ( C 0 , 1 + b 1 B ( s , T ) ) 2 α η θ 2 × exp ( ρ 0 T s [ 1 h ( B 0 , 1 + b 1 B ( 0 , u ) ) ] d u ) g ( B 0 , 1 + b 1 B ( s , T ) , 1 + b 1 B + b 1 A , t , s ) .
Now, substituting (36) and (37) into (35), (32) holds. The proof of Theorem 1 is complete. Next, we will present the numerical results of the yield spread. □
Until now, within the primary–secondary framework, we obtain the explicit pricing of defaultable bonds. Now, we pay attention to the yield spread, which is a function of maturity to explore the effect of counterparty risk on firm B’s bond price. In continuous-time modeling, we define the yield spread ψ between defaultable bond v and default-free bond p as
ψ ( t , T ) = 1 T t ln v ( t , T ) p ( t , T ) .
In the following, we present an example to discuss the term structure of yield spreads between primary firm A and secondary firm B.
Example 1.
Assume that β A = β B = 0 , then the yield spread of firm B’s bond is the corresponding martingale default probability averaged over T t . In the economy, if we ignore the counterparty risk and the short-term interest rate, mathematically speaking, b = 0 and b 1 B = 0 ; then, the yield spread will be the constant b 0 B . When b > 0 ( b 1 B = 0 ), firm B’s default is “accelerated" by firm A’s default with an increasing pattern in the yield spread. In contrast, a “deceleration" holds when b < 0 ( b 1 B = 0 ), with the yield spread being a decreasing function of maturity. These effects are illustrated in Figure 1 with some given parameters.
In contrast with the Markov model based on the ratings of Jarrow et al. [33], in our model, we ignore the interaction between market risk and credit risk. Usually investment-grade issuers’ credit spread curves are upward-sloping, while those of speculative-level issuers are tilted downward. However, when the possibility of default is affected by economy-wide risk factors, there will be additional aspects to this problem, as illustrated in Figure 1. Except for Panel (a), obviously, the shape of the default-free term structure and the coefficient b 1 A jointly determine the term structure of the credit spread.
When b 1 B 0 , r t follows (13); here, we assume that the jump’s sizes follow an exponential distribution with parameter ω > 0 . Since jump process is involved in our model, Equation (32) is far more complex for computing the yield spread compared to that of Figure 1. Futhermore, the shape of a contemporaneous default-free term structure also affects the yield spread. However, we expect the patterns exhibited in Figure 1 to persist. Namely, the yield spread will widen (narrow) over time with the default risk from the long (short) position.
In Figure 2, this results come from both economy-wide state variables and counterparty risk. We find that jump diffusion risk and counterparty risk have caused significant changes in the effects of different b 1 B on yield spread, while the CIR process dilutes the effect of b 1 B on yield spread, so the impact of b 1 B on yield spread is minimal. Besides these conjectures, in Figure 3, we could also get a way to extract information on the counterparty risk with empirical evidence. The effect of counterparty risk has caused the distortion of the term structure of the primary firms’ credit spreads (see Panel (d), (f) of Figure 4), which is basically the re-adjusted default-free term structure. In fact, during estimation, counterparty risk causes the derivation between the credit spread curve and the shape of the default-free term structure. This mainly comes from the simple one-factor term structure model as well as the default process, which only depends on the contemporaneous level of the spot rate. It is easy to extend our model to a multi-factor environment where other macro-economic factors can also affect default.
In Figure 5, we also get the way to extract information on the economy-wide risk with empirical evidence. The secondary firms’ possibility of default is affected by economy-wide risk factors, as is illustrated in Figure 4. Compared with Figure 4, interest rate is driven by the jump-diffusion Cox–Ingersoll–Ross (JCIR) process simultaneously in Figure 5. In addition, the default intensity is also correlated to the default states of counterparty firms. From Figure 5, we find that the larger the b 1 B , the larger the short-term yield spread; b 1 B is very sensitive to yield spread, and this phenomenon is completely in line with the facts.

4. Default Swap Rates

In this section, based on the results in Section 3, we derive the swap rate in the reduced-form model. We assume that firm C holds bonds from reference firm A, and in the meantime, is a CDS buyer from firm B for a compensation of the possible loss induced by firm A’s default. The maturity dates of the bond and the CDS are assumed to be T 1 and T 2 ( T 2 T 1 ), respectively. Here, we only consider firm A and B as being the primary and secondary party, respectively. The recovery rate of the bond issued by firm A is assumed to be zero, with its face value being 1 dollar. The swap rate could be obtained by the following theorem.
Theorem 2.
Assume that r t follows the stochastic differential Equation (13), while λ A and λ B satisfy (11) and (12), respectively. Then, the swap rate
c = V B ( 0 , T 2 ) e ( b 0 B + b 0 A ) T 2 g ( 0 , 1 + b 1 B + b 1 A , 0 , T 2 ) 0 T 2 V C ( 0 , s ) d s ,
where
g ( μ , k , t , T ) and V B ( 0 , T 2 ) , V C ( 0 , s ) are given by Lemma 1 and Theorem 1, respectively.
Proof. 
First, the time-0 market value of buyer C’s payment to B is
E [ 0 T 2 exp ( 0 s r u d u ) c 1 { τ C > s } d s ] = c 0 T 2 V C ( 0 , s ) d s ,
while that of B’s promised compensation in the case of A’s default is
E [ 1 { τ A T 2 } exp ( 0 T 2 r u d u ) 1 { τ B > T 2 } ] .
Thus, using the principle of no-arbitrage, we obtain
c = E [ exp ( 0 T 2 ( r u + λ u B ) d u ) ] E [ 1 { τ A > T 2 } exp ( 0 T 2 ( r u + λ u B ) d u ) ] 0 T 2 V C ( 0 , s ) d s .
A brief inspection of the expression in (42) confirms two intuitions. First, c is equal to zero when the reference asset A is default-free. Second, the swap rate increases with the default probabilities of A and C, and decreases with that of B. For more details, one can refer to Appendix F of [7]. Note that 0 T 2 V C ( 0 , s ) d s and the first term of the numerator of (42), which is equal to V B ( 0 , T 2 ) , can be derived by (31) and (32), respectively. Substituting (12) into the expectation term above, we obtain
E [ 1 { τ A > T 2 } exp ( 0 T 2 ( r u + λ u B ) d u ) ] = E [ 1 { τ A > T 2 } exp ( 0 T 2 ( b 0 B + ( 1 + b 1 B ) r u + b 1 { τ A u } ) d u ) ] = E [ E [ 1 { τ A > T 2 } | F T * ] exp ( b 0 B T 2 ( 1 + b 1 B ) R T 2 ) ] = E [ exp ( b 0 B T 2 ( 1 + b 1 B ) R T 2 0 T 2 λ u A d u ) ] ,
where the second equality involves the property of conditional expectation, the last one follows (2). Substituting (11) into (43), we obtain
E [ exp ( b 0 B T 2 ( 1 + b 1 B ) R T 2 0 T 2 λ u A d u ) ] = exp ( ( b 0 B + b 0 A ) T 2 ) g ( 0 , 1 + b 1 A + b 1 B , 0 , T 2 ) .
Now, substituting (44) into (42), (39) holds. The proof of Theorem 2 is complete. □
In conjunction with the pricing results of defaultable bonds, we obtain the explicit formula of default swap rates, which is an important indicator in the risk management market. With given parameters, we can obtain the relevant prices of defaultable bonds and swap rates. However, the method for evaluating the parameters is beyond the scope of our article and is left for future discussion. For an insight into the improvement of the results in our model, the numerical illustrations for swap rate are provided.
Example 2.
We assume that the jump’s sizes follow an exponential distribution with parameter ω > 0 , and firm C is not defaulted in the past with V C ( 0 , s ) = E [ e R s ] = g ( 0 , 1 , 0 , s ) . We analyze the dynamic relations between the default swap rate c and the maturity date T ( T 2 in Theorem 2) for different forms of risk-free interest rate, different mean positive jumps ω, and different counterparty risk levels b. The swap rate decreases with the maturity increasing by intuition.
Firstly, in panel (a) of Figure 6, in terms of the influence of risk-free interest rate forms chosen, although the swap rate decreases with the same trends, consistent with the intuition, there is a clear difference between the results of our model and the others. The swap rate, compared with applying the JCIR process, is overestimated in the other cases. These results clarify the essential modification of the swap rate when applying JCIR process for the risk-free interest rate.
Secondly, in panel (b) of Figure 6, we find that the default swap rate is notably affected by ω. For larger ω, generally, there will be a higher jump probability in the interest rate process, leading to a more unstable interest rate, which increases the firms’ default probability. Meanwhile, the swap rate increases with the default probabilities of firm A, but decreases with that of firm B. Moreover, firm B’s default intensities, empirically, play a decisive role in pricing the swap rate. Consequently, as a function of the maturity date, the larger the ω, the slower c decreases, which agrees with the trend obtained in panel (b). This emphasizes the necessity of considering the jump factor in the stochastic interest rate process.
Thirdly, according to the results of different b in panel (c) of Figure 6, we find that the default swap rate is significantly overpriced when ignoring the default correlation ( b = 0 ). Therefore, it is vital to consider the contagion. A successively higher b indicates an increasing default correlation between the swap seller, B, and the reference asset A. As b becomes very large, we expect the fair swap rate to drop to zero since the default swap is ineffective.
Finally, according to the results of different b 1 B in panel (d) of Figure 6, we discover that when the risk-free interest rate follows the same JCIR process, the default swap rate is affected by b 1 B to a very small extent. From a numerical perspective, in a short period of time, the larger the b 1 B value, the greater the rate of decline of the default swap rate.

5. Conclusions

In this paper, we derive the joint Laplace transform for the distribution of the vector process ( r t , R t ) , and using the joint Laplace transform, we obtain explicit pricing formulae of defaultable securities in the primary–secondary framework, where the counterparty relationship among these firms and economy-wide state variables (risk-free interest rate) follow the JCIR process, being able to influence their default probabilities. Such configuration is found to be vital for a better approach to the CDS market. We show that the interaction of market-wide risk factors with company-specific counterparty risks lead to various shapes of term structure for credit spreads. Within such a framework, the pricing of the CDS is shown. In the end, numerical graphics vividly illustrate the effect of exponential distributed jump sizes on the yield spread and CDS.
In addition, it is also of interest and is challenging to further explore other heavy-tailed distributions of the jump sizes, for instance, the Pareto distribution, the Gumbel distribution, and the Féchet distribution. However, the expression of the joint Laplace transform might be implicit.

Author Contributions

Methodology, W.L.; computation, L.T. All authors have read and agreed to the published version of the manuscript.

Funding

The work was sponsored by the Fundamental Research Funds for the Central Universities NO. 2232021D-29, NO. 2232021D-31 and NO. 2232023D-22.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

The authors thank the editor and the referees for their valuable comments and suggestions, which greatly improved the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jarrow, R.A.; Turnbull, S.M. Pricing options on derivative securities subject to credit risk. J. Finance 1995, 50, 53–86. [Google Scholar] [CrossRef]
  2. Duffie, D.; Kenneth, S. Modeling term structures of defaultable bonds. Rev. Financ. Stud. 1999, 12, 687–720. [Google Scholar] [CrossRef] [Green Version]
  3. Hull, J.C.; White, A. Valuing credit default swaps I: No counterparty default risk. J. Deriv. 2000, 8, 29–40. [Google Scholar] [CrossRef] [Green Version]
  4. Merton, R.C. On the pricing of corporate debt: The risk structure of interest rates. J. Finance 1974, 29, 449–470. [Google Scholar]
  5. Black, F.; Cox, J. Valuing corporate securities: Some effects of bond indentyre. J. Finance 1976, 31, 351–367. [Google Scholar] [CrossRef]
  6. Duffie, D.; Gârleanu, N. Risk and valuation of collateralized debt obligations. Financ. Anal. J. 2001, 57, 41–59. [Google Scholar] [CrossRef]
  7. Jarrow, R.A.; Yu, F. Counterparty risk and the pricing of defaultable securities. J. Finance 2001, 56, 1765–1799. [Google Scholar] [CrossRef]
  8. Giesecke, K.; Goldberg, L. Sequential defaults and incomplete information. J. Risk 2004, 7, 1–26. [Google Scholar] [CrossRef]
  9. Brigo, D.; Pallavicini, A.; Torresetti, R. Calibration of CDO tranches with the dynamical generalized Poission loss model. Risk 2007, 20, 70–75. [Google Scholar]
  10. Giesecke, K.; Goldberg, L.; Ding, X. A top down approach to multi-name credit. Oper. Res. 2011, 59, 283–300. [Google Scholar] [CrossRef]
  11. Ding, X.W.; Giesecke, K.; Tomecek, P.I. Time-changed birth processes and multiname credit derivatives. Oper. Res. 2009, 57, 990–1005. [Google Scholar] [CrossRef] [Green Version]
  12. Davis, M.; Lo, V. Infectious defaults. Quant. Financ. 2001, 1, 382–387. [Google Scholar] [CrossRef]
  13. Dong, Y.; Yuen, K.C.; Wang, G.; Wu, C. A reduced-form model for correlated defaults with regime-switching shot noise intensities. Methodol. Comput. Appl. Probab. 2016, 18, 459–486. [Google Scholar] [CrossRef]
  14. Chen, Y.; Xing, Y. Credit default swap pricing with counterparty risk in a reduced form model with a common jump process. Probab. Eng. Inf. Sci. 2022, 37, 275–293. [Google Scholar] [CrossRef]
  15. Dong, Y.; Yuen, K.C.; Wang, G. Valuation of CDS counterparty risk under a reduced-form model with regime-switching shot noise default intensities. Front. Math. China 2017, 12, 1085–1112. [Google Scholar] [CrossRef]
  16. Guo, J.; Wang, G. A reduced-form model with default intensities containing contagion and regime-switching Vasicek processes. Front. Math. China 2018, 13, 535–554. [Google Scholar] [CrossRef]
  17. Zhi, K.; Qian, X.; Xie, A. Valuation of kth-to-default credit-linked notes with counterparty risk in a reduced-form model. Commun. Stat.-Theory Methods 2021, 4, 1–24. [Google Scholar] [CrossRef]
  18. Liu, Y.; Hao, R.; Wang, Z. Pricing loan CDS with Vasicek interest rate under the contagious model. J. Math. Econ. 2016, 3, 416–430. [Google Scholar] [CrossRef] [Green Version]
  19. Vasicek, O. An equilibrium characterization of the term structure. J. Financ. Econ. 1977, 5, 177–188. [Google Scholar] [CrossRef]
  20. Leung, S.Y.; Kwork, Y.K. Credit default swap valuation with counterparty risk. Kyoto Econ. Rev. 2005, 74, 25–45. [Google Scholar]
  21. Li, W.; Guo, H.; Liang, J.; Brahim, A.B. The valuation of multi-counterparties CDS with credit rating migration. Appl. Math. Ser. B 2020, 35, 379–391. [Google Scholar] [CrossRef]
  22. Cox, J.; Ingersoll, J.; Ross, S. A theory of term structure of interest rates. Econometrica 1985, 53, 385–407. [Google Scholar] [CrossRef]
  23. Heston, S.L. A closed form solution for options with stochastic volatility with applications to bond and currency options. Rev. Financ. Stud. 1993, 6, 327–343. [Google Scholar] [CrossRef] [Green Version]
  24. Trutnau, G. Pathwise uniqueness of the squared Bessel and CIR processes with skew reflection on a deterministic time dependent curve. Stoch. Proc. Appl. 2011, 121, 1845–1863. [Google Scholar] [CrossRef] [Green Version]
  25. Bao, J.H.; Yuan, C.G. Long-term behavior of stochastic interest rate models with jumps and memory. Insur. Math. Econ. 2013, 53, 266–272. [Google Scholar] [CrossRef]
  26. Grasselli, M. The 4/2 stochastic volatility model: A unified approach for the heston and the 3/2 model. Math. Financ. 2017, 27, 1013–1034. [Google Scholar] [CrossRef]
  27. Ballestra, L.V.; Pacelli, G.; Radi, D. Modeling CDS spreads: A comparison of some hybrid approaches. J. Empir. Financ. 2020, 57, 107–124. [Google Scholar] [CrossRef]
  28. Artzner, P.; Delbaen, F. Default risk insurance and incomplete markets. Math. Financ. 1995, 5, 187–195. [Google Scholar] [CrossRef]
  29. Davis, M.H. Piecewise-deterministic Markov processes: A general class of non-difiusion stochastic models. J. R. Statist. Soc. Ser. B 1984, 46, 353–388. [Google Scholar] [CrossRef]
  30. Zhou, H. Jump-diffusion term structure and Itô conditional moment generator. SSRN Electron. J. 2001, 36. [Google Scholar] [CrossRef] [Green Version]
  31. Wu, Y.F. Jump-diffusion CIR model and its applications in credit risk. Hacet. Bull. Nat. Sci. Eng. Ser. B 2014, 43, 1095–1106. [Google Scholar] [CrossRef]
  32. Jacobsen, M. Point Process Theory and Applications: Marked Point and Piecewise Deterministic Processes; Springer Birkhäuser: Boston, MA, USA, 2006. [Google Scholar]
  33. Jarrow, R.A.; Lando, D.; Turnbull, S.M. A Markov model for the term structure of credit risk spreads. Rev. Financ. Stud. 1997, 10, 481–523. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Secondary firm B: term structure of yield spreads with b 0 B = 0.02 , b 1 B = 0 , b 0 A = 0.02 , b 1 A = 0 . This figure characterizes the term structure of firm B’s yield spread under different parameters b = −0.1, 0, 0.2, 0.5, 5.
Figure 1. Secondary firm B: term structure of yield spreads with b 0 B = 0.02 , b 1 B = 0 , b 0 A = 0.02 , b 1 A = 0 . This figure characterizes the term structure of firm B’s yield spread under different parameters b = −0.1, 0, 0.2, 0.5, 5.
Axioms 12 00782 g001
Figure 2. Secondary firm B: term structure of yield spreads with b 0 A = 0.02 , b 1 A = 0.01 , b 0 B = 0.02 , b = 1 and different b 1 B = 0.01 , 0 , 0.01 . Panel (a) exhibits a flat default-free term structure. Panel (b) assumes the risk-free interest rate follows the CIR process in which the parameters of the CIR model are taken from Wu [31] and Jarrow and Yu [7] with α = 0.05 , η = 0.5 , θ = 0.4 . In Panel (c), its difference with Panel (b) is that interest rate is driven by a jump process, while the parameters are the same, except that ω = 2.0 , ρ = 0.5 .
Figure 2. Secondary firm B: term structure of yield spreads with b 0 A = 0.02 , b 1 A = 0.01 , b 0 B = 0.02 , b = 1 and different b 1 B = 0.01 , 0 , 0.01 . Panel (a) exhibits a flat default-free term structure. Panel (b) assumes the risk-free interest rate follows the CIR process in which the parameters of the CIR model are taken from Wu [31] and Jarrow and Yu [7] with α = 0.05 , η = 0.5 , θ = 0.4 . In Panel (c), its difference with Panel (b) is that interest rate is driven by a jump process, while the parameters are the same, except that ω = 2.0 , ρ = 0.5 .
Axioms 12 00782 g002
Figure 3. Secondary firm B: term structure of yield spreads b 0 A = 0.02 , b 1 A = 0 , b 0 B = 0.02 , b 1 B = 0.1 with different value of b = 0.05 , 0 , 0.1 . The risk-free interest rate is assumed to follow the JCIR process in which parameters are chosen from [7,31] with α = 0.05 , η = 0.5 , θ = 0.4 , ω = 2.0 , ρ = 0.5 .
Figure 3. Secondary firm B: term structure of yield spreads b 0 A = 0.02 , b 1 A = 0 , b 0 B = 0.02 , b 1 B = 0.1 with different value of b = 0.05 , 0 , 0.1 . The risk-free interest rate is assumed to follow the JCIR process in which parameters are chosen from [7,31] with α = 0.05 , η = 0.5 , θ = 0.4 , ω = 2.0 , ρ = 0.5 .
Axioms 12 00782 g003
Figure 4. Primary firm A: term structure of yield spreads with initial interest rate r 0 = 0.05 , b 0 A = 0.02 and five different shapes of r t . The default intensity of firm A is presented in (11), where b 1 A is chosen to be different values {0.1, 0, −0.1}. Panel (a) exhibits a flat default-free term structure. In Panel (b), the interest rate is assumed to be a deterministic function of t with slope 0.2 % . In Panel (c), the interest rate is assumed to be a deterministic function of t with slope 0.1 % . Panel (d) assumes that the risk-free interest rate follows the CIR process in which the parameters of the CIR model are taken from Wu [31] and Jarrow and Yu [7] with α = 0.05 , η = 0.5 , θ = 0.4 . In Panel (e), its difference with Panel (d) is that interest rate is driven by a jump process, while the parameters are same, except that ω = 2.0 , ρ = 0.5 .
Figure 4. Primary firm A: term structure of yield spreads with initial interest rate r 0 = 0.05 , b 0 A = 0.02 and five different shapes of r t . The default intensity of firm A is presented in (11), where b 1 A is chosen to be different values {0.1, 0, −0.1}. Panel (a) exhibits a flat default-free term structure. In Panel (b), the interest rate is assumed to be a deterministic function of t with slope 0.2 % . In Panel (c), the interest rate is assumed to be a deterministic function of t with slope 0.1 % . Panel (d) assumes that the risk-free interest rate follows the CIR process in which the parameters of the CIR model are taken from Wu [31] and Jarrow and Yu [7] with α = 0.05 , η = 0.5 , θ = 0.4 . In Panel (e), its difference with Panel (d) is that interest rate is driven by a jump process, while the parameters are same, except that ω = 2.0 , ρ = 0.5 .
Axioms 12 00782 g004
Figure 5. Secondary firm B: term structure of yield spreads with b 0 A = 0.02 , b 1 A = 0 , b 0 B = 0.02 , b = 0.1 and a different value of b 1 B = 0.01 , 0.05 , 0.1 . The risk-free interest rate is assumed to follow the JCIR process in which the parameters for the JCIR model are chosen from [7,31] with α = 0.05 , η = 0.5 , θ = 0.4 , ω = 2.0 , ρ = 0.5 .
Figure 5. Secondary firm B: term structure of yield spreads with b 0 A = 0.02 , b 1 A = 0 , b 0 B = 0.02 , b = 0.1 and a different value of b 1 B = 0.01 , 0.05 , 0.1 . The risk-free interest rate is assumed to follow the JCIR process in which the parameters for the JCIR model are chosen from [7,31] with α = 0.05 , η = 0.5 , θ = 0.4 , ω = 2.0 , ρ = 0.5 .
Axioms 12 00782 g005
Figure 6. Credit default swap is correlated with the spot interest rate and counterparty risk. For the default swap rate as a function of the maturity date, see the text for detailed parameters. Panel (a): results of the models with three forms of risk-free interest rates as indicated in the panel at counterparty risk b=1.0. Panel (b): results of our model with different means of positive jumps ω at b=1.0. Panel (c): results of our model at different levels of counterparty risk. Panel (d): results of our model with different b 1 B . The basic parameters we used are α = 0.05 , η = 0.5 , θ = 0.4 , ω = 2.0 , ρ = 0.5 , r 0 = 0.05 , b 0 A = b 0 B = 0.02 , b 1 A = b 1 B = 0.01 , b = 1.0 , which follow [7,31]; ω , b, b 1 B are chosen to be different values in (bd).
Figure 6. Credit default swap is correlated with the spot interest rate and counterparty risk. For the default swap rate as a function of the maturity date, see the text for detailed parameters. Panel (a): results of the models with three forms of risk-free interest rates as indicated in the panel at counterparty risk b=1.0. Panel (b): results of our model with different means of positive jumps ω at b=1.0. Panel (c): results of our model at different levels of counterparty risk. Panel (d): results of our model with different b 1 B . The basic parameters we used are α = 0.05 , η = 0.5 , θ = 0.4 , ω = 2.0 , ρ = 0.5 , r 0 = 0.05 , b 0 A = b 0 B = 0.02 , b 1 A = b 1 B = 0.01 , b = 1.0 , which follow [7,31]; ω , b, b 1 B are chosen to be different values in (bd).
Axioms 12 00782 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lv, W.; Tian, L. Pricing of Credit Risk Derivatives with Stochastic Interest Rate. Axioms 2023, 12, 782. https://doi.org/10.3390/axioms12080782

AMA Style

Lv W, Tian L. Pricing of Credit Risk Derivatives with Stochastic Interest Rate. Axioms. 2023; 12(8):782. https://doi.org/10.3390/axioms12080782

Chicago/Turabian Style

Lv, Wujun, and Linlin Tian. 2023. "Pricing of Credit Risk Derivatives with Stochastic Interest Rate" Axioms 12, no. 8: 782. https://doi.org/10.3390/axioms12080782

APA Style

Lv, W., & Tian, L. (2023). Pricing of Credit Risk Derivatives with Stochastic Interest Rate. Axioms, 12(8), 782. https://doi.org/10.3390/axioms12080782

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop