Advances in Decision Sciences

Volume 2012 (2012), Article ID 123635, 12 pages

doi:10.1155/2012/123635

Research Article

Monounireducible Nonhomogeneous Continuous Time Semi-Markov Processes Applied to Rating Migration Models

Guglielmo D'Amico,¹ Jacques Janssen,² and Raimondo Manca³

¹Department of Pharmacy, University “G. d'Annunzio” of Chieti, Via dei Vestini 31, 66013 Chieti, Italy

²CESIAF, EURIA, University of Bretagne Occidentale, 6 Avenue L. Gorgeu, CS 93837, 29238 Brest Cedex 3, France

³MEMOTEF Department, University “La Sapienza” of Roma, Via del Castro Laurenziano 9, 00161 Roma, Italy

Received 3 April 2012; Accepted 11 September 2012

Academic Editor: C. D. Lai

Copyright © 2012 Guglielmo D'Amico et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Monounireducible nonhomogeneous semi- Markov processes are defined and investigated. The mono- unireducible topological structure is a sufficient condition that guarantees the absorption of the semi-Markov process in a state of the process. This situation is of fundamental importance in the modelling of credit rating migrations because permits the derivation of the distribution function of the time of default. An application in credit rating modelling is given in order to illustrate the results.

1. Introduction

Semi-Markov processes (SMPs) are a generalization of Markov processes in which the waiting time distributions before the occurrence of a transition are modelled by any kind of distribution function; see [15]. This means that, on the contrary of Markov processes, it is possible to use also no memoryless distributions which determine a duration effect.

The duration effect affirms that the time the system is in a state influences the system’s transition probability. One way to detect and quantify this effect, in a SMP, is by using backward and forward recurrence time processes associated to the SMP.

In [5, 10] general distributions of the transition probabilities of SMP with backward and forward times are investigated for discrete time nonhomogeneous and for continuous time homogeneous processes, respectively. In these papers a credit risk application is also described. In this paper a further generalization is presented. In fact, a duration-dependent semi-Markov model is presented for continuous time nonhomogeneous processes. This extension is motivated by theoretical reasons as well by the practical need of making an efficient rating migration model available. In the paper the mono-unreducible topological structure for nonhomogeneous continuous time semi-Markov processes (NHCTSMPs) is introduced and analysed. It represents a sufficient condition for the absorption of the process in the down states of the system. Then, it represents the ideal environment in which to treat credit risk models.

The paper is organized as follows. Section 2 gives definitions and notations on NHCTSMP and introduces backward and forward processes and their joint distributions together with that of the SMP. Section 3 studies the monounreducible topological structure and derives the distribution function of going into default state. Section 4 presents a credit risk model and a numerical example.

2. Continuous Time Non-Homogeneous Semi-Markov Process

Let us consider two sequences of random variables defined on a complete; filtered probability space as the following:

(i), representing the state at the th transition;(ii) representing the time of the th transition.

We suppose that is a Markov renewal process of kernel :

We also know that.

is the transition matrix of the embedded Markov chain .

Let us introduce the probability; the process will leave state within time :

It is possible to define the distribution function of the waiting time in each state , given that the state, successively, occupied is known

The related probabilities can be obtained by means of the following formula:

Denote by , the NHCTSMP is defined as . It represents for each waiting time the state occupied by the embedded Markov chain.

The transition probabilities of are defined in the following way:

They are obtained by solving the following system of integral equations:

where represents the Kronecker symbol.

The part of formula (2.7) gives the probability that the process does not have transitions up to time given that it starts in state at time .

The term considers the permanence of the system in state up to the time , where a transition in state occurs. After the transition, the system will move to state following one of all the possible trajectories going from state at time to state at time . All possible states and times are considered by the summation and the integration.

With the aim of considering a general duration-dependent model, which is of fundamental importance in many applicative domains such as credit rating modelling and reliability theory, we introduce the recurrence time processes. These processes were investigated by many authors; see, for example, [4–6].

Given , we define the following recurrence processes:

is called the non-homogeneous backward time (or age) process, and is the non-homogeneous forward time (or residual time) process.

The recurrence time processes complement to a Markov process with respect to , , then for any bounded -measurable function it results

The relation remains true if we interchange with ; see [7].

Figure 1 represents a trajectory of a SMP and the backward and forward processes. In Figure 1 we have and , a starting backward , a starting forward , a final backward , and a final forward .

Figure 1: Trajectory of a NHSMP with backward and forward times.

The transition probabilities of a SMP change in function of the values of the recurrence time processes. In fact the conditional waiting times distribution functions (2.4) can be of any type, and then, also no memoryless distributions can be used. In this case, the time length spent in the starting state (initial backward value) changes the transition probabilities, and the same effect is determined by a starting forward which imposes a constraint on the time of the next transition. The same argument holds for the final backward and forward values.

To understand how the transition probabilities of the SMP are perturbed by imposing constraints on the recurrence processes both at starting and arriving times, we report here some of the results obtained in [8].

Definition 2.1. For all , and for all such that , one defines the transition probabilities with initial and final backward and forward

The following result holds

Theorem 2.2. For all , and for all such that , one has

Proof. See [8].

Definition 2.3. For all , and for all such that , one defines the probabilities with initial backward-certain forward and final backward-forward

Theorem 2.4. For all , and for all such that , one has

Proof. See [8].

There are several particular cases of (2.14) and (2.11) which can be of interest. In Section 3, some of them will be studied in relation to the credit-rating problem also in the asymptotic case.

3. Mono-Unireducible Non-Homogeneous Semi-Markov Process

In this section, we introduce the mono-unireducible topological structure. In [9], it was studied the asymptotic behavior of a mono-unireducible discrete time non-homogeneous SMP. Here, we extend the results to the continuous time case.

Definition 3.1. Let be a matrix.

is mono-unireducible if:

(i) is a circuit matrix,(ii) is a nonnegative column vector with at least one positive element,(iii).

Remark 3.2. The product of a mono-unireducible matrix with a positive diagonal matrix is a mono-unireducible matrix.

Definition 3.3. A NHCTSMP is said mono-unireducible if:

(i) ,(ii) is a mono-unireducible matrix.

As it is well known, one defines the -fold convolution of by itself as:

Lemma 3.4. Let be a NHCTSMP with a continuously differentiable and mono-unireducible kernel . Then, for all such that

where , .

Proof. Given , since is mono-unireducible, is a mono-unireducible matrix. Moreover, for all such that , because is continuous, , . By setting we have

and by iteration we get

The matrices are all mono-uniredudible, and then, their product is a mono-unireducible matrix and for is a full matrix as a consequence of Theorem 2.1. in [9].

At this point, we can state the following.

Theorem 3.5. Let be a NHSMP with continuously differentiable and mono-unireducible kernel . Then, for all

Proof. Let then,

Then, we get

In this way, for all , . Moreover which implies, for all the existence of .

Let suppose that , ; then for all

By considering the following inequalities:

inside (3.8), we get

Since, by Lemma 3.4, for , is a full matrix, we set and . Then, by substitution in (3.10), we have

which implies

Since is arbitrary, nothing can ensure (3.12) for mono-unireducible semi-Markov process and consequently for all , .

Corollary 3.6. Let be a NHSMP with continuously differentiable and mono-unireducible kernel . Then for all , , such that it results that

Proof. Because we are interested in the asymptotic behaviour as , it is not a loss of generality to consider that .

Let , then

Now, from Theorem 3.5 we know that

consequently if , otherwise for , we have that

Remark 3.7. From the previous theorem it results that for all and for all is the distribution function of going into the down state . Moreover , for all and for all is the distribution function of going into the down state given that at time the process is in state with a backward and forward values of and , respectively.

4. The Credit Risk Model

SMP were proposed to describe rating dynamics for the first time by [10] with the purpose of construct efficient rating migration models. Efficient migration models are of interest for pricing rating sensitive derivatives (see [11, 12]), for the valuation of portfolio of defaulting bonds, for credit risk management and capital allocation; see [13, 14].

Here, we apply our non-homogeneous model to entity ratings history, instrument ratings history, and issue/maturity ratings history, respectively for each standard & poor's rated entity, for each Standard & Poor's rated instrument stock or bonds sold by an entity at particular time, and issue/maturity for the Global Issuers and Structured Finance (GI and FS) instruments that Standard & Poor's has rated, in the GI&SF sector formerly known as the Corporate Finance since 1922 to 2007, July, 16.

In the Standard & Poor’s case there are different classes of rating expressing the creditworthiness of the rated firm. The ratings are listed to form the following set of states:

The creditworthiness is highest for the rating , assigned to firm extremely reliable with regard to financial obligations and decrease towards the rating which expresses the occurrence of payment default on some financial obligation. A table showing the financial meaning of the Standard & Poor’s rating categories is reported in [15]. As a matter of example, the rating is assigned to firm susceptible to economic conditions having still good credit quality.

The state space is partitioned in up and down states as follows:

To dispose of a semi-Markov model it is necessary to construct the embedded Markov chain and to find the d.f. of waiting times . The embedded Markov chain, constructed by rating agency Standard & Poor's (S&P) real data, was given in [16] and is reported in Tables 1, 2, and 3 for different years.

Table 1: Embedded Markov chain for .

AAA
AA
A
BBB
BB
B
CCC
D

AAA
0.623
0.319
0.043
0.015
0.000
0.000
0.000
0.000

AA
0.090
0.650
0.250
0.000
0.010
0.000
0.000
0.000

A
0.017
0.123
0.714
0.132
0.014
0.000
0.000
0.000

BBB
0.000
0.011
0.365
0.581
0.043
0.000
0.000
0.000

BB
0.000
0.000
0.000
0.500
0.249
0.187
0.064
0.000

B
0.000
0.199
0.000
0.200
0.000
0.000
0.601
0.000

CCC
0.000
0.000
0.000
0.000
0.143
0.286
0.143
0.428

D
0.000
0.000
0.000
0.000
0.000
0.000
0.000
1.000

Table 2: Embedded Markov chain for year .

AAA
AA
A
BBB
BB
B
CCC
D

AAA
0.612
0.274
0.064
0.020
0.007
0.017
0.004
0.002

AA
0.104
0.623
0.225
0.024
0.004
0.006
0.012
0.002

A
0.026
0.138
0.625
0.137
0.041
0.031
0.001
0.001

BBB
0.029
0.051
0.195
0.540
0.131
0.039
0.011
0.003

BB
0.037
0.006
0.025
0.187
0.481
0.218
0.043
0.003

B
0.011
0.021
0.015
0.047
0.147
0.473
0.242
0.044

CCC
0.000
0.000
0.000
0.000
0.028
0.113
0.634
0.225

D
0.000
0.000
0.000
0.000
0.000
0.000
0.000
1.000

Table 3: Embedded Markov chain for .

AAA
AA
A
BBB
BB
B
CCC
D

AAA
0.747
0.139
0.074
0.034
0.004
0.002
0.000
0.000

AA
0.052
0.573
0.309
0.036
0.026
0.002
0.001
0.001

A
0.024
0.180
0.539
0.240
0.011
0.003
0.002
0.001

BBB
0.007
0.039
0.228
0.461
0.229
0.033
0.001
0.002

BB
0.005
0.024
0.048
0.426
0.307
0.175
0.012
0.003

B
0.002
0.014
0.016
0.096
0.371
0.341
0.150
0.010

CCC
0.000
0.010
0.015
0.070
0.120
0.525
0.165
0.095

D
0.000
0.000
0.000
0.000
0.000
0.000
0.000
1.000

Remark 4.1. We would outline that the Markov matrices that are given yearly in the S&P publications are usually diagonal dominant. Instead, the matrices that are presented in this paper do not have this property. The reason is that in a semi-Markov environment the transitions are taken into account only if the rating agency assesses a new rating. On the contrary, in the evaluation of the S&P transition Markov chain, if in a year there was no new rating evaluation for a firm, it is supposed that the firm is yet in the same state, and the rating agency in the construction of the transition matrix takes into account that there is a “virtual” transition in the same state. This implies that the number of virtual transitions is very high and the Markov chain becomes diagonally dominant. We think that this is another reason of the best behaviour of the semi-Markov environment with respect to the Markov one, because only the real rating transitions are considered.

By using the empirical estimator of cumulative distribution functions, we estimated the distributions . In Figure 2, we report only few cases.

By using relation (2.5), we recover the semi-Markov kernel, and the application of the theoretical arguments of previous sections allow us to compute the transition probabilities. Following the approach in [2, 10] we denote by

the reliability function with initial backward. It expresses the probability that the firm has never defaulted up to time given that at time it was in rating class , where it entered with last transition at time .

It is simple to realize that is the probability to be defaulted at time .

In Figure 3, we show the reliability function for rating class for different values of , , and . As it is possible to see, the reliabilities have a different behaviour depending on the values of . This is due to the nonhomogeneity of the rating process. The reliability exhibits also variability as a function of the backward process, then as a result, our model assigns different survival probabilities to firms having the same rating class but with different age in this state.

(a)

(b)

(c)

(d)

Figure 2: Conditional waiting time d.f. (dotted lines) and (continuous lines) for (a), (b), (c), and (d).

(a)

(b)

(c)

(d)

Figure 3: Reliability function for (dotted lines) and (continuous lines) for (a), (b), and (b), (d).

References

J. Janssen and R. Manca, Semi-Markov Risk Models for Finance, Insurance and Reliability, Springer, New York, NY, USA, 2007.

G. D'Amico, J. Janssen, and R. Manca, “Initial and final backward and forward discrete time non-homogeneous semi-Markov credit risk models,” Methodology and Computing in Applied Probability, vol. 12, no. 2, pp. 215–225, 2010.

G. D'Amico, J. Janssen, and R. Manca, “Semi-Markov reliability models with recurrence times and credit rating applications,” Journal of Applied Mathematics and Decision Sciences, vol. 2009, Article ID 625712, 17 pages, 2009.

J. Yackel, “Limit theorems for semi-Markov processes,” Transactions of the American Mathematical Society, vol. 123, pp. 402–424, 1966.

E. Çinlar, “Markov renewal theory,” Advances in Applied Probability, vol. 1, pp. 123–187, 1969.

N. Limnios and G. Oprişan, Semi-Markov Processes and Reliability Modelling, Birkhäuser Boston Inc., Boston, Mass, USA, 2001.

V. Korolyuk and A. Swishchuk, Semi-Markov Random Evolutions, Kluwer Academic Publishers, Dodrecht, The Netherlands, 1995.

G. D'Amico, J. Janssen, and R. Manca, “Duration dependent semi-Markov models,” Applied Mathematical Sciences, vol. 5, no. 41–44, pp. 2097–2108, 2011.

G. D'Amico, J. Janssen, and R. Manca, “The dynamic behaviour of non-homogeneous single-unireducible Markov and semi- Markov chains,” in Networks: Topology and Dynamic Lectures Notes in Economic and Mathematical Systems, pp. 195–211, Springer, New York, NY, USA, 2009.

G. D'Amico, J. Janssen, and R. Manca, “Homogeneous semi-Markov reliability models for credit risk management,” Decisions in Economics and Finance, vol. 28, no. 2, pp. 79–93, 2005.

A. Vasileiou and P.-C. G. Vassiliou, “An inhomogeneous semi-Markov model for the term structure of credit risk spreads,” Advances in Applied Probability, vol. 38, no. 1, pp. 171–198, 2006.

G. D'Amico, J. Janssen, and R. Manca, “Valuing credit default swap in a non-homogeneous semi-Markovian rating based model,” Computational Economics, vol. 29, no. 2, pp. 119–138, 2007.

D. Lando, Credit Risk Modeling, Princeton University Press, Princeton, NJ, USA, 2004.

S. Trueck and S. T. Rachev, Rating Based Modeling of Credit Risk, Academic Press, New York, NY, USA, 2009.

C. Bluhm, L. Overbeck, and C. Wagner, An introduction to Credit Risk Modeling, CRC Financial Mathematics Series, Chapman & Hall, Boca Raton, Fla, USA, 2002.

G. D'Amico, G. Di Biase, J. Janssen, and R. Manca, “Homogeneous and non-homogeneous semi-Markov backward credit risk migration models,” in Financial Hedging, chapter 1, Nova Science Publishers, New York, NY, USA, 2009.