A Markov Chain Copula Model for Credit Default Swaps with Bilateral Counterparty Risk

Reduced-form credit risk models are widely used in pricing and hedging credit derivatives. Generating default dependency is the key element in any such model. In this article, we use Markov copulae approach to model the dependence structure of defaults between the three obligors, one is the reference entity, another is the protection seller, the other is the protection buyer(the investor), so we can consider the bilateral counterparty risk of credit default swaps(CDS). In this Markov chain copula model, we obtain the explicit formulas of the CDS premium rates C ₁(T) (with unilateral counterparty risk) and C ₂(T) (with bilateral counterparty risk). And then we perform some numerical experiments to analyze the difference of the fair spreads between the unilateral case and the bilateral case.     

Some Results for Uniformizable Semi-Markov Processes

It is well known that finite Markov chains (M.C.s) in continuous time are uniformizable. That is, a finite M. C. in continuous time can be treated as an M. C. in discrete time with random Poisson transition epochs. In this paper, we see to what extent generalization of the uniformization to a class of semi-Markov Processes (S.M.P.s) is possible. A necessary condition under which S.M.P.s are uniformizable is provided. It is shown that, an S.M.P. with dwell-time distributions depending only on the current state is uniformizable if and only if the distributions are compound geometric distributions having the same base distribution. It is also shown that if the distributions are of generalized phase type then an S.M.P. being uniformizable implies that it is an M.C. in continuous time. Some properties that are shared by a uniformizable S.M.P. and the associated M.C. in discrete time are also discussed.     

基于贝叶斯方法的信用风险损失分布研究

现代商业银行进行经济资本配置时,采用的损失分布函数都存在严重的失真问题。运用贝叶斯方法,充分利用各种信息对正态分布形式的信用损失分布进行了修正,得到信用风险损失分布的优化模型,结果表明：修正后的信用风险损失分布具有较高的精度,从而为商业银行经济资本管理提供了一种很实用的管理工具。     

信用风险、融资担保和动态保费估值

基于市场化信息和风险中性的概念,度量了有担保贷款的边际违约概率和累积违约概率,确定了贷款担保风险的精算现值,根据市场信用价差的变化给出了动态保费每期的调整幅度,并利用数值模拟进行了担保费率的比较静态分析,最后根据实际的担保数据给出了动态保费的实证检验.结果显示,实际的违约支付非常接近于动态保费的估值,证明动态保费估值模型是一个简单、可行和实用的定价模型.     

知识共享对供应链信用风险传染的影响研究

将企业间知识共享这一因素引入供应链上下游企业博弈过程,在完全信息条件下研究三级供应链中零售商分别与分销商、制造商进行知识共享时知识共享量与企业信用风险水平和供应链信用风险传染强度之间的关系,并结合算例进行了数值分析.基于Stackelberg博弈模型的研究结果表明,供应链上的知识共享有助于降低供应链成员企业的信用风险;进行知识共享的企业间信用风险传染强度随着知识共享量的增加而增大;企业进行知识共享比不进行知识共享时受到其它企业信用风险传染的影响大.     

基于MCMC模拟的贝叶斯分层信用风险评估模型

缺少违约数据与债务人异质性是度量信用风险时面临的重要问题。贝叶斯模型中分层先验信息和马尔可夫链蒙特卡罗（MCMC）模拟方法的应用可以有效缓解数据缺失和测量误差问题,并能对债务人异质性进行评价和比较,从而避免低估风险。针对银行数据的模型拟合与模型诊断均展现了分层估计的适应性和灵活性,相关方法简洁清晰,利于国内风险分析人员采用。同时,涵盖宏观经济协变量的贝叶斯分层模型可以用于更加复杂的风险分析。     

中国P2P网络借贷信用风险的测量

分析中国P2P网络借贷信用风险的生成机理,借鉴商业银行典型个人信用风险测量理论基础与方法,构建基于涵盖借款人软信息的P2P网络借贷信用风险评级指标体系。通过向在银行从事信贷工作的专业人员发放调查问卷,确定具体指标,使用AHP(层次分析法)与DEMATEL(决策实验室法)相结合的方法综合确定指标体系权重,最后评级结果以百分制数值给出,并换算成具体的信用等级,有效解决了平台对借款人信用等级测量问题。     

Limit Theorems for Semi-Markov Queues and Their Applications

We begin by definition of semi-Markov flow and discussion of its properties. Asymptotic behavior of multi-server and single-server queueing systems is studied under assumption of time-compression or service time growth. The results obtained are used for calculation of large systems reliability. Statistical estimates of parameters involved are also provided.     

Goodness-of-fit Tests for Semi-Markov and Markov Survival Models with One Intermediate State

Survival data with one intermediate state are described by semi-Markov and Markov models for counting processes whose intensities are defined in terms of two stopping times T ₁< T ₂. Problems of goodness-of-fit for these models are studied. The test statistics are proposed by comparing Nelson–Aalen estimators for data stratified according to T ₁. Asymptotic distributions of these statistics are established in terms of the weak convergence of some random fields. Asymptotic consistency of these test statistics is also established. Simulation studies are included to indicate their numerical performance.     

Simulation/Regression Pricing Schemes for CVA Computations on CDO Tranches

We devise simulation/regression numerical schemes for pricing the CVA on CDO tranches, where CVA stands for Credit Valuation Adjustment, or price correction accounting for the defaultability of a counterparty in an OTC derivatives transaction. This is done in the setup of a continuous-time Markov chain model of default times, in which dependence between credit names is represented by the possibility of simultaneous defaults. The main idea of this article is to perform the nonlinear regressions which are used for computing conditional expectations, in the time variable for a given state of the model, rather than in the space variables at a given time in diffusive setups. This idea is formalized as a lemma which is valid in any continuous-time Markov chain model. It is then implemented on the targeted application of CVA computations on CDO tranches.     

A Bottom-Up Dynamic Model of Portfolio Credit Risk with Stochastic Intensities and Random Recoveries

In Bielecki et al. (2014a Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014a ). Dynamic hedging of portfolio credit risk in a markov copula model . J. Optimiz. Theor. Applic . doi: DOI 10.1007/s10957-013-0318-4 (forthcoming) .[Crossref] , [Google Scholar]), the authors introduced a Markov copula model of portfolio credit risk where pricing and hedging can be done in a sound theoretical and practical way. Further theoretical backgrounds and practical details are developed in Bielecki et al. (2014b Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014b ). A bottom-up dynamic model of portfolio credit risk - Part I: Markov copula perspective . In: Recent Adv. Fin. Eng. 2012 , World Scientific (preprint version available at http://dx.doi.org/10.2139/ssrn.1844574) . [Google Scholar],c) where numerical illustrations assumed deterministic intensities and constant recoveries. In the present paper, we show how to incorporate stochastic default intensities and random recoveries in the bottom-up modeling framework of Bielecki et al. (2014a Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014a ). Dynamic hedging of portfolio credit risk in a markov copula model . J. Optimiz. Theor. Applic . doi: DOI 10.1007/s10957-013-0318-4 (forthcoming) .[Crossref] , [Google Scholar]) while preserving numerical tractability. These two features are of primary importance for applications like CVA computations on credit derivatives (Assefa et al., 2011 Assefa , S. , Bielecki , T. R. , Crépey , S. , Jeanblanc , M. ( 2011 ). CVA computation for counterparty risk assessment in credit portfolios . In: Bielecki , T.R. , Brigo , D. , Patras , F. , Eds., Credit Risk Frontiers . Hoboken : Wiley/Bloomberg-Press . [Google Scholar]; Bielecki et al., 2012 Bielecki , T. R. , Crépey , S. , Jeanblanc , M. , Zargari , B. ( 2012 ). Valuation and Hedging of CDS counterparty exposure in a markov copula model . Int. J. Theoret. Appl. Fin. 15 ( 1 ): 1250004 .[Crossref] , [Google Scholar]), as CVA is sensitive to the stochastic nature of credit spreads and random recoveries allow to achieve satisfactory calibration even for “badly behaved” data sets. This article is thus a complement to Bielecki et al. (2014a Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014a ). Dynamic hedging of portfolio credit risk in a markov copula model . J. Optimiz. Theor. Applic . doi: DOI 10.1007/s10957-013-0318-4 (forthcoming) .[Crossref] , [Google Scholar]), Bielecki et al. (2014b Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014b ). A bottom-up dynamic model of portfolio credit risk - Part I: Markov copula perspective . In: Recent Adv. Fin. Eng. 2012 , World Scientific (preprint version available at http://dx.doi.org/10.2139/ssrn.1844574) . [Google Scholar]) and Bielecki et al. (2014c Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014c ). A bottom-up dynamic model of portfolio credit risk - Part II: Common-shock interpretation, calibration and hedging issues . Recent Adv. Fin. Eng. 2012 , World Scientific (preprint version available at http://dx.doi.org/10.2139/ssrn.2245130) . [Google Scholar]).     

An EM and a Stochastic Version of the EM Algorithm for Nonparametric Hidden Semi-Markov Models

The Hidden semi-Markov models (HSMMs) were introduced to overcome the constraint of a geometric sojourn time distribution for the different hidden states in the classical hidden Markov models. Several variations of HSMMs were proposed that model the sojourn times by a parametric or a nonparametric family of distributions. In this article, we concentrate our interest on the nonparametric case where the duration distributions are attached to transitions and not to states as in most of the published papers in HSMMs. Therefore, it is worth noticing that here we treat the underlying hidden semi-Markov chain in its general probabilistic structure. In that case, Barbu and Limnios (2008 Barbu , V. , Limnios , N. ( 2008 ). Semi-Markov Chains and Hidden Semi-Markov Models Toward Applications: Their Use in Reliability and DNA Analysis . New York : Springer . [Google Scholar]) proposed an Expectation–Maximization (EM) algorithm in order to estimate the semi-Markov kernel and the emission probabilities that characterize the dynamics of the model. In this article, we consider an improved version of Barbu and Limnios' EM algorithm which is faster than the original one. Moreover, we propose a stochastic version of the EM algorithm that achieves comparable estimates with the EM algorithm in less execution time. Some numerical examples are provided which illustrate the efficient performance of the proposed algorithms.     

Bivariate frailty model for the analysis of multivariate survival time

Because of limitations of the univariate frailty model in analysis of multivariate survival data, a bivariate frailty model is introduced for the analysis of bivariate survival data. This provides tremendous flexibility especially in allowing negative associations between subjects within the same cluster. The approach involves incorporating into the model two possibly correlated frailties for each cluster. The bivariate lognormal distribution is used as the frailty distribution. The model is then generalized to multivariate survival data with two distinguished groups and also to alternating process data. A modified EM algorithm is developed with no requirement of specification of the baseline hazards. The estimators are generalized maximum likelihood estimators with subject-specific interpretation. The model is applied to a mental health study on evaluation of health policy effects for inpatient psychiatric care.     

Credit risk assessment with Bayesian model averaging

Many credit risk models are based on the selection of a single logistic regression model, on which to base parameter estimation. When many competing models are available, and without enough guidance from economical theory, model averaging represents an appealing alternative to the selection of single models. Despite model averaging approaches have been present in statistics for many years, only recently they are starting to receive attention in economics and finance applications. This contribution shows how Bayesian model averaging can be applied to credit risk estimation, a research area that has received a great deal of attention recently, especially in the light of the global financial crisis of the last few years and the correlated attempts to regulate international finance. The paper considers the use of logistic regression models under the Bayesian Model Averaging paradigm. We argue that Bayesian model averaging is not only more correct from a theoretical viewpoint, but also slightly superior, in terms of predictive performance, with respect to single selected models.     

Bivariate location–scale models for regression analysis, with applications to lifetime data

Summary.  The literature on multivariate linear regression includes multivariate normal models, models that are used in survival analysis and a variety of models that are used in other areas such as econometrics. The paper considers the class of location–scale models, which includes a large proportion of the preceding models. It is shown that, for complete data, the maximum likelihood estimators for regression coefficients in a linear location–scale framework are consistent even when the joint distribution is misspecified. In addition, gains in efficiency arising from the use of a bivariate model, as opposed to separate univariate models, are studied. A major area of application for multivariate regression models is to clustered, 'parallel' lifetime data, so we also study the case of censored responses. Estimators of regression coefficients are no longer consistent under model misspecification, but we give simulation results that show that the bias is small in many practical situations. Gains in efficiency from bivariate models are also examined in the censored data setting. The methodology in the paper is illustrated by using lifetime data from the Diabetic Retinopathy Study.     

Continuous Time Wishart Process for Stochastic Risk

Risks are usually represented and measured by volatility–covolatility matrices. Wishart processes are models for a dynamic analysis of multivariate risk and describe the evolution of stochastic volatility–covolatility matrices, constrained to be symmetric positive definite. The autoregressive Wishart process (WAR) is the multivariate extension of the Cox, Ingersoll, Ross (CIR) process introduced for scalar stochastic volatility. As a CIR process it allows for closed-form solutions for a number of financial problems, such as term structure of T-bonds and corporate bonds, derivative pricing in a multivariate stochastic volatility model, and the structural model for credit risk. Moreover, the Wishart dynamics are very flexible and are serious competitors for less structural multivariate ARCH models.     

Testing for Bivariate Extreme Dependence Using Kendall's Process

Abstract. New tests for the hypothesis of bivariate extreme‐value dependence are proposed. All test statistics that are investigated are continuous functionals of either Kendall's process or its version with estimated parameters. The procedures considered are based on linear combinations of moments and on Cramér–von Mises distances. A suitably adapted version of the multiplier central limit theorem for Kendall's process enables the computation of asymptotically valid p‐values. The power of the tests is evaluated for small, moderate and large sample sizes, as well as asymptotically, under local alternatives. An illustration with a real data set is presented.     

First hitting probabilities for semi-Markov chains and estimation

We first consider a stochastic system described by an absorbing semi-Markov chain (SMC) with finite state space, and we introduce the absorption probability to a class of recurrent states. Afterwards, we study the first hitting probability to a subset of states for an irreducible SMC. In the latter case, a non-parametric estimator for the first hitting probability is proposed and the asymptotic properties of strong consistency and asymptotic normality are proven. Finally, a numerical application on a five-state system is presented to illustrate the performance of this estimator.     

Convergence Behavior of the em algorithm for the multivariate t -distribution

Iterative reweighting (IR) is a popular method for computing M-estimates of location and scatter in multivariate robust estimation. When the objective function comes from a scale mixture of normal distributions the iterative reweighting algorithm can be identified as an EM algorithm. The purpose of this paper is to show that in the special case of the multivariate t-distribution, substantial improvements to the convergence rate can be obtained by modifying the EM algorithm.     

A Bivariate Gamma Model for a Latent Degradation Process

A model is presented in this article based on a bivariate gamma process in which, the first component is latent and determines the failure time and the second represents the marker. This process is a more realistic model for a degradation process. After introducing the model, we obtain failure and survival probability distributions and discuss parametric and predictive inferences based on the Maximum Likelihood method and in a Bayesian setup.