1. Introduction
The migration process is the evolution of the credit quality of corporate bonds, corporate liabilities, corporate loans, etc. Credit rating models in credit risk management have flourished considerably in recent years.
These models known also as migration models represent the evolution of the rating of a company or a state. They do so in order to evaluate the default risk of a bond or a loan, or the term structure of both and constitutes an important issue for risk management and pricing. The commercial rating agencies such as Moody’s Investor Service and Standard and Poor’s et al. determine the credit classes or sometimes it is done internally. The modeling of the migration process is an important issue for risk management and pricing. Rating transition matrices are of particular interest for determining the economic capital figures like expected loss and VAR for credit portfolios, but also can be helpful as it comes to the pricing of more complex products in the credit industry.
The modeling of the evolution of credit migration started with a simple Markov chain model either in discrete or continuous time. Carty and Fons [
1], with data spanning from 1976 to 1993 from Moody’s Investor service, showed that the Weibull distribution was the appropriate one to represent the time spent in a credit class. Thus, in fact, they established using real data that the Markov chain model was not realistic. A generation of models followed dominated by duration analysis, and an important study in this area which could act also as a link to these class of models is Duffie et al. [
2].
More recently, inhomogeneous semi-Markov processes have been proposed for the migration process as more realistic models for the variability of the rating transition matrices by Vasileiou and Vassiliou [
3]. D’ Amico et al. [
4] proposed a homogeneous Markov renewal model for the evolution of the credit migration process. The difference of the semi-Markov approach as defined in Howard [
5] and the Markov renewal approach is basically on the basic parameters chosen to build the model. D’ Amico et al. [
6,
7] introduced non-homogeneous semi-Markov of reliability theory type models; in D’ Amico [
8], a semi-Markov maintenance model with imperfect repair at random times was introduced; in D’ Amico et al. [
9,
10], the backward and forward non-homogeneous semi-Markov process was introduced and studied to some extent.; in Vassiliou and Vasileiou [
11], an inhomogeneous semi-Markov process is introduced to study the asymptotic behavior of the survival probabilities; in Vassiliou [
12], fuzzy sets are introduced to face the problem that rating agencies disagree with the majority of their ratings. The new stochastic and mathematical problems created with the introduction of fuzzy sets are being answered mainly for the quasi-stationarity problem. The above publications of non-homogeneous semi-Markov and Markov renewal processes for credit risk was followed by quite a lot of literature by the same and other authors in credit risk and related subjects. For example, from recent years, see Huang [
13], D’ Amico et al. [
14,
15], D’ Amico [
16], Magni et al. [
17], Wu et al. [
18,
19], Puneet et al. [
20], De Blasis [
21]. In D’ Amico et al. [
22], bivariate semi-Markov processes are introduced for the pricing of Credit Default Swaps (CDS). The Credit Default Swap (CDS) is a bilateral agreement that transfers credit risk between two contractual parties, protection buyer (that faces credit risk from a third party), and protection seller. Finally, Vassiliou [
23] introduced the idea of the stochastic Market environment to express the changes in the general economy, which affect any industry in small or great amounts of turbulence.
The vast majority of the existing models for the migration process use a discrete time setting. The roots of this tendency is the fact that the most prominent risk management tools, such as J.P. Morgan’s Credit Metrics and McKinsey’s Credit Portfolio View, are built around estimates of rating migration probabilities. In essence, the estimates by these agencies and in the published academic literature use a discrete-time setting and rely on a “cohort” method that estimates the transition rates. In what follows, our semi-Markov chain and the related Markov renewal chain will be in a discrete-time setting. Lando and Skodeberg [
24] argued that the estimation of the transition probabilities in their proposed continuous time setting among other advantages is getting a better grip on the rare events. However, the results in the present are independent from the method used to estimate the transition probabilities that is the discrete time method used in Vasileiou et al. [
3] and the continuous time method proposed by Lando and Skodeberg [
24]. Note also that the continuous time method of estimation of the transition probabilities is easily transferred from the Markov chain model of Lando and Skodeberg [
24] in the present non-homogeneous semi-Markov chain and Markov renewal models while retaining the apparent advantages of the later models. The general problem of obtaining the results that follow in a continuous time setting altogether is an interesting and challenging mathematical problem for future research probably in the near future.
In
Section 2, we provide what is already known but necessary in what follows. The Market
of the savings account, the default free zero-coupon bond, and the defaultable zero-coupon bond are defined. We also introduce the Randon–Nikodym derivatives under the real probability, the equivalent martingale, and the forward martingale measures. Finally, the need to study the change to forward the martingale measure is shown for the evolution of defaultable bonds in the various grades. This is done by providing the price of defaultable bonds and their spreads as functions of the forward probabilities.
Section 3 introduces a new definition of the
-inhomogeneous semi-Markov process than the last in Vasileiou et al [
3]. The differences brought about into the new definition, although they are not apparently essential, are the ones needed in order that:
to state theorem 1, which establishes the change from the real world probabilities to forward probabilities in an inhomogeneous semi-Markov process;
to introduce the necessary definitions, theorems, and all the results and algorithms that follow in the present paper.
Section 4 introduces a new definition of the
-inhomogeneous Markov renewal process. Markov renewal processes have been used in many studies (see D’Amico et al. [
4,
6,
7,
8] as a model for the migration process using the real probabilities. In theorem 2, we state and prove that the change of measure under certain conditions retains the Markov renewal property. In addition, we provide under the same conditions functional relationships between the real world transition probabilities and the forward transition probabilities.
In
Section 5 and in Theorem 3, we provide a new martingale characterization, for the discrete time
-inhomogeneous semi-Markov process. This characterization apart from its general interest will prove to be very useful in
Section 7 and
Section 8.
In
Section 6 and in Theorem 4, we provide a new martingale characterization, for the discrete time
-inhomogeneous Markov renewal process.
We alter in
Section 7 the conditions of Theorem 1, in a way that the theorem will still be valid under the new conditions. The new conditions are more restrictive but still remain quite general. However, they will be useful and in Theorem 5 we establish new interesting simplified closed analytic relations among the forward probability measure transition probability sequences and the corresponding real world probability measure ones under the new conditions.
In
Section 8, we study the calibration of the
-inhomogeneous semi-Markov model. In
Section 8.1, we make the new risk premium assumptions and provide their consequences on the results of Theorem 5. In
Section 8.2, we introduce the real and forward entrance probability measures and establish stochastic difference equations among the real entrance probability measures. In addition, functional relationships are established between the forward entrance probability measures and the default free zero-coupon bonds and the defaultable zero-coupon bonds. In
Section 8.3, we propose Algorithm 8.1 for the evaluation of the needed forward entrance probability measures. The algorithm, although similar in nature and perception as the algorithm in Section 6 of Vasileiou and Vassiliou [
3], also serves the purpose of clarifying many details. In addition, by understanding in some depth the present Algorithm 8.1, we could go back to the algorithm in Section 6 of Vasileiou and Vassiliou [
3] and easily make some needed corrections.
In
Section 9, we provide an illustrative application of the inhomogeneous semi-Markov model in the evolution of the migration process in credit risk. The forward entrance probabilities are evaluated in a classical problem using representative data.
2. A General Discrete-Time Market Model
Let a complete filtered probability space with being the real-world probability measure which represents the actual probability of events in , the -algebra, and a filtration in where the variable t is discrete and represents the time, and is being generated by the economic assets to be modeled in the filtered probability space. Let the time set , where T is the terminal date of the economic assets, that is, the trading horizon. The points of are the admissible trading dates.
The Market. Let
be the discrete time stochastic processes which is the interest rate of the savings account in the market, and it is assumed to be adapted to the filtration.
The savings account of the markethas a value at timet given by
In addition, in the market, the default free zero-coupon bond whose price isfor is included. Let that k rating grades are distinguished by the rating agencies for a defaultable bond. Define by the price process of a defaultable bond of grate i for specific and , We call the market with the above economic assets the market , and we assume that it is perfect, i.e., all assets in market are perfectly divisible and the market is frictionless.
Let the column vector
represent the price process of the assets in market
:
By using as num
raire the savings account, we get the column vector
with
, and is called the discounted price process of the market.
Let
be an equivalent martingale measure equivalent with the real probability measure
. It is known then that
is such that the discounted bond price
, for
and all possible values of
T are a martingale under
and in relation with the filtration
Analogously, a
is an equivalent martingale measure for the market
, if it is an equivalent martingale measure for the discounted asset price process
. How to construct such equivalent martingale measures is known and could be found, for example, in Bingham and Kiesel [
25], Shreve [
26] and Vassiliou [
27].
We assume that the market
is viable and complete. Then, it is known (Vassiliou et al. [
23]) that the equivalent martingale measure
is unique for the market
and guarantees the existence of the forward martingale measure
. We have seen that
is an equivalent martingale measure for the discounted bond price
, hence it is known that
Now, let the Radon–Nikodym derivative
where the
-measurable random variable
is strictly positive
a.s. and
Then, the density process
follows a strictly positive martingale under
It is known that the forward measure
on
with
T the trading horizon or equivalently the forward martingale measure
on
is equivalent to
, and their Radon–Nikodym derivative is equal to
When restricted to the
algebra
, we get that
for every
From
and
, we find that
where the
measurable random variable
is strictly positive
-a.s. and
It is known that the density process
follows a strictly positive martingale under
Let
be the process of price at time t of a defaultable bond. Then, from Bielecki and Rutkowski [
28] or Vasileiou and Vassiliou [
3], we obtain
where
, which is equal to 1 if
and equals 0 in all other cases,
is the default time and
is
the recovery rate. Following Vasileiou and Vassiliou [
3] p. 176, we have that
It is apparent that it is important to estimate the probability of default under the forward probability measure. Among other reasons, it allows us to relax the assumption of constant interest rate and have instead an interest rate process
By definition, the forward rate of one time step of the risky bond at time
T, as seen from
, is given by
Analogously, for the default-free bond, we have
It follows that the credit spread process will be given by
Apparently, the above is true for any credit state i. We then have
From the above relations, it is apparent that it is necessary to calculate the probability . Therefore, we need to study the change to forward the martingale measure for the model we will use for the evolution of defaultable bonds in the various states.
3. The Inhomogeneous Semi-Markov Process and Change to
Forward Martingale Measure
In [
3], a different and novel definition of an inhomogeneous semi-Markov process appeared than the one in [
29]. The delicate essential differences in the cores of the two definitions were those that made possible the Proof of Theorem 4.1, p.179 in Vasileiou and Vassiliou [
3], which dealt with the change of measure from the real world probability measure to the forward probability measure in an inhomogeneous semi-Markov process. In the present section, we will amend the definition of an inhomogeneous semi-Markov process. The differences brought into the new definition, although they do not seem essential, are the ones that will allow the existence and use of Theorem 1 in the present, and the definitions, theorems, and all the results and algorithms that follow in the present paper.
Let be a stochastic process on the complete filtered probability space , with state space , where k is the the number of grades of a defaultable bond, and the default state. The value means that the defaultable bond enters state i at time t. Denote by the natural filtration of and assume . Now, let that represents the choice of movement in the interval [t,t+1) of the defaultable bond, given that it last entered its present grade at time t. Let the natural filtration of the process be denoted by and assume that V. The pair is a discrete-time -inhomogeneous semi-Markov process if what follows holds:
(
a) The process
is a non-homogeneous Markov chain adapted to
with respect to
under
, that is, if
is any function, then
Denote by
and
the transition probabilities with
(
b) According to (a), whenever a defaultable bond enters state i at time t, it chooses credit class j, and moves to it with
. Now, before performing the actual transition from credit class i to credit class j, it “holds” for a time
in state i. All
are positive, integer-valued random variables with a probability mass function
, where
The matrix
is of the form
This means that all elements of are assumed finite and .
(
c) Let
be any function, then
where
.
We call the inherent -inhomogeneous Markov process.
Our next goal is to examine the changes brought into Theorem 4.1 p. 179 in Vasileiou and Vassiliou [
3] by the changes, made in the present, in the definition of an inhomogeneous semi-Markov processes. Following the steps of the Proof of Theorem 4.1 in Vasileiou and Vassiliou [
3], it is possible to show that, as expected, the preservation of the
-inhomogeneous semi-Markov property remains, the relationships between the real-world probabilities
and probability mass functions
to the respective probabilities under
, are essentially the same, but the conditions under which both results hold changes accordingly. These results are stated in the following theorem, the proof of which is omitted due to the similarity with Theorem 4.1 p.179 in Vasileiou and Vassiliou [
3]:
Theorem 1. We assume the random variable be -measurable for any Consequently, for all values of t: for any function . In addition, we assume that is -measurable for any Equivalently, for some function . Assuming that follows a discrete-time -inhomogeneous semi-Markov process under , then a discrete-time -inhomogeneous semi-Markov process under also follows, and, in addition, 4. The Non-Homogeneous Markov Renewal Process and Change
to Forward Martingale Measure
The above definition of the
-non-homogeneous semi-Markov process was based on the definition of the non-homogeneous semi-Markov process given in Vassiliou and Papadopoulou (1992), which had its roots on the definition of the homogeneous semi-Markov process given in Howard [
5]. Pyke [
30] provided the definition and preliminary properties of the homogeneous Markov renewal process. The two stochastic processes are equivalent; however, the basic parameters defining the processes are quite different. Many authors have used homogeneous and inhomogeneous Markov renewal processes for the immigration process of defaultable bonds (see D’Amico et al. [
4,
6,
7,
8]). Therefore, there is a need to define an
-inhomogeneous Markov renewal process and study in detail the change to forward martingale measure of their founding parameters.
Let be a complete filtered probability space, and a stochastic process taking values in the state space and, representing the credit state, the defaultable bond enters at the n-th jump. Let the natural filtration of and assume that it is a subfiltration of . Let be the family of random variables taking values in such that expressing the time of the n-th jump . Let be the natural filtration of and allow that .
The stochastic process
is said to be a non-homogeneous Markov renewal process provided that, for every bounded or non-negative function
, we have
for
. Define by
The family of probabilities for , , is called a non-homogeneous Markov kernel.
The function
has all the properties of a distribution function except that
Indeed, one could easily find that
If
for some pair
, then
for all
x; we then define
. With this convention, we define
From
and
we get that
Remark 1. Note that, in words, the increments and in the inhomogeneous case, unlike the homogeneous, are not conditionally independent, given the Markov chain and the jumping times.
We will need the following Lemma from Musiela and Rutkowski [
31].
Lemma 1. (Abstract version of Bayes Formula). Let be a sub σ-field of the σ-field . Now, consider ψ a random variable integrable in relation to ; then, is a probability measure equivalent to withwith the random variable strictly positive -a.s.; moreover, is integrable with Then, With the use of the above Lemma, we prove the following basic theorem:
Theorem 2. Let be a complete filtered probability space and a non- homogeneous Markov renewal process with Markov kernel given by . In addition, let the random variable be measurable for all . Then, for every , The pair follows a -non-homogeneous Markov renewal process under and moreover where is the transition probability matrix of the forward probabilities.
Proof. Using Lemma 1, we fix
and, for any
, we get
Now,
actually shows that
is a
-measurable random variable. Since
, we get that
□
5. Martingale Characterization of the -Inhomogeneous Semi-Markov Process
Martingale characterization theorems have played traditionally a basic role in the theory of stochastic process—for example, in the methodology to obtain Laws of Large numbers and Central Limit theorems. In the present section, in the form of a theorem, we will provide a martingale characterization, for the discrete time
-inhomogeneous semi-Markov process
. The theorem will not be proved in the at most possible generalization, but in a form suitable for the purposes of
Section 7. In this respect, define the jump processes:
and the processes
the number of times the - inhomogeneous semi-Markov process entered state i and selected to move to state j in the time interval
the number of transitions of the -inhomogeneous semi-Markov process from i to j that occurred in the time interval remaining time m in state i.
Now, we arrive at the following relations for the above defined processes:
At this point, it is important to note that, for a discrete time non-homogeneous semi-Markov process , the sequences of transition probabilities and uniquely determine the process. In addition, together with the initial distribution, the above sequences provide sufficient information, to answer the important questions in the theory of semi-Markov processes and also the ones which proved to have great practical value in real problems.
Let be the set of all matrices with elements from In addition, let be the set of stochastic matrices. We now provide the following definition:
Definition 1. Let the sequence ; then, we say that it is an -sequence if the following conditions are satisfied:where are a matrix of ones. Theorem 3. Let there be - non-homogeneous semi-Markov process, with sequences of probabilities and . For any two states , the processesfor , following -martingales under . Conversely, letthen, the processesare -martingales under , then the sequences of matricesare equal, respectively, with the sequences of matrices In conclusion, the sequences of transition probabilitiesare the unique sequences for which the processesare -martingales under . Proof. It is sufficient to show:
We assume that
and, for
, the proof will be similar. For the first part of Equation
, we have
Now, from the second part of Equation
, we have
Conversely, since
is a martingale under
, we have that
From Equations
and
, we get
or
from which we conclude that
Since
is a martingale under
, we have that
From Equations
and
, we get
or
from which we conclude that
□
6. Martingale Characterization of the -Non-Homogeneous
Markov Renewal Process
We start by providing in the form of a theorem a martingale characterization, for the -non-homogeneous Markov renewal process . Define the process
the number of transitions of the -IMRP from i to j that occurred in the interval with sojourn time to state i being less than x.
In addition, define the jump processes
Then, we define the following stochastic process:
We will now prove the following martingale characterization theorem for a -inhomogeneous Markov renewal theorem.
Theorem 4. Let be a complete filtered probability space and a non-homogeneous Markov renewal process with Markov kernel given by . For any two states , the processfor follows -martingales under . Conversely, let the inhomogeneous Markov kernel ; then, ifis a -martingale under , then the inhomogeneous Markov kernel is equal to the inhomogeneous Markov kernel In conclusion, the inhomogeneous Markov kernel is the unique Markov kernel for which the processes are a -martingale under . Proof. It is sufficient to show:
We assume that
and for
the proof will be similar. We get that
Conversely, let
; then, we will prove that, if the process
is a
-martingale under
, then the transition probabilities
are equal to
.
Since the process
is a
-martingale under
, we have that
From
and
we get that
from which we get that
□
7. New Closed Analytic Functional Relationships between Forward and
Real World Probabilities of Transition
We will start by altering the conditions of Theorem 1, in a way in which the theorem will still be valid under the new conditions. The new conditions are more restrictive but still remain quite general. However, they will be useful in what follows, in order to provide interesting and simplified closed analytic relations, between the forward probability measure transition probabilities sequences and the corresponding real world probability measure ones.
Note that the conditions under which theorem 1 is valid were the following:
The random variable
be
-measurable for any
. As a result, we get that, for every such t, we have
for some function
.
The random variable
is
-measurable for all
Equivalently,
for some function
.
Let us now assume that the following Condition 1 holds:
Condition 1.The random variable is -measurable for any , and there exists a finite set, denoted by , such that the random variable admits the following representation: for some functions , and The random variable is -measurable for every and there exists a finite set, denoted , such that the random variable admits the following representation:for any and Remark 2.(i)It is apparent that the conditions of Theorem 1 are more general, in the sense that, when Condition 1 holds, then also that of Theorem 1 is satisfied. Conversely, the conditions of Theorem 1 could hold without the Condition 1 to be true. However, Condition 1 is quite general and the assumptions one has to make, during the calibration of the model (Vasileiou and Vassiliou [3]), are such that one would rather freely state that the two conditions are “almost” equivalent for our purposes.(ii)It is known (Vasileiou and Vassiliou [3]) that a discrete-time -homogeneous Markov chain is a discrete -time inhomogeneous semi-Markov process for which It is an immediate consequence then to check that, in the case of the -homogeneous Markov chain model, Condition 1 coincides with Condition in Bielecki and Rutkowski [28]. Theorem 5. Let be the -inhomogeneous semi-Markov process under . Suppose that Condition 1 holds. Then, we may choose as (a) of Condition 1then any chosen function should satisfy In addition, we may choose as (b) of Condition 1then any chosen function should satisfy Then, is a discrete time -inhomogeneous semi-Markov process under the forward probability measure , and, for every , we have the following: In addition, the chosen functions and should be such that
Proof. In order for Condition 1 be acceptable, it should satisfy two prerequisites. First,
should be a martingale process adapted to
and secondly
should be strictly positive. From Equation
and Theorem 5, we have
and so the process
is
-measurable, i.e., it is a function; let us say
where
From Equation
it is apparent that we may choose as
of Condition 4.1 the following:
Then, from
we get
and consequently
Thus, the stochastic process
as defined by
is a martingale under
In addition,
as defined by
should be strictly positive. It is easy to check that it is sufficient to have that
Now, from
for any
, we equivalently have
or, equivalently
from which we get
and finally
Thus, any chosen function should satisfy .
From Equation
since
, the process
is given by
Since Condition 1 holds, then the conditions of Theorem 5 are valid and thus we have
From Equation
, it is apparent that the choice of any function
apart from satisfying
, which guarantees the positivity of
, should be such that
from which we get that
Now, from Theorem 5, we have
and so the process
is a function of the form
where
It is apparent, from Equations
and
, that we may choose part b) of Condition 1 to be of the following form:
Thus, we get that
and consequently the process
is a martingale under
. In addition,
as defined by
should be strictly positive. It is easy to check that it is sufficient to have
For any
i,
j,
r, it is equivalent to have
Thus, the choice of
should be such that
Since Condition 1 holds, then the conditions of Theorem 5 are valid and thus we have
From Equation
it is apparent that any choice of a function
apart from satisfying Equation
should also satisfy the following relation:
which is equivalent with
From Equation
, since the initial information is included in the interval
, we may assume that
and then we get
It is apparent that the choices of
and
in addition to the already created restrictions, should also be such that
□
8. Model Calibration
The problem of calibration of the non-homogeneous semi-Markov model is an important one and will be resolved in the present section. The calibration firstly was discussed by Vasileiou and Vassiliou [
3] in
Section 6. In the present, we will resolve the calibration of the semi-Markov model under the new assumptions and the algorithm provided below. Although similar in nature and perception to the algorithm in
Section 6 of Vasileiou and Vassiliou [
3], it also serves the purpose of clarifying many details. In addition, by understanding in some depth the present Algorithm 8.1, we could go back in the algorithm in
Section 6 of Vasileiou and Vassiliou [
3] and easily make some needed corrections.
Available Data. Allow that is the present time. We assume that the available data span are for the time window . These data will consist of the following type:
(i) From historical data on credit migrations, we firstly estimate the matrix sequences
and
during our time window
. The main problem in estimation was that some issuers had at some time during the time window withdrawn from the rating process. This flow was treated as type III censoring (see Lee [
32]) as in biomedical data. Equivalently, it is also met as right censoring, in other studies, see for example McClean and Gribbin [
33,
34]. It is inherent that it is also assumed that the reason for the rating being withdrawn is not the understanding of possibility of a default. In the semi-Markov chain case, it is necessary to estimate the conditional density of duration in each grade and the duration in a grade before default. The same problem appears in manpower planning, where people move among the grades of an organization. For more details, see McClean et al. [
35,
36]. Given the above peculiarities of the data available and models used, the estimation methodology appeared in Subsection 6.1 of Vasileiou and Vassiliou [
3]: (ii) The observed values of
, for every
and
(iii) The observed values of , for every and
(iv) Observed values on recovery rates, , in the time window .
Our target is to estimate the sequences and , within that is, the time window. In fact, we do not need to calculate and separately as it is shown below. What is needed is to find the entrance probabilities defined below, which as it is proved in what follows are functions of the Hadamard products .
8.1. The Risk Premium Assumptions and Their Consequences
In order to calibrate the non-homogeneous semi-Markov chain models that are compatible with Theorem 5, we need to make the following assumptions, which are analogous with the case of the homogeneous Markov model.
Assumption A1.(The risk premium assumptions).
(i)
For any and , assume:(ii)
For any , and assume:The consequences of the Assumption 8.1 on the results of Theorem 7 are the following:
(a) For all
,
and for every
(b) For all
and
and for every
8.2. The Real and Forward Entrance Probability Measures
Allow the following real probability measures which we call the
entrance probabilities for the non-homogeneous semi-Markov chain:
and
. In addition, let
be the matrices of
entrance probabilities such that
where
, with
From Vasileiou and Vassiliou [
3] p.182, we get
with
for
and
with
for
From Theorem 1, we know that the change to forward measure preserves the semi-Markov structure, hence
Using
and
we get
and more generally we have
8.3. An Algorithm for Calculating the Forward
Entrance Probabilities
We are now in a position to propose an algorithm for evaluating the forward entrance probabilities which in summary is as follows:
Algorithm 8.1. (Risk Premium Algorithm).
Step 1. For all
and
using
we find
Since the forward measure preserves the semi-Markov structure using relation
we get
and, therefore, for all
and
using
and
we obtain
Now, from Theorem 5, we know that the forward measure preserves the semi-Markov process, hence, from
we get that
therefore, for all
, we find
and finally
Step 2. Using
and since, in Step 1 by relation
, we have evaluated
for all
t and all
i for
for
we obtain
From
since the change to forward measure preserves the semi-Markov process, we find that
Define by
the
-element of the matrix
, then
could be written as
and
Solve the system
and
to obtain
Hence, now we are in a position to evaluate
Since the change to forward measure preserves the semi-Markov property from
, we get
Now, from
and using
and
, we evaluate for
the entrance probabilities
. Hence, Step 2 is concluded by collecting or the evaluations done in the matrices:
We repeat the methodology indicated by Step 2 for evaluating all desired values of for
If we wish to find separate values of
and
for our time window
do as follows. Let, for every specific
t,
to be the maximum of m with
strictly positive, and define
From Theorem 5, the semi-Markov structure is preserved therefore
from which we get
and consequently
which provides the values
, and, consequently, the corresponding values
Now, from
and
, we get that
and
9. An Illustrative Application
In this section, we provide an illustrative application of the non-homogeneous semi-Markov chain in the evolution of migration process in credit risk. Credit rating classes are typically identified by a finite set of elements, which we call the set of credit rates. A good account of the rating systems of these agencies could be found in [
37]. Most rating systems involve both quantitative and qualitative information. In our semi-Markov model, the set of states are the credit grades, which, in many studies, are identified as
Following the Standard and Poor’s definitions, the above rating categories have the following meaning:
Aaa : A bond or an obligation is rated Aaa declaring that the obligor’s capacity to meet its financial commitment on the obligation is extremely strong.
Aa : A bond or an obligation is rated Aa when it differs from the highest rated obligations only by a small degree. The obligor’s capacity to meet its financial commitment on the obligation is very strong.
A : A bond or an obligation is rated A when it is more susceptible to the adverse changes in the economic conditions. The obligor’s capacity to meet its financial commitment on the obligation is still strong.
Baa : A bond or an obligation rated Baa exhibits sufficient protection. However, under adverse changes in the economic conditions, its capacity to meet its obligations will probably be weakened more.
Ba : A bond or an obligation rated Baa is considered to be a speculative issue. However, it is less vulnerable than other speculative issues.
B : A bond or an obligation rated B is more vulnerable to nonpayment than obligations rated Ba. However, currently, the obligor has the capacity to meet payments. Adverse changes in economic conditions most probably will impair that capacity.
Caa : A bond or an obligation rated Caa is currently highly vulnerable to nonpayment.
Therefore, there are seven grades and the default grade. In the present application for size purposes of the respective matrices, we will merge the grades into three and the default. If we do not do that, then we will have to provide a large number of matrices in the present application, which will take too much space. That is, the state space will be , where 4 is the default state, 1 is the grade with least risk involved and 3 the grade with the maximum risk involved. Our time window is The observed values of are provided by the market and for our present time window are the following:
Default free bonds for the time window.
The observed values of for and for our time window are the following:
Term structure of defaultable bonds for the time window.
GRADE 1:
GRADE 2:
GRADE 3:
In the above, the physical architecture of the defaultable bonds in the arbitrage free Market could be observed.
The value of is . The estimated transition probabilities for the inhomogeneous semi-Markov model available from historical data are the following:
(a) The transition probabilities sequence of the real world is
and are the following:
(b) The sequence of the real world duration probabilities is
We are interested for our purposes in finding the sequence of forward probabilities
for
and
Applying the data above in the Algorithm 8.1, we provide the values of the forward probabilities together with the relative real word entrance probabilities: