A nonparametric urn-based approach to interacting failing systems with an application to credit risk modeling

In this paper we propose a new nonparametric approach to interacting failing systems (FS), that is systems whose probability of failure is not negligible in a fixed time horizon, a typical example being firms and financial bonds. The main purpose when studying a FS is to calculate the probability of default and the distribution of the number of failures that may occur during the observation period. A model used to study a failing system is defined default model. In particular, we present a general recursive model constructed by the means of inter- acting urns. After introducing the theoretical model and its properties we show a first application to credit risk modeling, showing how to assess the idiosyncratic probability of default of an obligor and the joint probability of failure of a set of obligors in a portfolio of risks, that are divided into reliability classes

Moody's Correlated Binomial Default Distributions for Inhomogeneous Portfolios

This paper generalizes Moody's correlated binomial default distribution for homogeneous (exchangeable) credit portfolio, which is introduced by Witt, to the case of inhomogeneous portfolios. As inhomogeneous portfolios, we consider two cases. In the first case, we treat a portfolio whose assets have uniform default correlation and non-uniform default probabilities. We obtain the default probability distribution and study the effect of the inhomogeneity on it. The second case corresponds to a portfolio with inhomogeneous default correlation. Assets are categorized in several different sectors and the inter-sector and intra-sector correlations are not the same. We construct the joint default probabilities and obtain the default probability distribution. We show that as the number of assets in each sector decreases, inter-sector correlation becomes more important than intra-sector correlation. We study the maximum values of the inter-sector default correlation. Our generalization method can be applied to any correlated binomial default distribution model which has explicit relations to the conditional default probabilities or conditional default correlations, e.g. Credit Risk${}^{+}$, implied default distributions. We also compare some popular CDO pricing models from the viewpoint of the range of the implied tranche correlation.Comment: 29 pages, 17 figures and 1 tabl

Correlated Binomial Models and Correlation Structures

We discuss a general method to construct correlated binomial distributions by imposing several consistent relations on the joint probability function. We obtain self-consistency relations for the conditional correlations and conditional probabilities. The beta-binomial distribution is derived by a strong symmetric assumption on the conditional correlations. Our derivation clarifies the 'correlation' structure of the beta-binomial distribution. It is also possible to study the correlation structures of other probability distributions of exchangeable (homogeneous) correlated Bernoulli random variables. We study some distribution functions and discuss their behaviors in terms of their correlation structures.Comment: 12 pages, 7 figure

Correlation Structures of Correlated Binomial Models and Implied Default Distribution

We show how to analyze and interpret the correlation structures, the conditional expectation values and correlation coefficients of exchangeable Bernoulli random variables. We study implied default distributions for the iTraxx-CJ tranches and some popular probabilistic models, including the Gaussian copula model, Beta binomial distribution model and long-range Ising model. We interpret the differences in their profiles in terms of the correlation structures. The implied default distribution has singular correlation structures, reflecting the credit market implications. We point out two possible origins of the singular behavior.Comment: 16 pages, 7 figure

Infectious Default Model with Recovery and Continuous Limit

We introduce an infectious default and recovery model for N obligors. Obligors are assumed to be exchangeable and their states are described by N Bernoulli random variables S_{i} (i=1,...,N). They are expressed by multiplying independent Bernoulli variables X_{i},Y_{ij},Y'_{ij}, and default and recovery infections are described by Y_{ij} and Y'_{ij}. We obtain the default probability function P(k) for k defaults. Taking its continuous limit, we find two nontrivial probability distributions with the reflection symmetry of S_{i} \leftrightarrow 1-S_{i}. Their profiles are singular and oscillating and we understand it theoretically. We also compare P(k) with an implied default distribution function inferred from the quotes of iTraxx-CJ. In order to explain the behavior of the implied distribution, the recovery effect may be necessary.Comment: 13 pages, 7 figure

Nonlinear Parabolic Equations arising in Mathematical Finance

This survey paper is focused on qualitative and numerical analyses of fully nonlinear partial differential equations of parabolic type arising in financial mathematics. The main purpose is to review various non-linear extensions of the classical Black-Scholes theory for pricing financial instruments, as well as models of stochastic dynamic portfolio optimization leading to the Hamilton-Jacobi-Bellman (HJB) equation. After suitable transformations, both problems can be represented by solutions to nonlinear parabolic equations. Qualitative analysis will be focused on issues concerning the existence and uniqueness of solutions. In the numerical part we discuss a stable finite-volume and finite difference schemes for solving fully nonlinear parabolic equations.Comment: arXiv admin note: substantial text overlap with arXiv:1603.0387

Sampling from Archimedean copulas

We develop sampling algorithms for multivariate Archimedean copulas. For exchangeable copulas, where there is only one generating function, we first analyse the distribution of the copula itself, deriving a number of integral representations and a generating function representation. One of the integral representations is related, by a form of convolution, to the distribution whose Laplace transform yields the copula generating function. In the infinite-dimensional limit there is a direct connection between the distribution of the copula value and the inverse Laplace transform. Armed with these results, we present three sampling algorithms, all of which entail drawing from a one-dimensional distribution and then scaling the result to create random deviates distributed according to the copula. We implement and compare the various methods. For more general cases, in which an N-dimensional Archimedean copula is given by N-1 nested generating functions, we present algorithms in which each new variate is drawn conditional only on the value of the copula of the previously drawn variates. We also discuss the use of composite nested and exchangeable copulas for modelling random variates with a natural hierarchical structure, such as ratings and sectors for obligors in credit baskets.

Option hedging for small investors under liquidity costs

Following the framework of Cetin et al. (finance stoch. 8:311-341, 2004), we study the problem of super-replication in the presence of liquidity costs under additional restrictions on the gamma of the hedging strategies in a generalized black-scholes economy. We find that the minimal super-replication price is different from the one suggested by the black-scholes formula and is the unique viscosity solution of the associated dynamic programming equation. This is in contrast with the results of Cetin et al. (Finance Stoch. 8:311-341, 2004), who find that the arbitrage-free price of a contingent claim coincides with the Black-Scholes price. However, in Cetin et al. (Finance Stoch. 8:311-341, 2004) a larger class of admissible portfolio processes is used, and the replication is achieved in the L (2) approximating sense. JEL (C61 - G13 - D52)

Tracking bond indices in an integrated market and credit risk environment

The management of credit risky assets requires simulation models that integrate the disparate sources of credit and market risk, and suitable optimization models for scenario analysis. In this paper we integrate Monte Carlo simulation models for credit risk with scenario optimization, and develop a methodology for tracking broadly defined corporate bond indices. Testing of the models shows that the integration of the multiple risk factors improves significantly the performance of tracking models. Good tracking performance can be achieved by optimizing strategic asset allocation among broad classes of corporate bonds. However, extra value is generated with a tactical model that optimizes bond picking decisions as well. It is also shown that adding small corporate bond holdings in portfolios that track government bond indices improves the risk/return characteristics of the portfolios. The empirical results to substantiate the findings of this study are obtained by backtesting the model over a recent 30�month period.