CN109961362A - P2p平台信用风险动态评价方法和系统 - Google Patents
P2p平台信用风险动态评价方法和系统 Download PDFInfo
- Publication number
- CN109961362A CN109961362A CN201910127901.XA CN201910127901A CN109961362A CN 109961362 A CN109961362 A CN 109961362A CN 201910127901 A CN201910127901 A CN 201910127901A CN 109961362 A CN109961362 A CN 109961362A
- Authority
- CN
- China
- Prior art keywords
- platform
- probability
- submodel
- time
- default
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q40/00—Finance; Insurance; Tax strategies; Processing of corporate or income taxes
- G06Q40/03—Credit; Loans; Processing thereof
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Accounting & Taxation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Finance (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Mathematical Physics (AREA)
- Evolutionary Biology (AREA)
- General Engineering & Computer Science (AREA)
- Probability & Statistics with Applications (AREA)
- Development Economics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Software Systems (AREA)
- Economics (AREA)
- Marketing (AREA)
- Strategic Management (AREA)
- Technology Law (AREA)
- General Business, Economics & Management (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种P2P平台信用风险动态评价方法和系统,涉及互联网金融技术领域。其包括:S1、获取P2P平台的信用数据;S2、构建基于混合生存分析的动态评价模型;S3、利用面板数据估计模型参数;S4、预测P2P平台动态违约概率。本发明的一种P2P平台信用风险动态评价方法和系统,通过分别构建违约状态估计子模型和违约时间估计子模型来预测平台是否会发生违约以及何时会发生违约。在构建动态评价模型的同时,还设计了一种多重生存观测构建方法,使模型能够有效学习到平台信用特征面板数据中的有效信息,从而提升模型的预测效果,准确地预测平台的风险趋势,为P2P平台提供及时有效的内部风险控制,同时为投资人提供投资决策支持。
Description
技术领域
本发明涉及互联网金融技术领域,具体涉及一种P2P平台信用风险动态评价方法和系统。
背景技术
P2P(Peer-to-Peer)借贷指个体和个体间通过互联网P2P平台实现直接的借贷。P2P借贷作为新型的互联网金融模式,利用互联网为主要渠道,为借款人和投资人提供信息交互与信贷撮合等服务,具有门槛低,方便,快捷等优点。从2007年首次引入我国开始,P2P借贷得到飞速发展,P2P平台的数量和P2P借贷交易量都呈现大幅增长,极大程度弥补了传统借贷模式未能有效针对个体借款人提供小额借贷服务的不足。
但由于缺乏健全的法律法规、成熟的监管手段以及高度的信息不对称性,P2P平台跑路、经侦介入、提现困难以及网站关闭等风险事件频频发生。由于P2P网络借贷行业及信用风险问题频发,因此很多用户在使用P2P平台之前,都会对P2P平台的信用风险进行评价。目前对P2P平台的信用风险评价,一般是沿用传统的信用评分模型或通过专家评估。
然而,传统的信用评分方法是基于分类模型的违约风险预测,只能得到表现期或整个借款周期内的单一的违约概率,P2P平台的信用风险是变化的过程,仅仅依靠表现期或整个借款周期内的单一的违约概率无法准确评价P2P平台的风险趋势。专家评估则会受到人为因素干扰,导致评估结果出现误差。
发明内容
(一)解决的技术问题
针对现有技术的不足,本发明提供了一种P2P平台信用风险动态评价方法和系统,解决了现有技术中对P2P平台信用风险评估不准确的技术问题。
(二)技术方案
为实现以上目的,本发明通过以下技术方案予以实现:
本发明提供了一种P2P平台信用风险动态评价方法,该方法包括以下步骤:
S1、获取P2P平台的信用数据,所述信用数据包括:观测时间、观测时间内的信用特征变量、观测时间内的删失状态数据以及观测时间内的违约状态数据;所述观测时间内的删失状态数据表示在观测时间内未能观测到违约;
S2、构建基于混合生存分析的动态评价模型,包括:
S201、定义P2P平台的信用状态,包括定义所述步骤S1中的所述观测时间、所述删失状态数据以及所述违约状态数据;
S202、基于S201的信用状态构建P2P平台信用风险的动态评价模型;
S203、在S202的基础上构建违约状态估计子模型;
S204、在S202的基础上构建违约时间估计子模型;
S3、利用面板数据估计模型参数,包括:
S301、步骤S2中的模型在不同观测时间下,可观测到构建违约状态估计子模型以及构建违约时间估计子模型的信用特征变量,构成了P2P平台信用特征面板数据,针对不同观测时间,构建多重生存观测数据;
S302、基于步骤S301中所构建的多重生存观测数据,利用EM算法迭代估计步骤S2中动态评价模型的参数;
S4、预测P2P平台动态违约概率,通过步骤S1—步骤S3,得到训练好的动态评价模型,针对预测样本,预测P2P平台在不同时间下的生存概率;基于P2P平台在不同时间下的生存概率,可计算在不同时间下的违约概率;
其中,在步骤S202中:所述动态评价模型为:P2P平台在指定观察时间下的生存概率等于P2P平台非违约的概率与P2P平台在指定观察时间之后发生违约的概率之和;
在步骤S203中:所述违约状态估计子模型为:针对P2P平台的违约概率,利用逻辑回归构建违约状态估计子模型;
在步骤S204中:所述违约时间估计子模型为:针对P2P平台的条件生存概率,利用Cox比例风险回归构建违约时间估计子模型。
优选的,在步骤S201中,所述定义P2P平台的信用状态:
定义变量t为观测时间;
定义随机变量δ为删失状态数据指示变量,δ=0表示观测数据删失,即在时间t内未能观测到违约,δ=1表示观测数据未删失,即在t时间内观测到违约;
定义随机变量y为违约状态数据指示变量,y=0与表示P2P平台不会发生违约,y=1表示P2P平台会发生违约。
优选的,在步骤S202中,所述P2P平台信用风险的动态评价模型为:
S(t)=1-p+p*S(t|y=1)
其中:
S(t)为目标函数,是P2P平台在时间t的生存概率;
p=P(y=1)是P2P平台的违约概率;
S(t|y=1)=P(T>t|y=1)是假定P2P平台会发生违约下的条件生存函数;
目标函数S(t)=1-p+p*S(t|y=1)表达的含义是:P2P平台在时间t下的生存概率等于P2P平台非违约的概率与P2P平台在时间t之后发生违约的概率之和。
优选的,在步骤S203中,所述违约状态估计子模型为:
其中:
z=(1,z1,...,zn)为违约状态估计子模型的信用特征变量;
β=(β0,β1,...,βn)为信用特征变量的系数。
优选的,在步骤S204中,所述违约时间估计子模型为:
其中:
S0(t|y=1)为基准生存函数,是函数的非参数部分;
向量x=(x1,x2,…,xm)为违约时间估计子模型的信用特征变量;
向量α=(α1,α2,...,αm)为信用特征变量的系数。
优选的,在所述步骤S302中,包括:
步骤S3021:
利用删失变量作为初始违约概率期望值E(yi)(0)=δi;
估计违约概率估计子模型的初始参数β(0)以及违约时间估计子模型的初始参数α(0),S0 (0);
设定模型迭代的最大次数C,模型收敛函数g(·)=+∞以及收敛阈值G,迭代次数指示变量k=1;
若迭代次数k<C,且收敛函数g(·)>G,则进行以下迭代;
步骤S3022:
利用违约概率期望值E(yi)(k-1)计算违约概率估计子模型的极大似然估计函数:
其中:
LI为违约概率估计子模型的极大似然函数;
p(zi)为违约状态估计子模型结果;
得到该轮迭代下违约概率估计子模型的参数估计值:
β(k)=(β1,β2,…,βn)(k)
步骤S3023:
利用违约概率期望值E(yi)(k-1)计算违约时间估计子模型的极大似然估计函数:
其中:
LL为违约时间估计子模型的极大似然函数;
h(ti|y=1,xi)为条件生存函数S(ti|y=1,xi)所对应的风险函数;
S(ti|y=1,xi)是假定平台会发生违约下的条件生存函数;
得到该轮迭代下违约概率估计子模型的参数估计值:
α(k)=(α1,α2,…,αm)(k)
步骤S3024:
计算违约概率yi的期望值:
其中:
E(k)(yi|ti,δi,zi,xi)为k轮迭代下违约概率yi的期望值;
步骤S3025:
计算迭代收敛函数:
g(·)=||β(k)-β(k-1)||2 2+||α(k)-α(k-1)||2 2+||S0 (k)-S0 (k-1)||2 2
其中:
g(·)为迭代收敛函数;
β(k)为k轮迭代下违约状态估计子模型的参数值;
α(k),S0 (k)为k轮迭代下违约时间估计子模型的参数值;
步骤S3026:
更新迭代次数:
k=k+1
步骤S3022至S3026一直循环迭代,直至迭代达到最大次数(k=C)或迭代收敛函数小于设定的阈值(g(k)≤G);
得到动态评价函数的参数估计值:
β(*)=(β0,β1,...,βn)(*)
α(*)=(α1,α2,…,αm)(*)
S0 (*)(t|y=1)
其中:
β(*),α(*),S0 (*)(t|y=1)分别为迭代结束时违约状态估计子模型和违约时间估计子模型的参数值。
优选的,在步骤S4中,
预测P2P平台在不同时间下的生存概率表示为:
其表示:针对预测样本,预测P2P平台在不同时间下的生存概率;
违约概率的表示方法为:
P(t|xv,zv)=1-S(t|xv,zv)
其中:
xv为预测样本的违约状态估计子模型信用特征变量;
zv为违约时间估计子模型信用特征变量;
Ov=(xv,zv)表示预测样本。
本发明还提供了一种P2P平台信用风险动态评价系统,包括:
数据采集单元,用于获取P2P平台的信用数据;
建模单元,所述建模单元包括:
定义信用状态单元,用于对P2P平台在观测时间内的删失状态和违约状态进行定义;
构建动态评价模型单元,用于构建P2P平台信用风险的动态评价模型、违约状态估计子模型以及违约时间估计子模型;
模型参数处理单元,所述模型参数处理单元包括:
构建生存观测数据单元,用于针对不同观测时间下所观测到的构建违约状态估计子模型以及违约时间估计子模型的信用特征变量,构建多重生存观测数据;
迭代参数估计单元,基于构建生存观测数据单元中所构建的多重生存观测数据,利用EM算法迭代估计模型参数;
概率计算单元,所述概率计算单元包括:
预测生存概率单元,用于预测P2P平台在不同时间下的生存概率;
违约概率计算单元,通过P2P平台在不同时间下的生存概率,计算出P2P平台在不同时间下的违约概率。
(三)有益效果
本发明提供了一种P2P平台信用风险动态评价方法和系统,该方法在获取P2P平台在观测时间内的删失状态数据和违约状态数据后定义P2P平台的信用状态,能有效处理P2P平台观测中出现的删失状态数据;在定义信用状态的基础上构建基于混合生存分析的动态评价模型,建立评价模型后针对不同观测时间,构建多重生存观测数据,然后采用EM算法迭代估计评价模型的参数,使模型能够有效学习到面板数据中的有效信息;针对选取预测样本,预测P2P平台在不同时间下的生存概率,从而计算出模型的动态违约概率,准确地预测P2P平台的风险趋势。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1本发明实施例一种P2P平台信用风险动态评价方法的流程图;
图2本发明实施例中多重生存观测构建的示意图;
图3本发明实施例一种P2P平台信用风险动态评价系统的框图;
其中,图2中的竖轴P2P平台1到4分别代表了四类P2P平台,横轴表示观测时间。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本申请实施例通过提供一种P2P平台信用风险动态评价方法和系统,解决了现有技术中对P2P平台信用风险评估不准确的技术问题,不仅能预测P2P平台是否会发生风险,而且能预测P2P平台何时会发生风险,提供P2P平台在不同时间下的动态违约概率。
本申请实施例中的技术方案为解决上述技术问题,其发明构思在于:
基于背景所述,现有的对P2P平台的信用风险评估,一般是沿用传统的信用评分模型或通过专家评估。传统的信用评分模型只能得到表现期或整个借款周期内的单一的违约概率,P2P平台的信用风险是变化的过程,仅仅依靠表现期或整个借款周期内的单一的违约概率无法准确评价P2P平台的风险趋势,专家评估则会受到人为因素干扰,导致评估结果出现误差。可见,这两类方法对P2P平台的信用风险评价均不准确。
为解决上述问题,本申请提出了一种基于混合生存分析的动态评价方法,用于预测P2P平台的动态信用风险。通过分别构建违约状态估计子模型和违约时间估计子模型来预测平台是否会发生违约以及何时会发生违约。方法设计了一种多重生存观测构建方法,使模型能够有效学习到平台信用特征面板数据中的有效信息,从而提升模型的预测效果,准确地预测平台的风险趋势,为P2P平台提供及时有效的内部风险控制,同时为投资人提供投资决策支持。
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
本发明实施例提供了一种P2P平台信用风险动态评价方法和系统,如图1所示,包括以下步骤S1—步骤S4:
S1、获取P2P平台的信用数据,信用数据包括观测时间以及观测时间内的信用特征变量、删失状态数据与违约状态数据;观测时间内的删失数据表示在观测时间内未能观测到违约。
S2、构建基于混合生存分析的动态评价模型,其中动态评价模型包括:P2P平台信用风险的动态评价模型、违约状态估计子模型以及违约时间估计子模型,具体步骤如下:
S201、定义P2P平台的信用状态,包括定义步骤S1中的观测时间、删失状态数据以及所述违约状态数据;
S202、通过S201的信用状态构建P2P平台信用风险的动态评价模型:P2P平台在指定观察时间下的生存概率等于P2P平台非违约的概率与P2P平台在指定观察时间之后发生违约的概率之和;
S203、在S202的基础上构建违约状态估计子模型:所述违约状态估计子模型为:针对P2P平台的违约概率,利用逻辑回归构建违约状态估计子模型;
S204、在S202的基础上构建违约时间估计子模型:所述违约时间估计子模型为:针对P2P平台的条件生存概率,利用Cox比例风险回归构建违约时间估计子模型。
S3、利用面板数据估计模型参数,具体步骤包括:
S301、步骤S2中的模型在不同观测时间下,可观测到构建步骤S2中的违约状态估计子模型以及违约时间估计子模型的信用特征变量,构成了P2P平台信用特征面板数据,针对不同观测时间,构建多重生存观测数据;
S302、基于步骤S301中所构建的多重生存观测数据,利用EM算法迭代估计步骤S2中动态评价模型的参数,这里的EM算法是指的是最大期望算法(ExpectationMaximization Algorithm)。
S4、预测P2P平台动态违约概率,通过步骤S1—步骤S3,得到训练好的动态评价模型,针对预测样本,预测P2P平台在不同时间下的生存概率;基于P2P平台在不同时间下的生存概率,可计算在不同时间下的违约概率。
本发明另一个实施例还提供了一种P2P平台信用风险动态评价系统,如图3所示,其包括:数据采集单元、建模单元、模型参数处理单元以及概率计算单元。
其中,数据采集单元用于获取P2P平台的信用数据。
建模单元包括:
定义信用状态单元,用于将P2P平台在观测时间内的删失状态和违约状态进行定义;
构建动态评价模型单元,用于构建P2P平台信用风险的动态评价模型、违约状态估计子模型以及违约时间估计子模型。
模型参数处理单元包括:
构建生存观测数据单元,用于针对不同观测时间下所观测到的违约状态估计子模型以及违约时间估计子模型的信用特征变量,构建多重生存观测数据;
迭代法估计单元,基于构建生存观测数据单元中所构建的多重生存观测数据,利用EM算法迭代估计模型参数。
概率计算单元包括:
预测生存概率单元,用于预测P2P平台在不同时间下的生存概率;
违约概率计算单元,通过P2P平台在不同时间下的生存概率,计算出P2P平台在不同时间下的违约概率。
本发明实施例提供了一种P2P平台信用风险动态评价方法和系统,通过分别构建违约状态估计子模型和违约时间估计子模型,该方法不仅能预测P2P平台是否会发生风险,而且能预测P2P平台何时会发生风险,预测P2P平台在不同时间下的动态违约概率。同时,设计了一种多重生存观测构建方法,使模型能够有效学习到面板数据中的有效信息,从而提升模型的预测效果,准确地预测平台的风险趋势。
下面对个步骤进行详细的描述:
S1、获取P2P平台的信用数据,在具体实施时,这里的步骤S1中的获取信用数据可以通过多种方式实现,下面对其中一种实施方式进行具体介绍:利用网络爬虫技术获取观测时间内的信用特征变量、删失状态数据和违约状态数据,删失状态数据表示在观测时间内未能观测到违约。
S2、构建基于混合生存分析的动态评价模型,包括:
S201、定义P2P平台的信用状态,包括定义步骤S1中的观测时间、删失状态数据以及违约状态数据:在具体实施时,这里的步骤S201中的定义P2P平台的信用状态有多种方式,下面对其中一种实施方式进行具体介绍:定义变量t为观测时间;定义随机变量δ为删失状态数据指示变量,δ=0表示观测数据删失,即在时间t内未能观测到违约,δ=1表示观测数据未删失,即在时间t内观测到违约;定义随机变量y为违约状态数据指示变量,y=0与表示P2P平台不会发生违约,y=1表示P2P平台会发生违约。因此,在某个观测时间t下P2P平台在时间存在着三种信用状态,如表1所示。
表1 P2P平台的信用状态分类
δ | y | 信用状态描述 |
0 | 0 | 删失,P2P平台未观测到违约,且不会违约 |
0 | 1 | 删失,P2P平台未观测到违约,但最终会违约 |
1 | 1 | 未删失,P2P平台已经观测到违约 |
S202、通过S201的信用状态构建P2P平台信用风险的动态评价模型:P2P平台在指定观察时间下的生存概率等于P2P平台非违约的概率与P2P平台在指定观察时间之后发生违约的概率之和。在具体实施时,这里的步骤S202中的构建P2P平台信用风险的动态评价模型可以通过多种方式实现,下面对其中一种实施方式进行具体介绍,具体表现为:
S(t)=1-p+p*S(t|y=1)
其中:
S(t)为目标函数,是P2P平台在时间t的生存概率;
p=P(y=1)是P2P平台的违约概率;
S(t|y=1)=P(T>t|y=1)是假定P2P平台会发生违约下的条件生存函数;
目标函数S(t)=1-p+p*S(t|y=1)表达的含义是:P2P平台在时间t下的生存概率等于P2P平台非违约的概率与P2P平台在时间t之后发生违约的概率之和。
S203、在S202的基础上构建违约状态估计子模型:针对P2P平台的违约概率,利用逻辑回归构建违约状态估计子模型。在具体实施时,可以通过多种方式实现,下面对其中一种实施方式进行具体介绍,具体表现为:
其中:
z=(1,z1,…,zn)为违约状态估计模型的信用特征变量;
β=(β0,β1,…,βn)为信用特征变量的系数。
S204、在S202的基础上构建违约时间估计子模型:针对P2P平台的条件生存概率,利用Cox比例风险回归构建违约时间估计子模型。
具体表现为:
其中:
S0(t|y=1)为基准生存函数,是函数的非参数部分;
向量x=(x1,x2,…,xm)为违约时间估计模型的信用特征变量;
向量α=(α1,α2,…,αm)为信用特征变量的系数。
S3、基于面板数据的模型参数估计,包括:
S301、步骤S2中的模型在不同观测时间下,可观测到违约状态估计子模型以及违约时间估计子模型的信用特征变量,构成了P2P平台信用特征面板数据,针对不同观测时间,构建多重生存观测数据。在具体实施时,可以通过多种方式实现,下面对其中一种实施方式进行具体介绍:
在不同观测时间下,可观测到P2P平台违约状态估计子模型以及违约时间估计子模型的信用特征变量,构成了P2P平台信用特征面板数据。针对不同观测时间,构建多重生存观测数据。例如,如图2所示,P2P平台1到4分别代表了四类P2P平台,即P2P平台1在观测到违约,P2P平台2在观测到违约,P2P平台3在观测到违约,P2P平台4在观测周期内未观测到违约,因此在处删失。针对P2P平台1,在时刻开始观测,在时刻观测到P2P平台发生违约,则观测生存数据的信用状态为。针对P2P平台2,可以得到两条生存观测数据,第多重生存数据在时刻开始观测,在时刻观测到P2P平台发生违约,观测生存数据的信用状态为,第二条生存数据在时刻开始观测,在时刻观测到违约,观测生存数据的信用状态为。类似地,针对P2P平台3和P2P平台4可以分别得到三条生存观测数据,其信用状态如表2所示。
表2多重生存观测的信用状态(t,δ)
P2P平台ID | 观测1 | 观测2 | 观测3 |
1 | (1,1) | ||
2 | (2,1) | (1,1) | |
3 | (3,1) | (2,1) | (1,1) |
4 | (3,0) | (2,0) | (1,0) |
S302、基于步骤S301中所构建的多重生存观测数据,利用EM算法迭代估计步骤S2中动态评价模型的参数。在具体实施时,可以通过多种方式实现,下面对其中一种实施方式进行具体介绍,具体表现为:
步骤S3021:
利用删失变量作为初始违约概率期望值E(yi)(0)=δi;
估计违约概率估计子模型的初始参数β(0)以及违约时间估计子模型的初始参数α(0),S0 (0);
设定模型迭代的最大次数C,模型收敛函数g(·)=+∞以及收敛阈值G,迭代次数指示变量k=1;
若迭代次数k<C,且收敛函数g(·)>G,则进行以下迭代;
步骤S3022:
利用违约概率期望值E(yi)(k-1)计算违约概率估计子模型的极大似然估计函数:
其中:
LI为违约概率估计子模型的极大似然函数;
p(zi)为违约状态估计子模型结果。
得到该轮迭代下违约概率估计子模型的参数估计值:
β(k)=(β1,β2,…,βn)(k)
得到该轮迭代下违约概率估计子模型的参数估计值
β(k)=(β1,β2,…,βn)(k)
步骤S3023:
利用违约概率期望值E(yi)(k-1)计算违约时间估计子模型的极大似然估计函数:
其中:
LL为违约时间估计子模型的极大似然函数;
h(ti|y=1,xi)为条件生存函数S(ti|y=1,xi)所对应的风险函数;
S(ti|y=1,xi)是假定平台会发生违约下的条件生存函数。
得到该轮迭代下违约概率估计子模型的参数估计值:
α(k)=(α1,α2,...,αm)(k)
步骤S3024:
计算违约概率yi的期望值:
其中:
E(k)(yi|ti,δi,zi,xi)为k轮迭代下违约概率yi的期望值。
步骤S3025:
计算迭代收敛函数:
g(·)=||β(k)-β(k-1)||2 2+||α(k)-α(k-1)||2 2+||S0 (k)-S0 (k-1)||2 2
其中:
g(·)为迭代收敛函数;
β(k)为k轮迭代下违约状态估计子模型的参数值;
α(k),S0 (k)为k轮迭代下违约时间估计子模型的参数值。
步骤S3026:
更新迭代次数:
k=k+1
过程S3022至S3026一直循环迭代,直至迭代达到最大次数(k=C)或迭代收敛函数小于设定的阈值(g(k)≤G);
得到动态评价函数的参数估计值:
β(*)=(β0,β1,...,βn)(*)
α(*)=(α1,α2,,αm)(*)
S0 (*)(t|y=1)
β(*),α(*),S0 (*)(t|y=1)分别为迭代结束时违约状态估计子模型和违约时间估计子模型的参数值。
S4、预测P2P平台动态违约概率,在具体实施时,可以通过多种方式实现,下面对其中一种实施方式进行具体介绍,具体表示为:
其表示:针对预测样本,预测P2P平台在不同时间下的生存概率;
违约概率的表示方法为:
P(t|xv,zv)=1-S(t|xv,zv)
其中:
xv为预测样本的违约状态估计子模型信用特征变量;
zv为违约时间估计子模型信用特征变量;
Ov=(xv,zv)表示预测样本。
将上述实施例用于对现有的P2P平台进行评估,具体过程如下:
选取了396家P2P平台第三季度(2018年7月至2018年9月)的面板数据,其中正常运营P2P平台311家,违约P2P平台85家,其中7月违约的有52家,8月违约的有15家,9月违约的有18家。变量包括成交量,参考收益率,投资人数,借款人数,新投资人数,老投资人数,人均投资额,人均借款额,投资人总额,平均借款期限。由于变量表现出高度的偏态分布,针对变量进行了log(x+1)转换。
为了获得预测性能的无偏估计,选择了留一法交叉验证来估计模型的预测性能。针对模型的预测性能,选取了AUC,KS和H-measure三个评价指标对模型进行综合评价,其中AUC计算了ROC曲线下面的面积,其反应了模型对违约和非违约样本的综合预测能力,并考虑了样本不均衡的影响。KS计算了模型预测违约样本和非违约样本累计分布的最大差值,反应了模型对违约样本和非违约样本的分隔程度。H-measure克服了AUC指标在评价不同模型的预测性能中应用了不同损失分布的缺陷,利用beta分布明确了模型的误分类损失,从而获得了标准化的预测性能评估。
选择了三类模型与本发明中的模型(MSA)进行对比:
(1)基于截面数据的动态评价模型(MCM):利用截面数据,训练了本发明中的动态评价模型,其中违约状态估计子模型和违约时间估计子模型应用了逻辑回归和Cox比例风险回归。
(2)仅考虑何时违约的评价模型(COX):利用Cox比例风险回归模型构建了仅考虑何时违约的评价模型。
(3)仅考虑是否违约的评价模型(LR):利用利用逻辑回归构建了仅考虑是否违约的评价模型。
首先验证了多重生存观测方法的有效性,比较了MSA模型和MCM模型在以下三个场景的预测性能:
(1)P2P平台在前一个月发生违约的概率P(T≤1);
(2)P2P平台在前两个月发生违约的概率P(T≤2);
(3)P2P平台在前三个月发生违约的概率P(T≤3)。
表3动态违约概率估计
模型 | 时间 | 评分函数 | AUC | KS | H-measure |
MSA | 1 | 1-S(1) | 0.814 | 0.479 | 0.358 |
MCM | 0.739 | 0.358 | 0.251 | ||
MSA | 2 | 1-S(2) | 0.817 | 0.525 | 0.366 |
MCM | 0.757 | 0.421 | 0.266 | ||
MSA | 3 | 1-S(3) | 0.828 | 0.553 | 0.392 |
MCM | 0.776 | 0.492 | 0.298 |
模型违约概率的预测性能结果如表3所示。在三个时间下,动态评价模型的评分函数分别为1-S(1)、1-S(2)以及1-S(3)。留一法交叉验证的结果显示,在所有的时间下基于多重生存观测估计的MSA模型的预测效果均显著优于基于截面观测估计的MCM模型,且在AUC,KS以及H-measure三种预测性能评价指标下MSA的模型的效果均优于MCM模型,体现出MSA模型预测效果的稳健性。结果表明,相比于截面数据,P2P平台的面板数据中包含了更多违约风险判别的有用信息,本发明所提出的多重生存观测方法可以有效的利用P2P平台面板数据中的有用信息,从而提升P2P平台风险动态评价模型的违约预测性能。
验证了多重生存观测方法的有效性后,进而比较了MSA、MCM、COX以及LR模型在两个动态场景下的违约预测性能:
(1)P2P平台在第一个月未发生违约下,在第二个月发生违约的概率;
(2)P2P平台在前两个月未发生违约下,在第三个月发生违约的概率。由于逻辑回归模型无法预测P2P平台随时间的动态违约概率,为了比较逻辑回归与其他模型在动态违约预测中的性能,在不同时间分别构建了逻辑回归模型,并综合多个逻辑回归的预测结果来估计违约概率。
表4违约风险预测性能
模型 | 时间 | 评分函数 | AUC | KS | H-measure |
MSA | 2|1 | 1-S(2)/S(1) | 0.789 | 0.536 | 0.344 |
MCM | 1-S(2)/S(1) | 0.760 | 0.554 | 0.315 | |
COX | 1-S(2)/S(1) | 0.798 | 0.566 | 0.330 | |
LR | (p(2)-p(1))/(1-p(1)) | 0.678 | 0.496 | 0.253 | |
MSA | 3|2 | 1-S(3)/S(2) | 0.832 | 0.595 | 0.421 |
MCM | 1-S(3)/S(2) | 0.796 | 0.561 | 0.418 | |
COX | 1-S(3)/S(2) | 0.774 | 0.589 | 0.379 | |
LR | (p(3)-p(2))/(1-p(2)) | 0.821 | 0.588 | 0.413 |
模型的违约预测结果如表4所示。由留一法交叉验证的结果显示:
评价场景(2)下,本发明提出的MSA模型的动态违约预测性能显著优于其他三个模型,且在AUC,KS以及H-measure三种评价指标下均取得了最好的预测性能。
在评价场景(1)下,基于H-measure的标准化评价指标下,MAS模型取得了最高的预测性能,基于AUC和KS的评价指标下,COX模型取得了最高的预测性能。尽管在该场景下取得了较好的预测性能,但COX模型在场景(2)下的预测性能却低于其他模型,处于最低水平。同时,LR模型也表现出相似的情形,LR模型在场景(2)下取得较好的预测性能,仅次于MSA模型,但在场景(1)下却表现出最差的预测性能。由此可见,仅仅考虑是否违约(LR模型)或何时违约(COX模型)单方面因素,评价模型的动态预测性可能会出现波动性。巴塞尔新资本协议中明确指出银行或金融机构在开发内部评级系统时需考虑评价模型的稳定性,高质量的信用评价模型应能在不同时间下均取得良好的区分能力。基于留一法交叉验证的结果,本发明所提出的MSA模型在不同的时间下均取得了较好且稳定的违约判别性能,能够更为准确地预测P2P平台的风险趋势,因此可以更好地为P2P借贷投资人提供投资决策辅助,同时为P2P平台提供及时有效的内部风险控制。
综上所述,与现有技术相比,具备以下有益效果:
本发明提供了一种P2P平台信用风险动态评价方法和系统,该方法在获取P2P平台在观测时间内的删失状态数据和违约状态数据后定义平台的信用状态,能有效处理P2P平台观测中出现的删失数据;在定义信用状态的基础上构建基于混合生存分析的动态评价模型,建立评价模型后针对不同观测时间,构建多重生存观测数据,然后采用EM算法迭代估计评价模型的参数,使模型能够有效学习到面板数据中的有效信息;针对选取预测样本,预测P2P平台在不同时间下的生存概率,从而计算出模型的动态违约概率,准确地预测P2P平台的风险趋势。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (8)
1.一种P2P平台信用风险动态评价方法,其特征在于,该方法包括以下步骤:
S1、获取P2P平台的信用数据,所述信用数据包括:观测时间、观测时间内的信用特征变量、观测时间内的删失状态数据以及观测时间内的违约状态数据;所述观测时间内的删失状态数据表示在观测时间内未能观测到违约;
S2、构建基于混合生存分析的动态评价模型,包括:
S201、定义P2P平台的信用状态,包括定义所述步骤S1中的所述观测时间、所述删失状态数据以及所述违约状态数据;
S202、基于S201的信用状态构建P2P平台信用风险的动态评价模型;
S203、在S202的基础上构建违约状态估计子模型;
S204、在S202的基础上构建违约时间估计子模型;
S3、利用面板数据估计模型参数,包括:
S301、步骤S2中的模型在不同观测时间下,可观测到构建违约状态估计子模型以及构建违约时间估计子模型的信用特征变量,构成了P2P平台信用特征面板数据,针对不同观测时间,构建多重生存观测数据;
S302、基于步骤S301中所构建的多重生存观测数据,利用EM算法迭代估计步骤S2中动态评价模型的参数;
S4、预测P2P平台动态违约概率,通过步骤S1—步骤S3,得到训练好的动态评价模型,针对预测样本,预测P2P平台在不同时间下的生存概率;基于P2P平台在不同时间下的生存概率,可计算在不同时间下的违约概率;
其中,在步骤S202中:所述动态评价模型为:P2P平台在指定观察时间下的生存概率等于P2P平台非违约的概率与P2P平台在指定观察时间之后发生违约的概率之和;
在步骤S203中:所述违约状态估计子模型为:针对P2P平台的违约概率,利用逻辑回归构建违约状态估计子模型;
在步骤S204中:所述违约时间估计子模型为:针对P2P平台的条件生存概率,利用Cox比例风险回归构建违约时间估计子模型。
2.根据权利要求1所述的P2P平台信用风险动态评价方法,其特征在于,在步骤S201中,所述定义P2P平台的信用状态:
定义变量t为观测时间;
定义随机变量δ为删失状态数据指示变量,δ=0表示观测数据删失,即在时间t内未能观测到违约,δ=1表示观测数据未删失,即在t时间内观测到违约;
定义随机变量y为违约状态数据指示变量,y=0与表示P2P平台不会发生违约,y=1表示P2P平台会发生违约。
3.根据权利要求2所述的P2P平台信用风险动态评价方法,其特征在于,在步骤S202中,所述P2P平台信用风险的动态评价模型为:
S(t)=1-p+p*S(t|y=1)
其中:
S(t)为目标函数,是P2P平台在时间t的生存概率;
p=P(y=1)是P2P平台的违约概率;
S(t|y=1)=P(T>t|y=1)是假定P2P平台会发生违约下的条件生存函数;
目标函数S(t)=1-p+p*S(t|y=1)表达的含义是:P2P平台在时间t下的生存概率等于P2P平台非违约的概率与P2P平台在时间t之后发生违约的概率之和。
4.根据权利要求3所述的P2P平台信用风险动态评价方法,其特征在于,在步骤S203中,所述违约状态估计子模型为:
其中:
z=(1,z1,…,zn)为违约状态估计子模型的信用特征变量;
β=(β0,β1,…,βn)为信用特征变量的系数。
5.根据权利要求4所述的P2P平台信用风险动态评价方法,其特征在于,在步骤S204中,所述违约时间估计子模型为:
其中:
S0(t|y=1)为基准生存函数,是函数的非参数部分;
向量x=(x1,x2,…,xm)为违约时间估计子模型的信用特征变量;
向量α=(α1,α2,…,αm)为信用特征变量的系数。
6.根据权利要求5所述的P2P平台信用风险动态评价方法,其特征在于,在所述步骤S302中,包括:
步骤S3021:
利用删失变量作为初始违约概率期望值E(yi)(0)=δi;
估计违约概率估计子模型的初始参数β(0)以及违约时间估计子模型的初始参数α(0),S0 (0);
设定模型迭代的最大次数C,模型收敛函数g(·)=+∞以及收敛阈值G,迭代次数指示变量k=1;
若迭代次数k<C,且收敛函数g(·)>G,则进行以下迭代;
步骤S3022:
利用违约概率期望值E(yi)(k-1)计算违约概率估计子模型的极大似然估计函数:
其中:
LI为违约概率估计子模型的极大似然函数;
p(zi)为违约状态估计子模型结果;
得到该轮迭代下违约概率估计子模型的参数估计值:
β(k)=(β1,β2,…,βn)(k)
步骤S3023:
利用违约概率期望值E(yi)(k-1)计算违约时间估计子模型的极大似然估计函数:
其中:
LL为违约时间估计子模型的极大似然函数;
h(ti|y=1,xi)为条件生存函数S(ti|y=1,xi)所对应的风险函数;
S(ti|y=1,xi)是假定平台会发生违约下的条件生存函数;
得到该轮迭代下违约概率估计子模型的参数估计值:
α(k)=(α1,α2,…,αm)(k)
步骤S3024:
计算违约概率yi的期望值:
其中:
E(k)(yi|ti,δi,zi,xi)为k轮迭代下违约概率yi的期望值;
步骤S3025:
计算迭代收敛函数:
g(·)=||β(k)-β(k-1)||2 2+||α(k)-α(k-1)||2 2+||S0 (k)-S0 (k-1)||2 2
其中:
g(·)为迭代收敛函数;
β(k)为k轮迭代下违约状态估计子模型的参数值;
α(k),S0 (k)为k轮迭代下违约时间估计子模型的参数值;
步骤S3026:
更新迭代次数:
k=k+1
步骤S3022至S3026一直循环迭代,直至迭代达到最大次数(k=C)或迭代收敛函数小于设定的阈值(g(k)≤G);
得到动态评价函数的参数估计值:
β(*)=(β0,β1,…,βn)(*)
α(*)=(α1,α2,…,αm)(*)
S0 (*)(t|y=1)
其中:
β(*),α(*),S0 (*)(t|y=1)分别为迭代结束时违约状态估计子模型和违约时间估计子模型的参数值。
7.根据权利要求6所述的P2P平台信用风险动态评价方法,其特征在于,在步骤S4中,
预测P2P平台在不同时间下的生存概率表示为:
其表示:针对预测样本,预测P2P平台在不同时间下的生存概率;
违约概率的表示方法为:
P(t|xv,zv)=1-S(t|xv,zv)
其中:
xv为预测样本的违约状态估计子模型信用特征变量;
zv为违约时间估计子模型信用特征变量;
Ov=(xv,zv)表示预测样本。
8.一种P2P平台信用风险动态评价系统,其特征在于,包括:
数据采集单元,用于获取P2P平台的信用数据;
建模单元,所述建模单元包括:
定义信用状态单元,用于对P2P平台在观测时间内的删失状态和违约状态进行定义;
构建动态评价模型单元,用于构建P2P平台信用风险的动态评价模型、违约状态估计子模型以及违约时间估计子模型;
模型参数处理单元,所述模型参数处理单元包括:
构建生存观测数据单元,用于针对不同观测时间下所观测到的构建违约状态估计子模型以及违约时间估计子模型的信用特征变量,构建多重生存观测数据;
迭代参数估计单元,基于构建生存观测数据单元中所构建的多重生存观测数据,利用EM算法迭代估计模型参数;
概率计算单元,所述概率计算单元包括:
预测生存概率单元,用于预测P2P平台在不同时间下的生存概率;
违约概率计算单元,通过P2P平台在不同时间下的生存概率,计算出P2P平台在不同时间下的违约概率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910127901.XA CN109961362A (zh) | 2019-02-19 | 2019-02-19 | P2p平台信用风险动态评价方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910127901.XA CN109961362A (zh) | 2019-02-19 | 2019-02-19 | P2p平台信用风险动态评价方法和系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109961362A true CN109961362A (zh) | 2019-07-02 |
Family
ID=67023702
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910127901.XA Pending CN109961362A (zh) | 2019-02-19 | 2019-02-19 | P2p平台信用风险动态评价方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109961362A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110400013A (zh) * | 2019-07-22 | 2019-11-01 | 西北工业大学 | 一种基于多任务学习机制的移动应用消亡预测方法 |
CN110689427A (zh) * | 2019-10-12 | 2020-01-14 | 杭州绿度信息技术有限公司 | 一种基于生存分析的消费分期违约概率模型 |
CN110827143A (zh) * | 2019-11-08 | 2020-02-21 | 支付宝(杭州)信息技术有限公司 | 训练信用评分模型的方法、装置及设备 |
CN110929223A (zh) * | 2019-11-14 | 2020-03-27 | 中国科学技术大学 | 债权转让匹配方法、系统、设备及介质 |
CN112750022A (zh) * | 2020-12-16 | 2021-05-04 | 苏宁消费金融有限公司 | 一种消费金融的信用压缩方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004057503A2 (en) * | 2002-12-20 | 2004-07-08 | Accenture Global Services Gmbh | Quantification of operational risks |
CN102075356A (zh) * | 2010-12-31 | 2011-05-25 | 深圳市永达电子股份有限公司 | 一种网络风险评估方法与系统 |
CN106779755A (zh) * | 2016-12-31 | 2017-05-31 | 湖南文沥征信数据服务有限公司 | 一种网络电商借贷风险评估方法及模型 |
CN106920110A (zh) * | 2017-02-23 | 2017-07-04 | 国网安徽省电力公司合肥供电公司 | 一种电力客户信用和欠费风险的评价方法 |
CN107194803A (zh) * | 2017-05-19 | 2017-09-22 | 南京工业大学 | 一种p2p网贷借款人信用风险评估的装置 |
CN107657529A (zh) * | 2017-10-26 | 2018-02-02 | 陶冶 | Ifrs9的信用卡预期减值方法 |
CN108492001A (zh) * | 2018-02-13 | 2018-09-04 | 天津大学 | 一种用于担保贷款网络风险管理的方法 |
US20180308160A1 (en) * | 2016-02-01 | 2018-10-25 | Tencent Technology (Shenzhen) Company Limited | Risk assessment method and system |
-
2019
- 2019-02-19 CN CN201910127901.XA patent/CN109961362A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004057503A2 (en) * | 2002-12-20 | 2004-07-08 | Accenture Global Services Gmbh | Quantification of operational risks |
CN102075356A (zh) * | 2010-12-31 | 2011-05-25 | 深圳市永达电子股份有限公司 | 一种网络风险评估方法与系统 |
US20180308160A1 (en) * | 2016-02-01 | 2018-10-25 | Tencent Technology (Shenzhen) Company Limited | Risk assessment method and system |
CN106779755A (zh) * | 2016-12-31 | 2017-05-31 | 湖南文沥征信数据服务有限公司 | 一种网络电商借贷风险评估方法及模型 |
CN106920110A (zh) * | 2017-02-23 | 2017-07-04 | 国网安徽省电力公司合肥供电公司 | 一种电力客户信用和欠费风险的评价方法 |
CN107194803A (zh) * | 2017-05-19 | 2017-09-22 | 南京工业大学 | 一种p2p网贷借款人信用风险评估的装置 |
CN107657529A (zh) * | 2017-10-26 | 2018-02-02 | 陶冶 | Ifrs9的信用卡预期减值方法 |
CN108492001A (zh) * | 2018-02-13 | 2018-09-04 | 天津大学 | 一种用于担保贷款网络风险管理的方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110400013A (zh) * | 2019-07-22 | 2019-11-01 | 西北工业大学 | 一种基于多任务学习机制的移动应用消亡预测方法 |
CN110689427A (zh) * | 2019-10-12 | 2020-01-14 | 杭州绿度信息技术有限公司 | 一种基于生存分析的消费分期违约概率模型 |
CN110827143A (zh) * | 2019-11-08 | 2020-02-21 | 支付宝(杭州)信息技术有限公司 | 训练信用评分模型的方法、装置及设备 |
CN110929223A (zh) * | 2019-11-14 | 2020-03-27 | 中国科学技术大学 | 债权转让匹配方法、系统、设备及介质 |
CN110929223B (zh) * | 2019-11-14 | 2022-05-17 | 中国科学技术大学 | 债权转让匹配方法、系统、设备及介质 |
CN112750022A (zh) * | 2020-12-16 | 2021-05-04 | 苏宁消费金融有限公司 | 一种消费金融的信用压缩方法 |
CN112750022B (zh) * | 2020-12-16 | 2022-08-02 | 苏宁消费金融有限公司 | 一种消费金融的信用压缩方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109961362A (zh) | P2p平台信用风险动态评价方法和系统 | |
Pan et al. | Research on the competitive effect of local government's environmental expenditure in China | |
CN111191918A (zh) | 一种智能电网通信网的业务路由规划方法及装置 | |
CN112632765B (zh) | 一种加权和法与sem法结合的作战能力评估方法 | |
CN109118120B (zh) | 考虑水库调度方案可持续性利用的多目标决策方法 | |
Heinrich et al. | Diversity of firm sizes, complexity, and industry structure in the Chinese economy | |
CN111695754A (zh) | 一种电力物联网信息安全风险评估方法和装置 | |
CN111967954A (zh) | 一种资源归还增比确定方法、装置和电子设备 | |
CN112508679A (zh) | 一种小微企业贷款风险评估方法、装置及存储介质 | |
Cheng et al. | A decision-making framework for load rating planning of aging bridges using deep reinforcement learning | |
Harahap | Prediction of financial distress in foreign exchange banking firms using risk analysis, good corporate governance, earnings, and capital | |
Roveto et al. | Co-optimization of VaR and CVaR for data-driven stochastic demand response auction | |
Zhijun | RBF neural networks optimization algorithm and application on tax forecasting | |
CN111932106B (zh) | 一种有效实用的云制造任务与服务资源匹配方法 | |
CN116911994B (zh) | 对外贸易风险预警系统 | |
CN111382909A (zh) | 基于生存分析模型扩展坏样本的拒绝推断方法及相关设备 | |
CN110866694A (zh) | 一种电网建设项目财务评价系统及其方法 | |
CN115759711A (zh) | 面向负荷管理的需求响应执行效果综合评价方法和系统 | |
CN109978138A (zh) | 基于深度强化学习的结构可靠度抽样方法 | |
CN115629576A (zh) | 非侵入式柔性负荷聚合特性辨识与优化方法、装置及设备 | |
CN115829334A (zh) | 一种电网业务的风险评估方法及系统 | |
Finger et al. | Friendship between banks: An application of an actor-oriented model of network formation on interbank credit relations | |
Majdzadeh et al. | Determinants of private investment in Iran based on Bayesian Model Averaging | |
CN108053266A (zh) | 一种专利价值预估方法以及装置 | |
Ma et al. | Simulation analysis of the evolution of sustainable operation of transport infrastructure projects under government regulation based on prospect theory and BP neural network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |