Next Article in Journal
Exchange Market Liquidity Prediction with the K-Nearest Neighbor Approach: Crypto vs. Fiat Currencies
Next Article in Special Issue
Analysis of a Batch Arrival, Batch Service Queuing-Inventory System with Processing of Inventory While on Vacation
Previous Article in Journal
Variational Nonlinear Optimization in Fluid Dynamics: The Case of a Channel Flow with Superhydrophobic Walls
Previous Article in Special Issue
Analysis of Instantaneous Feedback Queue with Heterogeneous Servers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Homogeneous Semi-Markov and Markov Renewal Processes and Change of Measure in Credit Risk

by
P.-C.G. Vassiliou
Department of Statistical Sciences, University College London, Gower Street, London WC1E 6BT, UK
Mathematics 2021, 9(1), 55; https://doi.org/10.3390/math9010055
Submission received: 18 November 2020 / Revised: 18 December 2020 / Accepted: 22 December 2020 / Published: 29 December 2020
(This article belongs to the Special Issue Stochastic Modeling and Applied Probability)

Abstract

:
For a G -inhomogeneous semi-Markov chain and G -inhomogeneous Markov renewal processes, we study the change from real probability measure into a forward probability measure. We find the values of risky bonds using the forward probabilities that the bond will not default up to maturity time for both processes. It is established in the form of a theorem that the forward probability measure does not alter the semi Markov structure. In addition, foundation of a G -inhohomogeneous Markov renewal process is done and a theorem is provided where it is proved that the Markov renewal process is maintained under the forward probability measure. We show that for an inhomogeneous semi-Markov there are martingales that characterize it. We show that the same is true for a Markov renewal processes. We discuss in depth the calibration of the G -inhomogeneous semi-Markov chain model and propose an algorithm for it. We conclude with an application for risky bonds.

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 M d 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 G -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: i to state theorem 1, which establishes the change from the real world probabilities to forward probabilities in an inhomogeneous semi-Markov process; i i 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 G -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 G -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 G -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 G -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 ( Ω , G , Q , G t ) with Q being the real-world probability measure which represents the actual probability of events in Ω , G the σ -algebra, and G t a filtration in G 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 T = 0 , 1 , 2 , , T , where T is the terminal date of the economic assets, that is, the trading horizon. The points of T are the admissible trading dates.
The Market M d . Let r t t = 0 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 G t . The savings account of the market M d has a value at timet given by
B t = u = 0 t 1 1 + r u , with B 0 = 1 ;
In addition, in the market M d , the default free zero-coupon bond whose price is B t , T for t = 0 , 1 , , T is included. Let that k rating grades are distinguished by the rating agencies for a defaultable bond. Define by D i t , T the price process of a defaultable bond of grate i for specific T T and t = 0 , 1 , , T , i = 1 , 2 , , k . We call the market with the above economic assets the market M d , and we assume that it is perfect, i.e., all assets in market M d are perfectly divisible and the market is frictionless.
Let the column vector D t , T represent the price process of the assets in market M d :
D t , T = B t 0 , B 1 t , T , D 1 2 t , T , D 2 3 t , T , , D k k + 1 t , T .
By using as num e ´ raire the savings account, we get the column vector
D ˜ t , T = 1 , β t B 1 t , T , β t D 1 2 t , T , β t D 2 3 t , T , , β t D k k + 1 t , T ,
with β t = 1 / B t 0 , and is called the discounted price process of the market.
Let Q * be an equivalent martingale measure equivalent with the real probability measure Q . It is known then that Q * is such that the discounted bond price B ˜ t , T = B t , T / B t , for 0 t T and all possible values of T are a martingale under Q * and in relation with the filtration G t . Analogously, a Q * is an equivalent martingale measure for the market M d , if it is an equivalent martingale measure for the discounted asset price process D ˜ t , T . 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 M d is viable and complete. Then, it is known (Vassiliou et al. [23]) that the equivalent martingale measure Q * is unique for the market M d and guarantees the existence of the forward martingale measure Q T . We have seen that Q * is an equivalent martingale measure for the discounted bond price B ˜ t , T , hence it is known that
B ( t , T ) = B t E Q * ( B T 1 G t ) , t 0 , T
Now, let the Radon–Nikodym derivative
d Q * d Q G T = ψ T Q a . s . ,
where the G T -measurable random variable ψ T is strictly positive Q a.s. and E Q ( ψ T ) = 1 . Then, the density process ψ t = E Q ( ψ T G t ) , t = 0 , 1 , , T , follows a strictly positive martingale under Q . It is known that the forward measure Q T on ( Ω , G T ) with T the trading horizon or equivalently the forward martingale measure Q T on ( Ω , G T ) is equivalent to Q * , and their Radon–Nikodym derivative is equal to
d Q T d Q * = B T 1 E Q * ( B T 1 ) = 1 B T B ( 0 , T ) Q * a . s .
When restricted to the σ algebra G t , we get that
η t : = d Q T d Q * G t = E Q * 1 B T B ( 0 , T ) G t = B t , T B t B ( 0 , T )
for every t 0 , T .
From 6 and 7 , we find that
d Q T d Q G T = ψ T B T B ( 0 , T ) = θ T Q a . s .
where the G T measurable random variable θ T is strictly positive Q -a.s. and E Q ( θ T ) = 1 . It is known that the density process θ t = E Q ( θ T G t ) , t = 0 , 1 , , T , follows a strictly positive martingale under Q .
Let D δ t , T 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
D δ t , T = B t E Q * B T 1 δ 1 T τ + 1 T < τ G t ,
where 1 T τ , which is equal to 1 if T τ 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
D δ t , T = B t , T E Q T δ 1 T τ + 1 T < τ G t
= B t , T δ + 1 δ Q T τ > T G t .
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 r t . By definition, the forward rate of one time step of the risky bond at time T, as seen from t T , is given by
f D t , T = log D δ t , T + 1 D δ t , T .
Analogously, for the default-free bond, we have
f t , T = log B t , T + 1 B t , T .
It follows that the credit spread process will be given by
s t , T = f D t , T f t , T = log δ + 1 δ Q T τ > T G t δ + 1 δ Q T + 1 τ > T + 1 G t .
Apparently, the above is true for any credit state i. We then have
D i δ t , T = B t , T E Q T δ 1 T τ + 1 T < τ X t = i , i = 1 , 2 , , k .
From the above relations, it is apparent that it is necessary to calculate the probability Q T τ > T X t . 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 G 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 X t t = 0 be a stochastic process on the complete filtered probability space ( Ω , G , Q , G t ) , with state space K = 1 , 2 , , k , k + 1 , where k is the the number of grades of a defaultable bond, and k + 1 the default state. The value X t = i means that the defaultable bond enters state i at time t. Denote by F X the natural filtration of X t t = 0 and assume F X G . Now, let S t + 1 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 S t t = 0 be denoted by F S and assume that F S V F X G . The pair X t , S t + 1 t = 0 is a discrete-time G -inhomogeneous semi-Markov process if what follows holds:
(a) The process S t + 1 t = 0 is a non-homogeneous Markov chain adapted to X t t = 0 with respect to G under Q , that is, if f : K R is any function, then
E Q f S t + n G t = E Q f S t + n σ X t , t , n = 0 , 1 , 2 ,
Denote by p i j t = Q S t + 1 = j X t = i and P ˜ t = p i j t i , j K the transition probabilities with
P ˜ t = p 11 t p 12 t . . . p 1 k t p 1 , k + 1 t p 21 t p 22 t . . . p 2 k t p 2 , k + 1 t p k 1 t p k 2 t . . . p k k t p k , k + 1 t 0 0 . . . 0 1 = P t p k + 1 t 0 1
(b) According to (a), whenever a defaultable bond enters state i at time t, it chooses credit class j, and moves to it with p i j t = Q S t + 1 = j X t = i . Now, before performing the actual transition from credit class i to credit class j, it “holds” for a time ϖ i j t in state i. All ϖ i j t are positive, integer-valued random variables with a probability mass function h i j t , m , where
h i j t , m = Q ϖ i j t = m S t + 1 = j , X t = i , t , m = 0 , 1 , 2 ,
The matrix H ˜ t , m = h i j t , m i , j K is of the form
H ˜ t , m = h 11 t , m h 12 t , m . . . h 1 k t , m h 1 , k + 1 t , m h 21 t , m h 22 t , m . . . h 2 k t , m h 2 , k + 1 t , m . . . . . . . . . . . . . . . h k 1 t , m h k 2 t , m . . . h k k t , m h k , k + 1 t , m 0 0 . . . . . . 1 m = 1 =
= H t , m h k + 1 t , m 0 1 m = 1 .
This means that all elements of H ˜ t , m are assumed finite and H ˜ t , 0 = 0 .
(c) Let f : K × K × N * R be any function, then
E f X t , S t + 1 , ϖ X t , S t + 1 t G t + 1 =
E f X t , S t + 1 , ϖ X t , S t + 1 t σ X t V σ S t + 1 ,
where N * = 0 , 1 , 2 , .
We call S t + 1 t = 0 the inherent G -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 G -inhomogeneous semi-Markov property remains, the relationships between the real-world probabilities p i j t = Q S t + 1 = j X t = i and probability mass functions
h i j t , m = Q ϖ i j t = m S t + 1 = j , X t = i ,
to the respective probabilities under Q T , 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 θ t 1 θ t + 1 be σ X t V σ S t + 1 -measurable for any t = 0 , 1 , , T 1 . Consequently, for all values of t:
θ t 1 θ t + 1 = g t X t , S t + 1
for any function g t : K × K R . In addition, we assume that θ t 1 θ t + ϖ X t , S t + 1 t is σ X t V σ S t + 1 V σ ϖ X t , S t + 1 t -measurable for any t = 0 , 1 , , T 1 . Equivalently,
θ t + 1 1 θ t + ϖ X t , S t + 1 t = f t + 1 X t , S t + 1 , ϖ X t , S t + 1 t
for some function f t + 1 : K × K × N * R . Assuming that X t , S t + 1 t = 0 follows a discrete-time G -inhomogeneous semi-Markov process under Q , then a discrete-time G -inhomogeneous semi-Markov process under Q T also follows, and, in addition,
a F p i j t : = Q T S t + 1 = j X t = i = g t i , j p i j t
b F h i j t , m = Q T ϖ i j t = m S t + 1 = j , X t = i = f t + 1 i , j , m h i j t , m .

4. The G Non-Homogeneous Markov Renewal Process and Change to Forward Martingale Measure

The above definition of the G -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 G -inhomogeneous Markov renewal process and study in detail the change to forward martingale measure of their founding parameters.
Let ( Ω , G , Q , G t ) be a complete filtered probability space, and Y n n = 0 a stochastic process taking values in the state space K = 1 , 2 , , k , k + 1 and, representing the credit state, the defaultable bond enters at the n-th jump. Let F Y the natural filtration of Y n n = 0 and assume that it is a subfiltration of G . Let T n be the family of random variables taking values in R + = 0 , such that 0 T 0 T 1 . . . expressing the time of the n-th jump ( T n ) . Let F T be the natural filtration of T n n = 0 and allow that F T F Y G .
The stochastic process Y n , T n is said to be a non-homogeneous Markov renewal process provided that, for every bounded or non-negative function f : K R , we have
E Q f Y n , Y n + 1 , T n + 1 T n , T n G n =
E Q f Y n , Y n + 1 , T n + 1 T n , T n σ Y n σ T n ,
for n = 0 , 1 , 2 , . Define by
Q Y n + 1 = j , T n + 1 T n x Y n = i , T n = t = Q i , j , x , t .
The family of probabilities Q i , j , x , t for i , j K , x R + , T N is called a non-homogeneous Markov kernel.
The function x Q i , j , x , t has all the properties of a distribution function except that
p i j t = lim x Q i , j , x , t = Q i , j , , t .
Indeed, one could easily find that
p i j t 0 , j K p i j t = 1 .
If p i j ( t ) = 0 for some pair i , j , then Q i , j , x , t = 0 for all x; we then define Q i , j , x , t / p i j t = 1 . With this convention, we define
G i , j , x , t = Q i , j , x , t p i j t for i , j K , x R + , T N .
From 17 and 20 , we get that
G i , j , x , t = Q T n + 1 T n x Y n = i , Y n + 1 = j , T n = t .
Remark 1.
Note that, in words, the increments T 1 T 0 and T 2 T 1 in the inhomogeneous case, unlike the homogeneous, are not conditionally independent, given the Markov chain Y 0 , Y 1 , Y 2 and the T 0 , T 1 jumping times.
We will need the following Lemma from Musiela and Rutkowski [31].
Lemma 1.
(Abstract version of Bayes Formula). Let G be a sub σ-field of the σ-field G . Now, consider ψ a random variable integrable in relation to Q ; then, Q T is a probability measure equivalent to Q with
d Q T d Q G T = η T , Q - a . s . ,
with the random variable η T strictly positive Q -a.s.; moreover, η T is Q integrable with E Q η T = 1 . Then,
E Q T ψ G t = E Q ψ η G t E Q ψ G t .
With the use of the above Lemma, we prove the following basic theorem:
Theorem 2.
Let ( Ω , G , Q , G t ) be a complete filtered probability space and Y n , T n a non- homogeneous Markov renewal process with Markov kernel Q i , j , x , t given by 4.2 . In addition, let the random variable θ T n 1 θ T n + 1 be σ Y n σ Y n + 1 σ T n σ ( T n + 1 T n ) measurable for all T n = 0 , 1 , , T 1 . Then, for every T n ,
θ T n 1 θ T n + 1 = F T n ( Y n , Y n + 1 , T n + 1 T n ) .
The pair Y n , T n n = 0 follows a G -non-homogeneous Markov renewal process under Q T and moreover
F Q i , j , x , t = F t i , j , x Q i , j , x , t ,
where F Q i , j , x , t is the transition probability matrix of the forward probabilities.
Proof. 
Using Lemma 1, we fix T n and, for any j K , we get
Q Y n + 1 = j , T n + 1 T n x G T n = E Q 1 Y n + 1 = j , T n + 1 T n x G T n
= E Q θ T 1 Y n + 1 = j , T n + 1 T n x G T n E Q θ T G T n
= θ T is a martingale under Q and G T n - measurable
= E Q θ T n 1 θ T 1 Y n + 1 = j , T n + 1 T n x G T n
= Tower property
= E Q E Q θ T n 1 θ T 1 Y n + 1 = j , T n + 1 T n x G T n + 1 G T n
= E Q θ T n 1 E Q θ T G T n + 1 1 Y n + 1 = j , T n + 1 T n x G T n
= E Q θ T n 1 θ T n + 1 1 Y n + 1 = j , T n + 1 T n x G T n
= E Q F T n ( Y n , Y n + 1 , T n + 1 T n ) 1 Y n + 1 = j , T n + 1 T n x G T n
= E Q f ( Y n , Y n + 1 , T n + 1 T n , T n ) G T n
= ( by 4.1 )
= E Q f ( Y n , Y n + 1 , T n + 1 T n , T n ) σ Y n σ T n .
Now, 22 actually shows that Q Y n + 1 = j , T n + 1 T n x G T n is a σ Y n σ T n -measurable random variable. Since σ Y n σ T n G T n , we get that
Q Y n + 1 = j , T n + 1 T n x G T n =
= Q Y n + 1 = j , T n + 1 T n x σ Y n σ T n .
Moreover, we have
F Q i , j , x , t = Q Y n + 1 = j , T n + 1 T n x Y n = i , T n = t
= E Q F T n ( Y n , Y n + 1 , T n + 1 T n ) 1 Y n + 1 = j , T n + 1 T n x Y n = i , T n = t
= F t ( i , j , x ) Q i , j , x , t .

5. Martingale Characterization of the G -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 G -inhomogeneous semi-Markov process X t , S t + 1 t = 0 . 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:
H t i = 1 X t = i , Υ t + 1 j = 1 S t + 1 = j and Φ t + 1 i j = 1 X t = i , S t + 1 = j .
and the processes
C 0 , t i j : the number of times the G - inhomogeneous semi-Markov process X t , S t + 1 t = 0 entered state i and selected to move to state j in the time interval 0 , t .
Z 0 , t i j m : the number of transitions of the G -inhomogeneous semi-Markov process X t , S t + 1 t = 0 from i to j that occurred in the time interval 0 , t remaining time m in state i.
Now, we arrive at the following relations for the above defined processes:
C 0 , t i j = u = 0 t 1 H u i Υ u + 1 j , Z 0 , t i j m = u = 0 t m Φ u + 1 i j H u + m j = u = 0 t m H u i Υ u + 1 j H u + m j .
At this point, it is important to note that, for a discrete time G non-homogeneous semi-Markov process X t , S t + 1 t = 0 , the sequences of transition probabilities P ˜ t and H ˜ t , m t , m = 0 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 M n , m be the set of all n × m matrices with elements from R . In addition, let SM n be the set of n × n stochastic matrices. We now provide the following definition:
Definition 1.
Let the sequence H t , m t , m M n , n ; then, we say that it is an H -sequence if the following conditions are satisfied:
H t , m 0 and m = 0 H t , m = U , f o r e v e r y t , m N *
where U M n , n are a matrix of ones.
Theorem 3.
Let there be X t , S t + 1 t = 0 G - non-homogeneous semi-Markov process, with sequences of probabilities P ˜ t t = 0 and H ˜ t , m t , m = 0 . For any two states i , j K , the processes
C 0 , t i j = C 0 , t i j u = 0 t 1 p i j u H u i a n d Y 0 , t i j m = Z 0 , t i j m u = 0 t m h i j u , m Φ u + 1 i j ,
for t N * , following G -martingales under Q . Conversely, let
P t t = 0 SM k + 1 a n d l e t H t , m t , m = 0 a n H - s e q u e n c e ,
then, the processes
C 0 , t i j = C 0 , t i j u = 0 t 1 p i j u H u i a n d Y 0 , t i j m = Z 0 , t i j m u = 0 t m h i j u , m Φ u + 1 i j ,
are G -martingales under Q , then the sequences of matrices
P ˜ t t = 0 a n d H ˜ t , m t , m = 0 ,
are equal, respectively, with the sequences of matrices
P t t = 0 a n d H t , m t , m = 0 .
In conclusion, the sequences of transition probabilities
P ˜ t t = 0 a n d H ˜ t , m t , m = 0 ,
are the unique sequences for which the processes
C 0 , t i j t = 0 a n d Y 0 , t i j m t , m = 0 ,
are G -martingales under Q .
Proof. 
It is sufficient to show:
E Q C 0 , t + 1 i j G t = C 0 , t i j and E Q Y 0 , t + 1 i j m G t = Y 0 , t i j m .
We assume that i j and, for i = j , the proof will be similar. For the first part of Equation 15 , we have
E Q C 0 , t + 1 i j C 0 , t i j G t =
= E Q u = 0 t H u i Υ u + 1 j u = 0 t p i j u H u i u = 0 t 1 H u i Υ u + 1 j + u = 0 t 1 p i j u H u i G t
= E Q H t i Υ t + 1 j p i j t H t i G t
= E Q H t i Υ t + 1 j p i j t σ X t = 0 .
Now, from the second part of Equation 28 , we have
E Q Y 0 , t + 1 i j m Y 0 , t i j m G t =
= E Q Z 0 , t + 1 i j m u = 0 t m + 1 h i j u , m Φ u + 1 i j Z 0 , t i j m + u = 0 t m h i j u , m Φ u + 1 i j G t
= E Q u = 0 t m + 1 Φ u + 1 i j H u + m j u = 0 t m + 1 h i j u , m Φ u + 1 i j u = 0 t m Φ u + 1 i j H u + m j + u = 0 t m h i j u , m Φ u + 1 i j G t
= E Q Φ t m + 2 i j H t + 1 j h i j t m + 1 , m Φ t m + 2 i j G t
= E Q H t m + 1 i Υ t m + 2 j H t + 1 j h i j t m + 1 , m G t
= E Q H t m + 1 i Υ t m + 2 j H t + 1 j h i j t m + 1 , m σ X t V σ S t + 1 = 0 .
Conversely, since C 0 , t i j is a martingale under Q , we have that
E Q C 0 , t + 1 i j C 0 , t i j G t = 0 .
Consequently,
E Q u = 0 t H u i Υ u + 1 j u = 0 t p i j u H u i u = 0 t 1 H u i Υ u + 1 j + u = 0 t 1 p i j u H u i G t
= E Q H t i Υ t + 1 j p i j t H t i G t = 0 .
From Equations 29 and 31 , we get
E Q H t i Υ t + 1 j p i j t H t i H t i Υ t + 1 j + p i j t H t i G t = 0 ,
or
E Q p i j t H t i p i j t H t i G t = E Q p i j t H t i p i j t H t i σ X t = 0 ,
from which we conclude that
P t = P ˜ t for every t N * .
Since Y 0 , t i j m is a martingale under Q , we have that
E Q Y 0 , t + 1 i j m Y 0 , t i j m G t = 0 .
Consequently,
E Q Z 0 , t + 1 i j m u = 0 t m + 1 h i j u , m Φ u + 1 i j Z 0 , t i j m + u = 0 t m h i j u , m Φ u + 1 i j G t
= E Q Φ t m + 2 i j H t + 1 j h i j t m + 1 , m Φ t m + 2 i j G t
= E Q H t m + 1 i Υ t m + 2 j H t + 1 j h i j t m + 1 , m G t = 0
From Equations 30 and 32 , we get
E Q Φ t m + 2 i j H t + 1 j h i j t m + 1 , m Φ t m + 2 i j Φ t m + 2 i j H t + 1 j + h i j t m + 1 , m Φ t m + 2 i j G t = 0
or
E Q Φ t m + 2 i j h i j t m + 1 , m h i j t m + 1 , m G t
= E Q H t m + 1 i Υ t m + 2 j h i j t m + 1 , m h i j t m + 1 , m G t
= E Q H t m + 1 i Υ t m + 2 j h i j t m + 1 , m h i j t m + 1 , m σ X t 1 V σ S t = 0 ,
from which we conclude that
H ˜ t , m = H t , m for every t , m N * .

6. Martingale Characterization of the G -Non-Homogeneous Markov Renewal Process

We start by providing in the form of a theorem a martingale characterization, for the G -non-homogeneous Markov renewal process Y n , T n n = 0 . Define the process
R 0 , T n i j x : the number of transitions of the G -IMRP from i to j that occurred in the interval 0 , T n with sojourn time to state i being less than x.
In addition, define the jump processes
Γ T n j = 1 Y T n = j , D r a = 1 T r + 1 T r = a .
Now, we have that
R 0 , T n i j x = a = 1 x r = 0 n 1 Γ T r i Γ T r + 1 j D r a .
Then, we define the following stochastic process:
R 0 , T n i j x = R 0 , T n i j x r = 0 n 1 Q i , j , x , t Γ T r i Γ T r + 1 j .
We will now prove the following martingale characterization theorem for a G -inhomogeneous Markov renewal theorem.
Theorem 4.
Let ( Ω , G , Q , G t ) be a complete filtered probability space and Y n , T n a non-homogeneous Markov renewal process with Markov kernel Q i , j , x , t given by 17 . For any two states i , j K , the process
R 0 , T n i j x = R 0 , T n i j x r = 0 n 1 Q i , j , x , T r Γ T r i Γ T r + 1 j ,
for T n N * follows G -martingales under Q . Conversely, let the inhomogeneous Markov kernel Q i , j , x , t ; then, if
R 0 , T n i j x = R 0 , T n i j x r = 0 n 1 Q i , j , x , T r Γ T r i Γ T r + 1 j ,
is a G -martingale under Q , then the inhomogeneous Markov kernel Q i , j , x , t is equal to the inhomogeneous Markov kernel Q i , j , x , t . In conclusion, the inhomogeneous Markov kernel Q i , j , x , t is the unique Markov kernel for which the processes R 0 , T n i j x are a G -martingale under Q .
Proof. 
It is sufficient to show:
E Q R 0 , T n + 1 i j x G T n = R 0 , T n i j x
We assume that i j and for i = j the proof will be similar. We get that
E Q R 0 , T n + 1 i j x R 0 , T n i j x G T n =
E Q ( a = 1 x r = 0 n Γ T r i Γ T r + 1 j D r a r = 0 n Q i , j , x , T r Γ T r i Γ T r + 1 j
a = 1 x r = 0 n 1 Γ T r i Γ T r + 1 j D r a + r = 0 n 1 Q i , j , x , T r Γ T r i Γ T r + 1 j G T n )
= E Q a = 1 x Γ T n i Γ T n + 1 j D r a Q i , j , x , T n Γ T n i Γ T n + 1 j G T n
= E Q Γ T n i Γ T n + 1 j a = 1 x D r a Q i , j , x , T n G T n
= E Q Γ T n i Γ T n + 1 j a = 1 x D r a Q i , j , x , T n σ Y n σ T n = 0 .
Conversely, let Q i , j , x , T n M k + 1 ; then, we will prove that, if the process
R 0 , T n i j x = R 0 , T n i j x r = 0 n 1 Q i , j , x , T r Γ T r i Γ T r + 1 j ,
is a G -martingale under Q , then the transition probabilities Q i , j , x , T n are equal to Q i , j , x , T n .
Since the process R 0 , T n i j x is a G -martingale under Q , we have that
E Q R 0 , T n + 1 i j x R 0 , T n i j x G T n = 0 .
Consequently,
E Q ( a = 1 x r = 0 n Γ T r i Γ T r + 1 j D r a r = 0 n Q i , j , x , T r Γ T r i Γ T r + 1 j
a = 1 x r = 0 n 1 Γ T r i Γ T r + 1 j D r a + r = 0 n 1 Q i , j , x , T r Γ T r i Γ T r + 1 j G T n )
= E Q Γ T n i Γ T n + 1 j a = 1 x D r a Q i , j , x , T n G T n
= E Q Γ T n i Γ T n + 1 j a = 1 x D r a Q i , j , x , T n σ Y n σ T n = 0 .
From 35 and 36 , we get that
E Q Γ T n i Γ T n + 1 j Q i , j , x , T n Q i , j , x , T n σ Y n σ T n = 0 ,
from which we get that
Q i , j , x , T n = Q i , j , x , T n .

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:
( a ) The random variable θ t 1 θ t + 1 be σ X t V σ S t + 1 -measurable for any t = 0 , 1 , , T 1 . As a result, we get that, for every such t, we have
θ t 1 θ t + 1 = g t X t , S t + 1
for some function g t : K × K R .
( b ) The random variable θ t + 1 1 θ t + ϖ X t , S t + 1 is σ X t V σ S t + 1 V σ ϖ X t , S t + 1 -measurable for all t = 0 , 1 , , T 1 . Equivalently,
θ t + 1 1 θ t + ϖ X t , S t + 1 = f t + 1 X t , S t + 1 , ϖ X t , S t + 1
for some function f t + 1 : K × K × N * R .
Let us now assume that the following Condition 1 holds:
Condition 1. ( a ) The random variable θ t 1 θ t + 1 is σ X t V σ S t + 1 -measurable for any t = 0 , 1 , , T 1 , and there exists a finite set, denoted by A , such that the random variable θ t 1 θ t + 1 admits the following representation:
θ t 1 θ t + 1 = 1 + a A g ^ t a X t g ˜ t + 1 a X t , S t + 1 .
for some functions g ^ t a : K R , and g ˜ t + 1 a : K × K R . ( b ) The random variable θ t 1 θ t + ϖ X t , S t + 1 is σ X t V σ S t + 1 V σ ϖ X t , S t + 1 -measurable for every t = 0 , 1 , , T 1 and there exists a finite set, denoted B , such that the random variable θ t 1 θ t + ϖ X t , S t + 1 t admits the following representation:
θ t + 1 1 θ t + ϖ X t , S t + 1 = 1 + b B f ^ t + 1 X t , S t + 1 f ˜ t + 1 X t , S t + 1 , ϖ X t , S t + 1
for any f ^ t + 1 : K × K R and f ˜ t + 1 : K × K × N * .
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 G -homogeneous Markov chain is a discrete -time G inhomogeneous semi-Markov process for which
p i j t = p i j f o r a l l i , j K a n d f o r a l l t = 0 , 1 , 2 , ,
h i j t , 1 = h i j 1 = 1 f o r a l l t = 0 , 1 , 2 , ,
h i j t , m = 0 f o r m > 1 a n d f o r a l l t = 0 , 1 , 2 , .
It is an immediate consequence then to check that, in the case of the G -homogeneous Markov chain model, Condition 1 coincides with Condition B . 2 in Bielecki and Rutkowski [28].
Theorem 5.
Let X t , S t + 1 t = 0 be the G -inhomogeneous semi-Markov process under Q . Suppose that Condition 1 holds. Then, we may choose as (a) of Condition 1
θ t 1 θ t + 1 = 1 + ν = 1 k + 1 l = 1 k + 1 g ^ t ν l X t Δ C 0 , t ν l ,
then any chosen function g ^ t i j i should satisfy
1 p i j t + l j k + 1 g ^ t i l i p i l t 1 p i j t g ^ t i j i > l j k + 1 g ^ t i l i p i l t 1 1 p i j t .
In addition, we may choose as (b) of Condition 1
θ ^ t + 1 1 θ ^ t + m = 1 + ϖ = 1 m ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l X t , S t + 1 Δ Y 0 , t + ϖ 1 ν l ϖ ,
then any chosen function f ^ t + 1 i j i , j should satisfy
1 h i j t , m h i j t , m 1 ϖ = 1 m h i j t , ϖ f ^ t + 1 i j i , j > 1 ϖ = 1 m h i j t , ϖ 1 .
Then, X t , S t + 1 t = 0 is a discrete time G -inhomogeneous semi-Markov process under the forward probability measure Q T , and, for every t = 0 , 1 , 2 , , T 1 , we have the following:
( a ) θ t = 1 + u = 1 t θ u 1 ν = 1 k + 1 l = 1 k + 1 g ^ u ν l X u g ˜ u + 1 ν l X u , S u + 1 .
( b ) F p i j t = p i j t 1 + g ^ t i j i l = 1 k + 1 p i l t g ^ t i l i .
( c ) F h i j t , r = h i j t , r 1 + f ^ t + 1 i j i , j m = 1 r f ^ t + 1 i j i , j h i j t m + 1 , m .
( d ) θ ^ t = 1 + u = 1 t 1 θ ^ u ϖ = 1 2 ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l i , j Δ Y 0 , t + ϖ 2 ν l ϖ .
In addition, the chosen functions g ^ t i j i and f ^ t + 1 i j i , j should be such that θ ^ t = θ t + 1 .
Proof. 
In order for Condition 1 be acceptable, it should satisfy two prerequisites. First, θ t should be a martingale process adapted to G t and secondly θ t should be strictly positive. From Equation 25 and Theorem 5, we have
Δ C 0 , t ν l = C 0 , t + 1 ν l C 0 , t ν l = u = 0 t H u ν Υ u + 1 l u = 0 t p ν l u H u ν u = 0 t 1 H u ν Υ u + 1 l + u = 0 t 1 p ν l u H u ν
= H t ν Υ t + 1 l p ν l t H t ν = 1 X t = ν , S t + 1 = l p ν l t 1 X t = ν
and so the process Δ C 0 , t ν l is σ X t V σ S t + 1 -measurable, i.e., it is a function; let us say g ˜ t + 1 ν l X t , S t + 1 , where
g ˜ t + 1 ν l i , j = 1 X t = ν i 1 S t + 1 = l j p ν l t
Now, we have
E Q ν = 1 k + 1 l = 1 k + 1 g ^ t ν l X t Δ C 0 , t ν l G t = . ν = 1 k + 1 l = 1 k + 1 g ^ t ν l X t E Q Δ C 0 , t ν l G t
= ν = 1 k + 1 l = 1 k + 1 g ^ t ν l X t E Q H t ν Υ t + 1 l p ν l t H t ν G t
= ν = 1 k + 1 l = 1 k + 1 g ^ t ν l X t E Q H t ν Υ t + 1 l p ν l t H t ν X t = ν = 0
From Equation 38 , it is apparent that we may choose as ( a ) of Condition 4.1 the following:
θ t 1 θ t + 1 = 1 + ν = 1 k + 1 l = 1 k + 1 g ^ t ν l X t Δ C 0 , t ν l .
Then, from 41 , we get
E Q θ t 1 θ t + 1 G t = due to 40 = 1 ,
and consequently
E Q θ t + 1 G t = θ t .
Thus, the stochastic process θ t as defined by 41 is a martingale under Q . In addition, θ t as defined by 41 should be strictly positive. It is easy to check that it is sufficient to have that
ν = 1 k + 1 l = 1 k + 1 g ^ t ν l X t Δ C 0 , t ν l > 1 .
Now, from 43 for any i , j K , we equivalently have
ν = 1 k + 1 l = 1 k + 1 g ^ t ν l i 1 X t = ν i 1 S t + 1 = l p ν l t > 1 ,
or, equivalently
l = 1 k + 1 g ^ t i l i 1 S t + 1 = l p ν l t > 1 ,
from which we get
g ^ t i j i l = 1 k + 1 g ^ t i l i p i l t > 1 ,
and finally
g ^ t i j i > l j k + 1 g ^ t i l i p i l t 1 1 p i j t .
Thus, any chosen function g ^ t i j i should satisfy 44 .
From Equation 41 , since θ 0 = 1 , the process θ t is given by
θ t = 1 + u = 1 t θ u 1 ν = 1 k + 1 l = 1 k + 1 g ^ u 1 ν l X u 1 Δ C 0 , u 1 ν l .
Since Condition 1 holds, then the conditions of Theorem 5 are valid and thus we have
F p i j t = p i j t 1 + ν = 1 k + 1 l = 1 k + 1 g ^ t ν l i 1 X t = ν i 1 S t + 1 = l j p ν l t
= p i j t 1 + l = 1 k + 1 g ^ t i l i 1 S t + 1 = l j p i l t
= p i j t 1 + g ^ t i j i l = 1 k + 1 p i l t g ^ t i l i .
From Equation 46 , it is apparent that the choice of any function g ^ t i j i , apart from satisfying 44 , which guarantees the positivity of F p i j t , should be such that
p i j t 1 + g ^ t i j i l = 1 k + 1 p i l t g ^ t i l i 1 ,
from which we get that
g ^ t i j i 1 p i j t + l j k + 1 g ^ t i l i p i l t 1 p i j t .
Now, from Theorem 5, we have
Δ Y 0 , t + m 1 ν l m = Y 0 , t + m ν l m Y 0 , t + m 1 ν l m
= Y 0 , t + m ν l m u = 0 t h ν l u , m Φ u + 1 ν l Y 0 , t + m 1 ν l m + u = 0 t 1 h ν l u , m Φ u + 1 ν l
= u = 0 t Φ u + 1 ν l H u + m l u = 0 t h ν l u , m Φ u + 1 ν l u = 0 t 1 Φ u + 1 ν l H u + m l + u = 0 t 1 h ν l u , m Φ u + 1 ν l
= H t ν Υ t + 1 l H t + m l h ν l t , m
= 1 X t = ν , S t + 1 = l 1 ϖ ν l t = m h ν l t , m .
and so the process Δ Y 0 , t + m 1 ν l m is a function of the form
f ˜ t + 1 X t , S t + 1 , ϖ X t , S t + 1 t
where
f ˜ t + 1 i , j , r = 1 X t = ν , S t + 1 = l i , j 1 ϖ ν l t = m r h ν l t , r .
It is apparent, from Equations 48 and 49 , that we may choose part b) of Condition 1 to be of the following form:
θ ^ t + 1 1 θ ^ t + m = 1 + ϖ = 1 m ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l X t , S t + 1 Δ Y 0 , t + ϖ 1 ν l ϖ .
Now, we have
E Q θ ^ t + 1 1 θ ^ t + m G t + 1 = 1 + ϖ = 1 m ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l X t , S t + 1 E Q Δ Y 0 , t + ϖ 1 ν l ϖ G t + 1
= 1 + ϖ = 1 m ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l X t , S t + 1 E Q H t ν Υ t + 1 l H t + ϖ l h ν l t , ϖ G t + 1
= 1 + ϖ = 1 m ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l X t , S t + 1
× E Q H t ν Υ t + 1 l H t + ϖ l h ν l t , ϖ X t = ν , S t + 1 = l = 1
Thus, we get that
E Q θ ^ t + m G t + 1 = θ ^ t + 1
and consequently the process θ ^ t is a martingale under Q . In addition, θ ^ t as defined by 50 should be strictly positive. It is easy to check that it is sufficient to have
ϖ = 1 m ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l X t , S t + 1 Δ Y 0 , t + ϖ 1 ν l ϖ > 1 .
For any i, j, r, it is equivalent to have
1 + ϖ = 1 m ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l i , j 1 X t = ν , S t + 1 = l i , j 1 ϖ ν l t = ϖ r h ν l t , ϖ
= 1 + ϖ = 1 m f ^ t + 1 i j i , j 1 ϖ i j t = ϖ r h i j t , ϖ
= 1 + f ^ t + 1 i j i , j ϖ = 1 m f ^ t + 1 i j i , j h i j t , ϖ > 0 .
Thus, the choice of f ^ t + 1 i j i , j should be such that
f ^ t + 1 i j i , j > 1 ϖ = 1 m h i j t , ϖ 1
Since Condition 1 holds, then the conditions of Theorem 5 are valid and thus we have
F h i j t , m = h i j t , m 1 + ϖ = 1 m ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l i , j Δ Y 0 , t + ϖ 1 ν l ϖ
= h i j t , m ( 1 + ϖ = 1 m ν = 1 k + 1 l = 1 k + 1 f ^ t + 1 ν l i , j 1 X t = ν , S t + 1 = l i , j
× 1 ϖ ν l t = m r h ν l t , ϖ )
= h i j t , m 1 + ϖ = 1 m f ^ t + 1 i j i , j 1 ϖ i j t = m r h i j t , ϖ
= h i j t , m 1 + f ^ t + 1 i j i , j ϖ = 1 m f ^ t + 1 i j i , j h i j t , ϖ .
From Equation 55 , it is apparent that any choice of a function f ^ t + 1 i j i , j apart from satisfying Equation 23 should also satisfy the following relation:
h i j t , m 1 + f ^ t + 1 i j i , j ϖ = 1 m f ^ t + 1 i j i , j h i j t , ϖ 1 ,
which is equivalent with
f ^ t + 1 i j i , j 1 h i j t , m h i j t , m 1 ϖ = 1 m h i j t , ϖ .
From Equation 50 , since the initial information is included in the interval [ 0 , 1 ) , we may assume that θ ^ 1 = 1 and then we get
θ ^ t = 1 + u = 1 t 1 θ ^ u ϖ = 1 2 ν = 1 k + 1 l = 1 k + 1 f ^ u ν l i , j Δ Y 0 , u + ϖ 2 ν l ϖ .
It is apparent that the choices of g ^ t i j i and f ^ t + 1 i j i , j , in addition to the already created restrictions, should also be such that
θ t + 1 = θ ^ t for every t = 0 , 1 , , T .

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 T * is the present time. We assume that the available data span are for the time window 0 , T * . These data will consist of the following type:
(i) From historical data on credit migrations, we firstly estimate the matrix sequences P ˜ t t = 0 T * 1 and H ˜ t , m t , m = 0 T * 1 during our time window 0 , T * . 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 B t , T , for every t = 0 , 1 , , T * and T = t + 1 , , T * .
(iii) The observed values of D i δ t , T , for every i = 1 , 2 , , k , t = 1 , 2 , , T * , and T = t + 1 , , T * .
(iv) Observed values on recovery rates, δ , in the time window 0 , T * .
Our target is to estimate the sequences F P ˜ t t = 0 T * 1 and F H ˜ t , m t , m = 0 T * 1 , within 0 , T * , that is, the time window. In fact, we do not need to calculate F P ˜ t t = 0 T * 1 and F H ˜ t , m t , m = 0 T * 1 separately as it is shown below. What is needed is to find the entrance probabilities E ˜ s , t defined below, which as it is proved in what follows are functions of the Hadamard products F P t F H t , m .

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 t 0 , T * 1 and i , j = 1 , 2 , , k , assume:
g t i j X t = g t i X t .
(ii) For any t 0 , T * 1 , i = 1 , 2 , , k , and j = 1 , 2 , , k + 1 , assume:
f t + 1 i j X t , S t + 1 = f t + 1 i X t .
The consequences of the Assumption 8.1 on the results of Theorem 7 are the following:
(a) For all i = 1 , 2 , , k , j = 1 , 2 , , k + 1 and for every t 0 , T * 1
F p i j t = p i j t 1 + g t i i l i k + 1 p i l t g t i i = p i j t 1 + p i i t g t i i
= g ˜ t i p i j t .
(b) For all i , j = 1 , 2 , , k and j = k + 1 and for every t 0 , T * 1 ,
F h i j t , r = h i j t , r 1 + f t + 1 i i 1 m = 1 r h i j t m + 1 , m
= f ˜ t i , r h i j t , r .

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:
Q the bond enters credit class j at time t it entered credit class i at time s
= Q X t = j X s = i = e i j s , t .
and E s , t = e i j ( s , t ) , i , j = 1 , 2 , , k . In addition, let E ˜ s , t be the matrices of entrance probabilities such that
E ˜ s , t = E s , t e k + 1 s , t 0 1 t s = 1 ,
where e k + 1 s , t = e 1 k + 1 s , t , e 2 k + 1 s , t , , e k k + 1 s , t , with
e i k + 1 s , t = Q X t = k + 1 X s = i for i = 1 , 2 , , k .
From Vasileiou and Vassiliou [3] p.182, we get
E s , t = m = 0 t s P s H s , m E s + m , t for s < t ,
with E s , s = I for s = 0 , 1 , , and
e k + 1 s , t = m = 0 t s P s H s , m e k + 1 s + m , t + p k + 1 s h k + 1 s , t s ,
with e k + 1 s , s = 0 for s = 0 , 1 , T * 1 .
From Theorem 1, we know that the change to forward measure preserves the semi-Markov structure, hence
Q T τ > T X t = i = 1 Q T τ T X t = i = 1 s = t T F e i , k + 1 s , T .
Using 9 and 67 , we get
F e k + 1 t , t + 1 = B t , t + 1 D i δ t , t + 1 1 δ B t , t + 1 for every t 0 , T * 1 ,
and more generally we have
B t , t + n D i δ t , t + n 1 δ B t , t + n = k = 0 n 1 F e k + 1 t , t + 1 for every t 0 , T * 1 .

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 t = 0 , 1 , , T * 1 , and i = 1 , 2 , , k , using 68 , we find
F e i , k + 1 t , t + 1 = B t , t + 1 D i δ t , t + 1 1 δ B t , t + 1 .
Since the forward measure preserves the semi-Markov structure using relation 66 , we get
F e k + 1 t , t + 1 = F p k + 1 t F h t , 1 ,
and, therefore, for all t = 0 , 1 , , T * 1 , and i = 1 , 2 , , k using 70 and 71 , we obtain
g ˜ t i f ˜ t i , 1 = B t , t + 1 D i δ t , t + 1 1 δ B t , t + 1 p i , k + 1 t h i , k + 1 t , 1 .
Now, from Theorem 5, we know that the forward measure preserves the semi-Markov process, hence, from 65 , we get that
F E t , t + 1 = F P t F H t , 1 ;
therefore, for all t = 0 , 1 , , T * 1 , we find
F E t , t + 1 = g ˜ t 1 f ˜ t 1 , 1 p 11 t h 11 t , 1 . . . g ˜ t 1 f ˜ t 1 , 1 p 1 k t h 1 k t , 1 . . . . . . . . . g ˜ t k f ˜ t k , 1 p k 1 t h k 1 t , 1 . . . g ˜ t k f ˜ t k , 1 p k 1 t h k 1 t , 1 ,
and finally
F E ˜ t , t + 1 = F E t , t + 1 F e k + 1 t , t + 1 0 1 .
Step 2. Using 69 and since, in Step 1 by relation 70 , we have evaluated F e i , k + 1 t + 1 , t + 2 for all t and all i for t = 0 , 1 , . . , T * 2 , for i = 1 , 2 , , k , we obtain
F e i , k + 1 t , t + 2 = B t , t + 2 D i t , t + 2 1 δ B t , t + 2 F e i , k + 1 t + 1 , t + 2
From 66 , since the change to forward measure preserves the semi-Markov process, we find that
F e k + 1 t , t + 2 = F P t F H t , 1 F e k + 1 t + 1 , t + 2 + F p k + 1 t F h k + 1 t , 2
Define by B i j the i j -element of the matrix B , then 77 could be written as
F e i , k + 1 t , t + 2 = l = 1 k F P t F H t , 1 i l F e l , k + 1 t + 1 , t + 2 +
F p i , k + 1 t F h i , k + 1 t , 2 ,
and
F e i , k + 1 t , t + 2 = l = 1 k F P t F H t , 1 i l F e l , k + 1 t + 1 , t + 2
+ g ˜ t i f ˜ t i , 2 p i , k + 1 t h i , k + 1 t , 2
Solve the system 76 and 79 to obtain
g ˜ t i f ˜ t i , 2 for i = 1 , 2 , , k and t = 0 , 1 , , T * 1 .
Hence, now we are in a position to evaluate
F P t F H t , t + 2 =
g ˜ t 1 f ˜ t 1 , 2 p 11 t h 11 t , 2 . . . g ˜ t 1 f ˜ t 1 , 2 p 1 k t h 1 k t , 2 . . . . . . . . . g ˜ t k f ˜ t k , 2 p k 1 t h k 1 t , 2 . . . g ˜ t k f ˜ t k , 2 p k 1 t h k 1 t , 2
Since the change to forward measure preserves the semi-Markov property from 65 , we get
F E t , t + 2 = F P t F H t , t + 1 E t + 1 , t + 2 + F P t F H t , t + 2 .
Now, from 82 and using 81 and 75 , we evaluate for t = 0 , 1 , , T * 2 the entrance probabilities F E t , t + 2 . Hence, Step 2 is concluded by collecting or the evaluations done in the matrices:
F E ˜ t , t + 1 = F E t , t + 1 F e k + 1 t , t + 1 0 1 .
We repeat the methodology indicated by Step 2 for evaluating all desired values of F E ˜ t , t + n for n = 1 , 2 , , T * t .
If we wish to find separate values of g ˜ t i and f ˜ t i , m for our time window 0 , T * , do as follows. Let, for every specific t, M i j t to be the maximum of m with h i j t , m strictly positive, and define
h i j t , M i j t = 1 m = 1 M i j t 1 h i j t , m .
From Theorem 5, the semi-Markov structure is preserved therefore
m = 1 F h i j t , m = 1 for every t 0 , T * 1 and i , j = 1 , 2 , , k ,
from which we get
m = 1 M i j t f ˜ t i , m h i j t , m = 1 ,
and consequently
g ˜ t i = m = 1 M i j t g ˜ t i f ˜ t i , m h i j t , m ,
which provides the values g ˜ t i i , and, consequently, the corresponding values
f ˜ t i , 1 , f ˜ t i , 2 , , f ˜ t i , M i j t .
Now, from 60 and 61 , we get that
g t i i = g ˜ t i 1 p i i t ,
and
f t + 1 i j i = f ˜ t i , m 1 1 m = 1 r h i j t m + 1 , m .

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
Aaa , Aa , A , Baa , Ba , B , Caa , d e f a u l t .
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 8 × 8 matrices in the present application, which will take too much space. That is, the state space will be S = 1 , 2 , 3 , 4 , 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 0 , 3 . The observed values of B ( t , T ) are provided by the market and for our present time window are the following:
Default free bonds for the time window 0 , 3 .
B 0 , 0 = 1 .
B 0 , 1 = 0.9 , B 1 , 1 = 1 .
B 0 , 2 = 0.87 , B 1 , 2 = 0.888 , B 2 , 2 = 1 .
B 0 , 3 = 0.8 , B 1 , 3 = 0.84 , B 2 , 3 = 0.86 , B ( 3 , 3 ) = 1 .
The observed values of D i δ t , T for i = 1 , 2 , 3 and for our time window are the following:
Term structure of defaultable bonds for the time window 0 , 3 .
GRADE 1: D 1 δ 0 , 0 = 1 .
D 1 δ 0 , 1 = 0.894 , D 1 δ 1 , 1 = 1 .
D 1 δ 0 , 2 = 0.81 ,     D 1 δ 1 , 2 = 0.882 , D 1 δ 2 , 2 = 1 .
D 1 δ 0 , 3 = 0.73 ,     D 1 δ 1 , 3 = 0.78 ,   D 1 δ 2 , 3 = 0.845 , D 1 δ 3 , 3 = 1 .
GRADE 2: D 2 δ 0 , 0 = 1 .
D 2 δ 0 , 1 = 0.884 , D 2 δ 1 , 1 = 1 .
D 2 δ 0 , 2 = 0.78 ,     D 2 δ 1 , 2 = 0.875 , D 2 δ 2 , 2 = 1 .
D 2 δ 0 , 3 = 0.70 ,     D 2 δ 1 , 3 = 0.75 ,   D 2 δ 2 , 3 = 0.82 , D 2 δ 3 , 3 = 1 .
GRADE 3: D 3 δ 0 , 0 = 1 .
D 3 δ 0 , 1 = 0.875 , D 3 δ 1 , 1 = 1 .
D 3 δ 0 , 2 = 0.72 ,     D 3 δ 1 , 2 = 0.85 , D 3 δ 2 , 2 = 1 .
D 3 δ 0 , 3 = 0.65 ,     D 3 δ 1 , 3 = 0.71 ,   D 3 δ 2 , 3 = 0.78 , D 3 δ 3 , 3 = 1 .
In the above, the physical architecture of the defaultable bonds in the arbitrage free Market could be observed.
The value of δ is 0.3 . 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 P t t = 0 2 and are the following:
P 0 = 0.7 0.15 0.05 0.01 0.10 0.65 0.1 0.15 0.08 0.12 0.55 0.25 0 0 0 1 , P 1 = 0.72 0.12 0.06 0.1 0.08 0.67 0.11 0.14 0.04 0.10 0.53 0.33 0 0 0 1 ,
P 2 = 0.66 0.14 0.07 0.13 0.09 0.60 0.12 0.19 0.05 0.11 0.50 0.34 0 0 0 1 .
(b) The sequence of the real world duration probabilities is
For t = 0
H 0 , 1 = 0.25 0.2 0.23 0.15 0.18 0.21 0.24 0.17 0.16 0.14 0.25 0.21 0 0 0 1 , H 0 , 2 = 0.18 0.17 0.2 0.30 0.15 0.18 0.21 0.33 0.14 0.12 0.18 0.36 0 0 0 0 ,
H 0 , 3 = 0.128 0.144 0.125 0.35 0.147 0.096 0.144 0.38 0.128 0.063 0.125 0.041 0 0 0 0 .
For t = 1
H 1 , 1 = 0.24 0.18 0.22 0.16 0.19 0.20 0.21 0.24 0.18 0.16 0.23 0.28 0 0 0 1 , H 1 , 2 = 0.19 0.21 0.24 0.22 0.18 0.22 0.17 0.23 0.18 0.21 0.20 0.19 0 0 0 0 ,
For t = 2
H 2 , 1 = 0.2 0.12 0.21 0.23 0.17 0.18 0.20 0.27 0.16 0.15 0.20 0.35 0 0 0 1 .
We are interested for our purposes in finding the sequence of forward probabilities F E t , t + m for t = 0 , 1 , 2 and m = 1 , 2 , 3 . 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:
F E 0 , 1 = 0.12 0.02 0.01 0.01 0.02 0.13 0.02 0.03 0.01 0.01 0.10 0.40 0 0 0 1 , E 0 , 1 = 0.18 0.03 0.01 0.02 0.02 0.14 0.03 0.03 0.02 0.02 0.14 0.06 0 0 0 1
F E 1 , 2 = 0.08 0.01 0.01 0.01 0.01 0.09 0.02 0.02 0.01 0.02 0.08 0.06 0 0 0 1 , E 1 , 2 = 0.15 0.03 0.02 0.04 0.02 0.14 0.02 0.06 0.01 0.02 0.11 0.11 0 0 0 1
F E 2 , 3 = 0.11 0.02 0.01 0.03 0.02 0.14 0.03 0.07 0.01 0.02 0.12 0.13 0 0 0 1 , E 2 , 3 = 0.13 0.02 0.02 0.03 0.02 0.11 0.02 0.05 0.01 0.02 0.11 0.11 0 0 0 1
F E 0 , 2 = 0.37 0.08 0.03 0.09 0.04 0.33 0.06 0.13 0.01 0.02 0.13 0.11 0 0 0 1 , E 0 , 2 = 0.15 0.03 0.01 0.04 0.02 0.14 0.03 0.06 0.01 0.02 0.11 0.10 0 0 0 1
F E 1 , 3 = 0.14 0.04 0.02 0.03 0.06 0.35 0.04 0.09 0.01 0.03 0.13 0.09 0 0 0 1 , E 1 , 3 = 0.24 0.05 0.04 0.03 0.06 0.22 0.05 0.04 0.02 0.04 0.17 0.08 0 0 0 1
F E 0 , 3 = 0.13 0.05 0.01 0.05 0.02 0.11 0.03 0.06 0.01 0.02 0.06 0.08 0 0 0 1 , E 0 , 3 = 0.15 0.04 0.02 0.05 0.04 0.11 0.03 0.07 0.02 0.02 0.11 0.12 0 0 0 1

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No real data were used, only representative data.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Carty, L.V.; Fons, J.S. Measuring changes in corporate credit quality. J. Fixed Income 1994, 4, 27–41. [Google Scholar] [CrossRef]
  2. Duffie, D.; Saita, L.; Wang, K. Multiperiod corporate default prediction with stochastic covariates. J. Financ. Econ. 2007, 83, 635–665. [Google Scholar] [CrossRef] [Green Version]
  3. Vasileiou, A.; Vassiliou, P.-C.G. An inhomogeneous semi-Markov model for the term Stucture of credit risk spreads. Adv. Appl. Prob. 2006, 38, 171–198. [Google Scholar] [CrossRef]
  4. D’Amico, G.; Janssen, J.; Manca, R. Homogeneous semi-Markov reliability models for credit risk management. Des. Econ. Finan. 2005, 28, 79–93. [Google Scholar] [CrossRef]
  5. Howard, D.T. Dynamic Probabilistic Systems. Vol. I and II; John Wiley: Hoboken, NJ, USA, 1971. [Google Scholar]
  6. D’Amico, G.; Janssen, J.; Manca, R. Valuing credit default swap in a non-homogeneous semi-Markov rating based model. Comput. Econom. 2007, 29, 119–138. [Google Scholar] [CrossRef]
  7. D’Amico, G.; Janssen, J.; Manca, R. Discete time non-homogeneous semi Markov reliability transition credit risk models and the default distributions functions. Comput. Econom. 2011, 38, 465–481. [Google Scholar] [CrossRef]
  8. D’Amico, G. A semi-Markov maintainance model with credit rating application. IMA J. Managem. Math. 2009, 20, 51–58. [Google Scholar] [CrossRef]
  9. D’Amico, G.; Janssen, J.; Manca, R. Initial and final backward and forward discrete time non-homogeneous semi-Markov credit risk models. Methodol. Comput. Appl. Prob. 2010, 12, 215–225. [Google Scholar] [CrossRef]
  10. D’Amico, G.; Janssen, J.; Manca, R. Downward migration credit risk problem: A non-homogeneous backward semi-Markov reliability approach. J. Oper. Res. Soc. 2016, 67, 393–401. [Google Scholar] [CrossRef]
  11. Vassiliou, P.-C.G.; Vasileiou, A. Asymptotic behaviour of the survival probabilities in an inhomogeneous semi-Markov model for the migration process in credit risk. Linear Algebra Its Appl. 2013, 438, 2880–2903. [Google Scholar] [CrossRef]
  12. Vassiliou, P.-C.G. Fuzzy semi-Markov migration process in credit risk. Fuzzy Sets Syst. 2013, 223, 39–58. [Google Scholar] [CrossRef]
  13. Huang, J. Review on term structure. In Proceedings of the 6th International Conference on Electronic, Mechanical, Information and Management Society, Shenyang, China, 1–3 April 2016; Volume 40, pp. 1278–1281. [Google Scholar]
  14. D’Amico, G.; Manca, R.; Salvi, G. Bivariate semi-Markov reward chain and credit spreads. IMA J. Manag. Math. 2016, 27, 529–556. [Google Scholar] [CrossRef]
  15. D’Amico, G.; Lika, A.; Petroni, F. Risk management of pention fund: A model for salary evolution. Int. J. Financ. Stud. 2019, 7, 44. [Google Scholar] [CrossRef] [Green Version]
  16. D’Amico, G. Stochastic dividend discount model: Risc and return. Markov Process. Relat. Fields 2017, 23, 349–376. [Google Scholar]
  17. Magni, C.; Magnoli, S.; Marchioni, A. Rating firms and sensitivity analysis. J. Oper. Res. 2019, 71, 1940–1958. [Google Scholar] [CrossRef]
  18. Wu, L.; Zhuang, Y.; Li, W. A new default intensity model with Fuzziness and hesitation. Int. J. Comput. Intell. Syst. 2016, 9, 340–350. [Google Scholar] [CrossRef] [Green Version]
  19. Wu, L.; Mei, X.; Sun, J. A new default probability calculation formula and its application under uncertain enviroments. Discret. Dyn. Nat. Soc. 2018, 2018, 3481863. [Google Scholar] [CrossRef] [Green Version]
  20. Puneet, P.; Selvamuthu, D.; D’Amico, G.; Manca, R. Portfolio optimization of credit risky bonds: A semi-Markov approach. Financ. Innov. 2020, 6, 25. [Google Scholar]
  21. De Blassis, R. The price leadership share: A new measure of price discovery in financial markets. Ann. Financ. 2020, 16, 381–405. [Google Scholar] [CrossRef]
  22. D’Amico, G.; Manca, R.; Salvi, G. Bivariate semi-Markov process for counterparty credit risk. Commun. Stat. Theory Methods 2014, 43, 1503–1522. [Google Scholar] [CrossRef]
  23. Vassiliou, P.-C.G. Semi-Markov migration process in a stochastic market in credit risk. Linear Algebra Its Appl. 2014, 450, 13–43. [Google Scholar] [CrossRef]
  24. Lando, D.; Skodeberg, T.M. Analyzing rating transitions and rating drift with continuous observations. J. Bank. Financ. 2002, 26, 423–444. [Google Scholar] [CrossRef] [Green Version]
  25. Bingham, N.H.; Kiesel, R. Risk-Neutral Valuation, Pricing and Hedging Financial Derivatives; Springer: London, UK, 1998. [Google Scholar]
  26. Shreve, S.E. Stochastic Calculus for Finance, Vol II, Continuous-Time Models; Springer: New York, NY, USA, 2004. [Google Scholar]
  27. Vassiliou, P.-C.G. Discrete-Time Asset Pricing Models in Applied Stochastic Finance; John Wiley: Hoboken, NJ, USA, 2010. [Google Scholar]
  28. Bielecki, T.R.; Rutkowski, M. Credit Risk: Modeling, Valuation, and Hedging; Springer: Berlin, Germany, 2002. [Google Scholar]
  29. Vassiliou, P.-C.G.; Papadopoulou, A.A. Nonhomogeneous semi-Markov systems and maintainability of the state sizes. J. Appl. Prob. 1992, 29, 519–534. [Google Scholar] [CrossRef]
  30. Pyke, R. Markov Renewal Processes with Finitely Many States. Ann. Math. Stat. 1961, 32, 1243–1259. [Google Scholar] [CrossRef]
  31. Musiela, M.; Rutkowski, M. Martingale Methods in Financial Modelling; Springer Science: Berlin/Heidelberg, Germany, 1997. [Google Scholar]
  32. Lee, E.T. Statistical Methods for Survival Data Analysis, 2nd ed.; John Wiley: Hoboken, NJ, USA, 2003. [Google Scholar]
  33. McClean, S.I.; Gribbin, J.O. Estimation for incomplete manpower data. Appl. Stoch. Models Data Anal. 1987, 3, 13–25. [Google Scholar] [CrossRef]
  34. McClean, S.I.; Gribbin, J.O. A non-parametric competing risks model for manpower planning. Appl. Stoch. Models Data Anal. 1991, 7, 327–341. [Google Scholar] [CrossRef]
  35. McClean, S.I.; Montgomery, E.; Ugwuowo, F. Inhomogeneous continuous-time Markov and semi-Markov manpower models. Appl. Stoch. Models Data Anal. 1997, 13, 191–198. [Google Scholar] [CrossRef]
  36. McClean, S.I.; Montgomery, E. Estimation for semi-Markov manpower models in a stochastic enviroment. In Semi-Markov Models and Applications; Janssen, J., Limnios, N., Eds.; Kluwer: Alphen aan den Rijn, The Netherland, 1999; pp. 219–227. [Google Scholar]
  37. Crouchy, M.; Galai, D.; Mark, R. Prototype risk rating system. J. Bank. Financ. 2001, 25, 27–41. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vassiliou, P.-C.G. Non-Homogeneous Semi-Markov and Markov Renewal Processes and Change of Measure in Credit Risk. Mathematics 2021, 9, 55. https://doi.org/10.3390/math9010055

AMA Style

Vassiliou P-CG. Non-Homogeneous Semi-Markov and Markov Renewal Processes and Change of Measure in Credit Risk. Mathematics. 2021; 9(1):55. https://doi.org/10.3390/math9010055

Chicago/Turabian Style

Vassiliou, P.-C.G. 2021. "Non-Homogeneous Semi-Markov and Markov Renewal Processes and Change of Measure in Credit Risk" Mathematics 9, no. 1: 55. https://doi.org/10.3390/math9010055

APA Style

Vassiliou, P. -C. G. (2021). Non-Homogeneous Semi-Markov and Markov Renewal Processes and Change of Measure in Credit Risk. Mathematics, 9(1), 55. https://doi.org/10.3390/math9010055

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