R语言用CPV模型的房地产信贷信用风险的度量和预测

复制粘贴侠
• 阅读 801

全文链接:http://tecdat.cn/?p=30401

原文出处:拓端数据部落公众号

本文基于 CPV 模型, 对房地产信贷风险进行了度量与预测。我们被客户要求撰写关于CPV模型的研究报告。结果表明, 该模型在度量和预测房地产信贷违约率方面具有较好的效果。

CPV 模型的基本原理和框架

CPV 模型是一个用于分析贷款组合风险和收益的多因素模型, 它根据失业率、长期利率、GDP 增长率、汇率、政府支出和总储蓄率等宏观因素, 借助于经济计量工具和蒙特卡罗技术, 将每个国家不同行业中不同等级的违约概率和转移概率相联系, 进而计算出风险价值( 迪迪埃、皮罗特, 2005) 。

指标及样本数据选择

当实际运用 CPV 模型时, 宏观经济因素个数必须达到 3 个以上时该模型才具有一定的有效性 ( 估计有效性及预测有效性) 。因此, 本文分别从国家宏观经济、房地产行业状况、房地产企业状况三个方面选择出三个宏观经济因素指标, 运用 CPV 模型评估房地产信贷的信用风险。1. 综合领先指标(Composite Leading Indicator) 。经济合作与发展组织 (OECD) 的综合领先指标(CompositeLeading Indicator, CLI) 被认为是预测全球经济变动趋势的良好指标, 它是指一系列引导经济由增长至衰退的循环的相关经济指标和经济变量的加权平均数, 主要用来预测整体经济的转变情况, 预测未来数月的经济发展趋势:

###违约率再0到1之间 ,因此对其进行标准化

cr=(fbzdata$NA.-min(fbzdata$NA.))/(max(fbzdata$NA.)-min(fbzdata$NA.))

##绘制时间序列图

plot(diff(cr),type="l")

 R语言用CPV模型的房地产信贷信用风险的度量和预测

从上图可以看到数据呈现出稳定的波动趋势,再均值周围回摆,因此差分后数列稳定。

CPV模型

#逻辑回归结果

 

##滞后期的选择

#MA(1)

ma1=arma(x=cr,order=c(0,1))
summary(ma4)#查看相关系数,Std. Error  t value,AIC值

滞后期的选择根据Akaike Information Criterion(AIC)准则确定。

 R语言用CPV模型的房地产信贷信用风险的度量和预测

根据回归结果显示CLI、CRECI和ECI的P值显示的数值表明该估计的有效性。 从估计出的敏感系数可以看出,CLI、和ECI的系数(Coefficient)为正, CRECI的系数(Coefficient)为负,并且ECI的系数值是CLI系数值的2倍左右所以, 当CRECI的值上升的时候,Y值也上升 (违约率下降), 它们之间是反向变化。并且, 综合领先指标ECI相对于其他宏观经济指标来说, 影响力度最大。

残差值分析

plot(CPVmod$residuals,type="l")#残差值

lines(CPVmod$fitted.values,type="l")#拟合值

lines(cr,type="l")#实际值

abline(h=0)#添加0基准线

 R语言用CPV模型的房地产信贷信用风险的度量和预测

根据上面的模型估计结果, 可以得出Y的残差值、实际值和拟合值的趋势图。如图1所示。

根据图1的趋势图, 可以清楚地看出, Y的拟合值和实际值的曲线几乎完全重合, 这也说明该模型很好地拟合了样本数据, 也跟上述结果表明的事实相符合。

另外, 残差的相关系数和偏相关系数如图2所表示。从图2可以看出, 模型的残差不存在序列自相关。

残差自相关和偏相关分析

acf(CPVmod$residuals)#自相关系数

pacf(CPVmod$residuals)#偏相关系数

 R语言用CPV模型的房地产信贷信用风险的度量和预测

 R语言用CPV模型的房地产信贷信用风险的度量和预测

BG高阶自相关系数检验

 R语言用CPV模型的房地产信贷信用风险的度量和预测

进一步用布—戈弗雷检验(Breusch-Godfrey, BG检验)来检验残差的序列自相关性的情况, 从上面的检验结果,可以看出, P值较大, 根据BG高阶自相关系数检验原理, 该检验结果接受原假设, 即上述模型的残差不存在自相关性, 也进一步证实了前面的自相关系数和偏自相关系数的检验结果。

综上所述, 此模型检验达到预期要求, 较好地反映了宏观经济实际状况和各宏观经济变量之间的对应关系。

宏观经济变量的趋势


lines(fbzdata$CLI,col="red",type="l")

 

lines(fbzdata$CERCI,col="green",type="l")

 R语言用CPV模型的房地产信贷信用风险的度量和预测

违约率拟合值和实际值的对比

plot(CPVmod$fitted.values,type="l")#拟合值

lines(cr,type="l")#实际值

 R语言用CPV模型的房地产信贷信用风险的度量和预测

预测信贷违约率

pre=predict(ma7, 50)#预测值

plot(unlist(pre),type="l")

 R语言用CPV模型的房地产信贷信用风险的度量和预测

仍选择综合领先指标(CLI), 国房景气指数(CRECI)以及企业景气指数(ECI)三个宏观经济变量, 运用CPV模型的估计结果对未来2年的每个月的房地产信贷违约率进行预测,

为了直观地表示出违约率的变化趋势及其与CLI、CREIC、ECI三个变量之间的关系, 分别做出趋势图。从中可以看出:在11月份, 虽然ECI上升, 但是由于CRECI和CLI下降的共同作用, 对于房地产信贷来说经济状况趋势还是下降, 因此违约率略有上升。但是从3月份开始, 由于宏观经济、房地产行业和企业三个层面的景气好转, 房地产信贷违约率逐渐下降。这个预测趋势结果仅作为参考, 这里主要是为了说明模型的预测功能对于商业银行抵抗信贷风险的重要意义。


 R语言用CPV模型的房地产信贷信用风险的度量和预测

最受欢迎的见解

1.在python中使用lstm和pytorch进行时间序列预测

2.python中利用长短期记忆模型lstm进行时间序列预测分析

3.Python用RNN循环神经网络:LSTM长期记忆、GRU门循环单元、回归和ARIMA对COVID-19新冠疫情新增人数时间序列

4.Python TensorFlow循环神经网络RNN-LSTM神经网络预测股票市场价格时间序列和MSE评估准确性

5.r语言copulas和金融时间序列案例

6.R 语言用RNN循环神经网络 、LSTM长短期记忆网络实现时间序列长期利率预测

7.Matlab创建向量自回归(VAR)模型分析消费者价格指数 (CPI) 和失业率时间序列

8.r语言k-shape时间序列聚类方法对股票价格时间序列聚类

9.R语言结合新冠疫情COVID-19股票价格预测:ARIMA,KNN和神经网络时间序列分析

点赞
收藏
评论区
推荐文章
blmius blmius
3年前
MySQL:[Err] 1292 - Incorrect datetime value: ‘0000-00-00 00:00:00‘ for column ‘CREATE_TIME‘ at row 1
文章目录问题用navicat导入数据时,报错:原因这是因为当前的MySQL不支持datetime为0的情况。解决修改sql\mode:sql\mode:SQLMode定义了MySQL应支持的SQL语法、数据校验等,这样可以更容易地在不同的环境中使用MySQL。全局s
皕杰报表之UUID
​在我们用皕杰报表工具设计填报报表时,如何在新增行里自动增加id呢?能新增整数排序id吗?目前可以在新增行里自动增加id,但只能用uuid函数增加UUID编码,不能新增整数排序id。uuid函数说明:获取一个UUID,可以在填报表中用来创建数据ID语法:uuid()或uuid(sep)参数说明:sep布尔值,生成的uuid中是否包含分隔符'',缺省为
Jacquelyn38 Jacquelyn38
4年前
2020年前端实用代码段,为你的工作保驾护航
有空的时候,自己总结了几个代码段,在开发中也经常使用,谢谢。1、使用解构获取json数据let jsonData  id: 1,status: "OK",data: 'a', 'b';let  id, status, data: number   jsonData;console.log(id, status, number )
Stella981 Stella981
3年前
PhoneGap设置Icon
参考:http://cordova.apache.org/docs/en/latest/config\_ref/images.html通过config.xml中的<icon标签来设置Icon<iconsrc"res/ios/icon.png"platform"ios"width"57"height"57"densi
Stella981 Stella981
3年前
KVM调整cpu和内存
一.修改kvm虚拟机的配置1、virsheditcentos7找到“memory”和“vcpu”标签,将<namecentos7</name<uuid2220a6d1a36a4fbb8523e078b3dfe795</uuid
Wesley13 Wesley13
3年前
mysql设置时区
mysql设置时区mysql\_query("SETtime\_zone'8:00'")ordie('时区设置失败,请联系管理员!');中国在东8区所以加8方法二:selectcount(user\_id)asdevice,CONVERT\_TZ(FROM\_UNIXTIME(reg\_time),'08:00','0
Wesley13 Wesley13
3年前
00:Java简单了解
浅谈Java之概述Java是SUN(StanfordUniversityNetwork),斯坦福大学网络公司)1995年推出的一门高级编程语言。Java是一种面向Internet的编程语言。随着Java技术在web方面的不断成熟,已经成为Web应用程序的首选开发语言。Java是简单易学,完全面向对象,安全可靠,与平台无关的编程语言。
Stella981 Stella981
3年前
Django中Admin中的一些参数配置
设置在列表中显示的字段,id为django模型默认的主键list_display('id','name','sex','profession','email','qq','phone','status','create_time')设置在列表可编辑字段list_editable
Wesley13 Wesley13
3年前
MySQL部分从库上面因为大量的临时表tmp_table造成慢查询
背景描述Time:20190124T00:08:14.70572408:00User@Host:@Id:Schema:sentrymetaLast_errno:0Killed:0Query_time:0.315758Lock_
Python进阶者 Python进阶者
1年前
Excel中这日期老是出来00:00:00,怎么用Pandas把这个去除
大家好,我是皮皮。一、前言前几天在Python白银交流群【上海新年人】问了一个Pandas数据筛选的问题。问题如下:这日期老是出来00:00:00,怎么把这个去除。二、实现过程后来【论草莓如何成为冻干莓】给了一个思路和代码如下:pd.toexcel之前把这
美凌格栋栋酱 美凌格栋栋酱
5个月前
Oracle 分组与拼接字符串同时使用
SELECTT.,ROWNUMIDFROM(SELECTT.EMPLID,T.NAME,T.BU,T.REALDEPART,T.FORMATDATE,SUM(T.S0)S0,MAX(UPDATETIME)CREATETIME,LISTAGG(TOCHAR(