随着我国人口老龄化的加剧和住房率的不断提高,以房养老作为一种新型养老金融产品,受到了社会各界的广泛关注。是否获得足够的养老金决定了老年人是否拥有更好的生活,这凸显了控制潜在风险的必要性。因此,建立有效的以房养老风险评估体系是养老保险行业不可或缺的。以住房反向抵押贷款为例进行模型分析,与租房养老金和租房养老金相比,住房反向抵押贷款的风险介于普通风险和略高风险之间,应采取措施控制和最小化这些风险。 With the aggravation of population aging and the continuous improvement of housing rate in China, house-for-pension as a new type of pension financial products, has been attracted widespread attention from all walks of life. However, better life of the elderly depends muchon sufficient pension, which emphasize the necessity of evaluating and preventing the potential risks. Hence, it is necessary to establish an efficient risk evaluation system of House-for-pension in China. Taking Reverse Mortgage, House-renting Pension and House-selling Pension as examples for model analysis, it was concluded that the risk of Reverse Mortgage is 3.7912, the risk of House-renting Pension is 3.4164, the risk of House-selling Pension is 3.7357. All models of House-for-pension products in China are slightly high, thus measures should be taken to control and minimize the potential risks.
随着我国人口老龄化的加剧和住房率的不断提高,以房养老作为一种新型养老金融产品,受到了社会各界的广泛关注。是否获得足够的养老金决定了老年人是否拥有更好的生活,这凸显了控制潜在风险的必要性。因此,建立有效的以房养老风险评估体系是养老保险行业不可或缺的。以住房反向抵押贷款为例进行模型分析,与租房养老金和租房养老金相比,住房反向抵押贷款的风险介于普通风险和略高风险之间,应采取措施控制和最小化这些风险。
AHP模糊综合评价,以房养老,定量风险评估,风险评估
Lin Qiu
School of Management, Shanghai University of Engineering Science, Shanghai
Received: May 18th, 2023; accepted: Jun. 20th, 2023; published: Jun. 29th, 2023
With the aggravation of population aging and the continuous improvement of housing rate in China, house-for-pension as a new type of pension financial products, has been attracted widespread attention from all walks of life. However, better life of the elderly depends muchon sufficient pension, which emphasize the necessity of evaluating and preventing the potential risks. Hence, it is necessary to establish an efficient risk evaluation system of House-for-pension in China. Taking Reverse Mortgage, House-renting Pension and House-selling Pension as examples for model analysis, it was concluded that the risk of Reverse Mortgage is 3.7912, the risk of House-renting Pension is 3.4164, the risk of House-selling Pension is 3.7357. All models of House-for-pension products in China are slightly high, thus measures should be taken to control and minimize the potential risks.
Keywords:AHP Fuzzy Comprehensive Evaluation, House-for-Pension, Quantitative Risk Assessment, Risk Evaluation
Copyright © 2023 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
1999年,占世界老年人口五分之一的中国进入老龄化社会,成为世界上老龄人口最多的国家。截至2021年底,中国有2.6736亿60岁以上老年人,占总人口的18.9% [
风险评估可以定义为组织信息以支持管理过程中做出的风险决策的系统过程 [
图1. 典型风险评估流程概述
作为风险评估的开始部分,风险识别是系统地使用信息来识别和记录与风险问题或问题描述相关的危险 [
风险分析是对与设想危险相关的风险的估计。这是一个定性或定量的过程,将事故发生的可能性和由此产生的危害的严重程度联系起来 [
从不同主体(借款人和贷款机构)的角度出发,许多学者开始了以房养老风险评估的研究,并成功地取得了许多研究成果。房价、寿命和利率风险是养老贷款机构面临的三个主要问题。最近的文献研究了影响养老贷款机构住房的各种因素。例如,在一个随机多周期模型中,Shao等人(2015)研究了房价和寿命风险对住房养老保险定价和风险状况的综合影响 [
几项研究考察了HECM定价模型所解释的与住房换养老金产品相关的风险。这些风险包括利率、房价和寿命风险。Kogure等人(2014)还提出了一个考虑利率、房价和寿命风险的房屋养老价值贝叶斯多变量模型 [
通过对上述反向抵押的相关研究可以看出,随着人口老龄化的不断加剧和人均可支配收入的增加,以房养老在中国已经成为可能。以房养老可以有效增加老年人的可支配收入,并改善多层次老年保险制度。然而,现有关于以房养老风险的文献大多是从金融机构的角度出发,而不是从老年人、金融机构和抵押物的多维度角度系统地评估以房养老的风险。因此,本文从借款人、贷款人、抵押物和社会的角度对以房养老进行了全面、系统的风险评估。
本研究基于AHP-fuzzy综合评价模型进行
AHP-fuzzy综合评估的两个主要原则是:1) 质量风险评估应基于科学知识,并最终与患者的保护相联系;2) 质量风险管理过程的工作水平、形式和文件应与风险水平相称。
作者主要关注中国住房养老所带来的风险,该模型将住房养老分为4个部分:借款人风险、贷款人风险、抵押物风险和社会风险。表1显示了中国以房养老的其他具体量化风险评估。
风险因素 | 一级指标 | 二级指标 |
---|---|---|
以房养老风险评估的构成要素(A) | 借款人风险 B1 | 传统伦理风险Z1 |
通货膨胀风险Z2 | ||
成本风险Z3 | ||
老年财务风险Z4 | ||
贷款人风险 B2 | 长寿风险Z5 | |
道德风险Z6 | ||
运营风险Z7 | ||
信任风险Z8 | ||
抵押物风险 B3 | 房屋价格波动风险Z9 | |
房屋产权风险Z10 | ||
社会风险 B4 | 利率风险Z11 | |
政策风险Z12 | ||
或然风险Z13 |
表1. 以房养老的风险因素
AHP-fuzzy综合评价模型可分为以下步骤:
1) 指标集和评估集的定义
根据模糊集理论,必须综合考虑多个相关因素,才能做出恰当、不矛盾和逻辑一致的判断。
索引集 U = ( U 1 , U 2 , ⋯ , U n ) , i ∈ ( 1 , 2 , 3 , ⋯ , S ) ;评估集 V = ( 1 , 2 , 3 , 4 , 5 ) 。
指标集U包含公认的风险要素,评估集V为专家提供评估基础。高风险可替换为5;略高风险4;常见风险3;略低风险2;低风险1。
2) 通过AHP确定权重
为每个较低级别构造一组成对比较矩阵。进行成对比较,使行 i ( i = 1 , 2 , 3 , 4 , ⋯ , n ) 中的属性相对于列表示的每个属性n进行排序。成对比较是根据哪个元素支配另一个元素(即基于元素的相对重要性)进行的。然后将这些判断表示为表2和表3中所示的整数值1~9。
数字比率 | 内涵 |
---|---|
1 | 因素i和因素j相等 |
3 | 因素i比因素j稍微重要一些 |
5 | 因素i显然比因素j更重要 |
7 | 因素i比因素j重要得多 |
9 | 因素i比因素j极其重要 |
2,4,6,8 | 中间值 |
倒数 | a i :因素i和因素j相比, a i j :因素j和因素i相比 a j i = 1 / a i j |
表2. 判断矩阵中要素的值
A | B1 | B2 | …… | Bn |
---|---|---|---|---|
B1 | 1 | a12 | …… | a1n |
B2 | a21 | 1 | …… | a2n |
…… | …… | …… | …… | …… |
Bn | an1 | an2 | …… | 1 |
表3. 判断矩阵
如果 a i j = 1 / a j i 和 a i j = B i / B j , a i j > 0 ,矩阵A是一个正互易矩阵 A a i j * a j k = a i k , i , j , k = 1 , 2 , ⋯ , n ,并且如果在矩阵中,则矩阵A是完全一致的矩阵。对于完全一致的矩阵A,存在唯一的非零特征值,并使用相应的特征向量来计算权重向量W。最后,为了避免人为误差和不同因素的矛盾,进行一致性测试,直到满足满意的条件。
Satty建议使用一致性指数(CI)和一致性比率(CR)来检查与比较矩阵相关的一致性。
CI = λ − n n − 1
如果CI = 0,则矩阵A具有完全一致性。如果CI ≠ 0,但CI接近0,则矩阵A具有令人满意的一致性。为了有效测量CI,Saaty引入RI进行计算(表4)。他随机构建了500个比较矩阵 A 1 , A 2 , ⋯ , A 500 ,得到了结果 CI 1 , CI 2 , ⋯ , CI 500
n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 |
---|---|---|---|---|---|---|---|---|---|---|---|
RI | 0 | 0 | 0.58 | 0.90 | 1.12 | 1.24 | 1.32 | 1.41 | 1.45 | 1.49 | 1.51 |
表4. 随机一致性指标RI
RI = CI 1 + CI 2 + ⋯ + CI 500 500 = ( λ 1 + λ 2 + ⋯ + λ 500 ) / 500 − n n − 1
最后我们得到了等式: CR = CI RI ,当 CR < 0.1 ,不一致程度被接受时。换句话说,当一致性测试通过时,归一化特征向量可以被视为权重向量。
3) 模糊评价矩阵的确定 R i j = k / n :
R = [ r 11 r 12 ⋯ r 1 n r 21 r 22 ⋯ r 2 n ⋮ ⋮ ⋱ ⋮ r m 1 r m 2 ⋯ r m n ] 。
4) 结果
AHP-fuzzy综合评价的结果使用以下公式计算:
E = W * C 。
其中W是对应于风险要素(1类)的权重矩阵;C是与风险要素(1级)相对应的评估得分。C由下式表示:
C = B * V 。
其中V是百分变换矩阵;
V = [ 1 2 3 4 5 ] 。
B是对应于每个风险要素(2级)的准确得分矩阵。B由下式表示:
B = A * R 。
其中A是对应于风险要素的权重矩阵;R是灰色决策矩阵。
经过一系列计算,我们得到了AHP模糊综合评估的结果,然后按照风险等级的顺序以百分的方式转换结果,并确保辅料的风险等级(表5)。
分数 | 5 | 4 | 3 | 2 | 1 |
---|---|---|---|---|---|
风险等级 | 风险很高 | 风险较高 | 风险一般 | 风险较低 | 风险很低 |
表5. 风险等级
本文以我国三种以房养老为例,说明AHP-fuzzy综合评价模型的应用。
住房反向抵押贷款是指老年人将自有房屋抵押给特定机构,机构根据房屋估价、预期利率和老年人预期寿命等,一次性或者分期给付给老年人养老金,在老年人去世之后回收房屋予以变现 [
租房养老是指不转移老人自持物业的产权,而是将房产抵押给金融机构,由金融机构将房产出租,获取的收益按比例归属老人一部分,待老人百年之后其子女或继承人可以接力还贷款,也可以选择拍卖变现 [
售房养老是指老年人在壮年时期通过储蓄存款和按揭贷款的方式购买属于自己的房屋,在退休之后将该房屋的产权出售给特定机构并保留房屋使用权,特定机构分期或是一次性给付给老年人养老金,在老年人去世之后回收房屋予以变现 [
1) 指标集和评估集的定义
如上所述,我们从指标集和评估集的定义入手,以识别风险因素,并为专家提供评估依据。
指标集 U = ( U 1 , U 2 , ⋯ , U i ) ; i ∈ ( 1 , 2 , 3 , ⋯ , S ) ;评估集 V = ( 1 , 2 , 3 , 4 , 5 )
2) 权重向量A的测定
根据方法一节中介绍的方法,确定权重向量A并进行一致性测试。基于综合评价理论,确定评价矩阵具有可接受的不一致性。因此,矩阵是理性的。
比较矩阵:
A i = [ 1 4.038 3.781 3.990 0.248 1 2.800 3.562 0.264 0.357 1 3.562 0.251 0.281 0.281 1 ] ;
A 1 i = [ 1.000 3.762 3.171 3.514 0.266 1.000 2.110 2.410 0.315 0.474 1.000 2.314 0.285 0.415 0.432 1.000 ] ;
A 2 i = [ 1.000 3.810 3.781 3.848 0.263 1.000 2.705 2.539 0.264 0.370 1.000 2.571 0.260 0.394 0.389 1.000 ] ;
A 3 i = [ 1.000 2.600 0.385 1.000 ] ;
A 4 i = [ 1.000 1.908 1.967 0.524 1.000 2.505 0.508 0.399 1.000 ] ; i = 1 , 2 , 3 , 4.
利用Excel计算特征值和相应的特征向量
w = [ 0.5267 0.2423 0.1468 0.0739 ] ;
W 1 = [ 0.5230 0.2216 0.1576 0.0977 ] ;
W 2 = [ 0.5608 0.2382 0.1456 0.0918 ] ;
W 3 = [ 0.7222 0.2778 ] ;
W 4 = [ 0.4801 0.3383 0.1816 ] 。
一致性检验:
CR = 0.282 , CR 1 i = 0.0653 , CR 2 i = 0.0880 , CR 3 i = 0 , CR 4 i = 0.0760 ,结果 CR < 0.1 ,通过一致性检验。
3) 模糊评价矩阵的确定:
根据15位专家给出的分数(表6~8),构建了模糊评价矩阵 D = ( d i j ) m * n , m = 15 , n = 13 。
风险因素 | 风险很低 | 风险较低 | 风险一般 | 风险较高 | 风险很高 |
---|---|---|---|---|---|
Z1 | 0 | 2 | 4 | 5 | 4 |
Z2 | 0 | 2 | 8 | 5 | 0 |
Z3 | 0 | 0 | 10 | 5 | 0 |
Z4 | 1 | 0 | 2 | 8 | 4 |
Z5 | 0 | 1 | 3 | 9 | 2 |
Z6 | 0 | 1 | 2 | 12 | 0 |
Z7 | 1 | 0 | 4 | 8 | 2 |
Z8 | 0 | 0 | 1 | 9 | 5 |
Z9 | 0 | 0 | 3 | 10 | 2 |
Z10 | 1 | 1 | 8 | 5 | 0 |
Z11 | 0 | 2 | 2 | 11 | 0 |
Z12 | 1 | 2 | 3 | 9 | 0 |
Z13 | 1 | 5 | 5 | 4 | 0 |
表6. 专家对于住房反向抵押贷款的打分
风险因素 | 风险很低 | 风险较低 | 风险一般 | 风险较高 | 风险很高 |
---|---|---|---|---|---|
Z1 | 2 | 2 | 5 | 4 | 2 |
Z2 | 0 | 4 | 6 | 5 | 0 |
Z3 | 0 | 0 | 8 | 7 | 0 |
Z4 | 1 | 1 | 2 | 8 | 3 |
Z5 | 0 | 1 | 2 | 10 | 2 |
Z6 | 0 | 0 | 8 | 7 | 0 |
Z7 | 1 | 0 | 5 | 7 | 2 |
Z8 | 0 | 1 | 4 | 8 | 2 |
Z9 | 0 | 0 | 3 | 11 | 1 |
Z10 | 1 | 4 | 5 | 4 | 1 |
Z11 | 0 | 2 | 4 | 7 | 2 |
Z12 | 1 | 2 | 5 | 6 | 1 |
Z13 | 1 | 7 | 3 | 4 | 0 |
表7. 专家对于租房养老的打分
风险因素 | 风险很低 | 风险较低 | 风险一般 | 风险较高 | 风险很高 |
---|---|---|---|---|---|
Z1 | 1 | 1 | 1 | 7 | 5 |
Z2 | 0 | 2 | 4 | 9 | 0 |
Z3 | 0 | 1 | 6 | 8 | 0 |
Z4 | 1 | 0 | 4 | 7 | 3 |
Z5 | 0 | 0 | 2 | 11 | 2 |
Z6 | 0 | 1 | 3 | 10 | 1 |
Z7 | 1 | 1 | 4 | 8 | 1 |
Z8 | 0 | 0 | 6 | 8 | 1 |
Z9 | 0 | 1 | 1 | 9 | 4 |
Z10 | 1 | 2 | 6 | 6 | 0 |
Z11 | 0 | 2 | 4 | 8 | 1 |
Z12 | 1 | 2 | 5 | 6 | 1 |
Z13 | 2 | 6 | 4 | 3 | 0 |
表8. 专家对于售房养老的打分
D = [ d 11 d 12 ⋯ d 1 n d 21 d 22 ⋯ d 2 n ⋮ ⋮ ⋱ ⋮ d m 1 d m 2 ⋯ d m n ] 。
4) 灰色等级的确定、隶属函数和灰色决策的构建
矩阵和一系列计算后的结果
根据“方法论”一节中介绍的方法,15名无量纲处理专家给出了分数,然后得出如下结果:
R 1 i = [ 0.0000 0.1333 0.2667 0.3333 0.2667 0.0000 0.1333 0.5333 0.3333 0.0000 0.0000 0.0000 0.6667 0.3333 0.0000 0.0667 0.0000 0.1333 0.5333 0.2667 ] ;
R 2 i = [ 0.1353 0.0667 0.0000 0.1333 0.5333 0.0575 0.0000 0.0667 0.2000 0.6000 0.0351 0.0000 0.0667 0.1333 0.8000 0.0221 0.0667 0.0000 0.2667 0.5333 ] ;
R 3 i = [ 0.0000 0.0000 0.0667 0.6000 0.3333 0.0000 0.0000 0.2000 0.6667 0.1333 ] ;
R 4 i = [ 0.0667 0.0667 0.5333 0.3333 0.0000 0.0000 0.1333 0.1333 0.7333 0.0000 0.0667 0.1333 0.2000 0.6000 0.0000 ] , i = 1 , 2 , 3 , 4.
R 1 = B 1 = [ 0.0065 0.0993 0.3758 0.3529 0.1655 ] ;
R 2 = B 2 = [ 0.1029 0.0463 0.0272 0.1770 0.6465 ] ;
R 3 = B 3 = [ 0.0000 0.0000 0.1037 0.6185 0.2778 ] ;
R 4 = B 4 = [ 0.0441 0.1013 0.3375 0.5171 0.0000 ] ;
B = [ 0.0303 0.0723 0.2524 0.3658 0.2792 ] ;
C = B * V = [ 0.0065 0.0993 0.3758 0.3529 0.1655 0.1029 0.0463 0.0272 0.1770 0.6465 0.0000 0.0000 0.1037 0.6185 0.2778 0.0441 0.1013 0.3375 0.5171 0.0000 ] * [ 1 2 3 4 5 ] = [ 3.5716 4.2179 4.1741 3.3275 ] ;
E = W * C = 3.7912 。
反向抵押的最终得分为3.7912,表明反向抵押介于普通风险和略高风险之间。
5) 租房养老和售房养老的风险评估
鉴于该风险评估模型的广泛应用,我们还尝试将租房养老和售房养老作为对照进行评分。详细的处理程序如下:
R 1 i = [ 0.1333 0.1333 0.3333 0.2667 0.1333 0.0000 0.2667 0.4000 0.3333 0.0000 0.0000 0.0000 0.5333 0.4667 0.0000 0.0667 0.0667 0.1333 0.5333 0.2000 ] ;
R 2 i = [ 0.0000 0.0667 0.1333 0.6667 0.1333 0.0000 0.0000 0.5333 0.4667 0.0000 0.0667 0.0000 0.3333 0.4667 0.1333 0.0000 0.0667 0.2667 0.5333 0.1333 ] ;
R 3 i = [ 0.0000 0.0000 0.2000 0.7333 0.0667 0.0667 0.2667 0.3333 0.2667 0.0667 ] ;
R 4 i = [ 0.0000 0.1333 0.2667 0.4667 0.1333 0.0667 0.1333 0.3333 0.4000 0.0667 0.0667 0.4667 0.2000 0.2667 0.0000 ] ;
R 1 = B 1 = [ 0.0763 0.1354 0.3601 0.3390 0.0893 ] ;
R 2 = B 2 = [ 0.0094 0.0420 0.2652 0.5808 0.1027 ] ;
R 3 = B 3 = [ 0.0185 0.0741 0.2370 0.6037 0.0667 ] ;
R 4 = B 4 = [ 0.0347 0.1939 0.2771 0.4078 0.0866 ] ;
B = [ 0.0482 0.1078 0.3124 0.4426 0.0890 ] ;
C = B * V = [ 0.0763 0.1354 0.3601 0.3390 0.0893 0.0094 0.0420 0.2652 0.5808 0.1027 0.0185 0.0741 0.2370 0.6037 0.0667 0.0347 0.1939 0.2771 0.4078 0.0866 ] * [ 1 2 3 4 5 ] = [ 3.2297 3.7255 3.6259 3.3178 ] ;
E = W * C = 3.4164 。
R 1 i = [ 0.0667 0.0667 0.0667 0.4667 0.3333 0.0000 0.1333 0.2667 0.6000 0.0000 0.0000 0.0667 0.4000 0.5333 0.0000 0.0667 0.0000 0.2667 0.4667 0.2000 ] ;
R 2 i = [ 0.0000 0.0000 0.1333 0.7333 0.1333 0.0000 0.0667 0.2000 0.6667 0.0667 0.0667 0.0667 0.2667 0.5333 0.0667 0.0000 0.0000 0.4000 0.5333 0.0667 ] ;
R 3 i = [ 0.0000 0.0667 0.0667 0.6000 0.2667 0.0667 0.1333 0.4000 0.4000 0.0000 ] ;
R 4 i = [ 0.0000 0.1333 0.2667 0.5333 0.0667 0.0667 0.1333 0.3333 0.4000 0.0667 0.1333 0.4000 0.2667 0.2000 0.0000 ] ;
R 1 = B 1 = [ 0.0414 0.0749 0.1831 0.5067 0.1939 ] ;
R 2 = B 2 = [ 0.0094 0.0247 0.1910 0.6722 0.1027 ] ;
R 3 = B 3 = [ 0.0185 0.0852 0.1593 0.5444 0.1926 ] ;
R 4 = B 4 = [ 0.0468 0.1817 0.2892 0.4177 0.0546 ] ;
B = [ 0.0306 0.0721 0.1894 0.5469 0.1610 ] ;
C = B * V = [ 0.0414 0.0749 0.1831 0.5067 0.1939 0.0094 0.0247 0.1910 0.6722 0.1027 0.0185 0.0852 0.1593 0.5444 0.1926 0.0468 0.1817 0.2892 0.4177 0.0546 ] * [ 1 2 3 4 5 ] = [ 3.7368 3.8343 3.8074 3.2616 ] ;
E = W * C = 3.7357 。
因此,经过处理后,租房养老金得分为3.4164,房屋出售养老金得分为3.7357。这反映了中国人更喜欢租房养老模式,而不是住房反向抵押和售房养老。
在中国,住房养老保险的量化风险评估是一项具有挑战性的任务。这主要是由于以下三个原因:1) 数量性质;2) 以房养老的复杂性;以及3) 决策过程中涉及的不确定性。由于三个原因,大多数现有的以房养老风险评估方法都受到限制。层次分析法可以处理定性和定量变量,并解决第一个挑战。然而,这种方法对于解决复杂抵押的实际应用中的不确定性不是最佳的。因此,我们在模型中引入了模糊逻辑,这是一种广泛使用的将不确定因素纳入决策过程的方法。通过层次分析法确定各指标评价体系的权重,构建层次分析法模糊综合评价,并结合模糊综合评价形成选择(决策)模型。与其他模糊AHP模型相比,AHP模糊综合评价在计算能力方面简单,因为大多数非对称计算是使用EXCEL计算的。新的模型足够容易处理以房养老中激励因素的模糊性和复杂性。遗憾的是,这种AHP模糊综合评价存在一些局限性和弱点。风险可以通过几种方法来识别,包括头脑风暴、德尔福技术和访谈 [
住房养老一直是许多国家关注的焦点。本文采用AHP-fuzzy综合评价模型对以房养老风险进行定量评价。该评估系统可用于指导许多国家的住房换养老金风险管理,并提供有关住房反向抵押、租房养老和售房养老的有用信息。
邱 琳. 我国基于AHP-Fuzzy模糊综合评判的以房养老风险评估Risk Assessment of House-for-Pension with AHP-Fuzzy Comprehensive Evaluation in China[J]. 运筹与模糊学, 2023, 13(03): 2241-2253. https://doi.org/10.12677/ORF.2023.133224
https://doi.org/10.1201/9781420004137
https://doi.org/10.3109/03639045.2012.699065
https://doi.org/10.1108/09566160210431088
https://doi.org/10.1111/j.1539-6924.1985.tb00159.x
https://doi.org/10.2991/ieesasm-16.2016.96
https://doi.org/10.1016/j.insmatheco.2015.03.026
https://doi.org/10.1057/gpp.2011.26
https://doi.org/10.1016/j.habitatint.2019.102073
https://doi.org/10.1016/j.insmatheco.2021.02.008
https://doi.org/10.1016/j.jjie.2004.03.003
https://doi.org/10.1016/j.jfs.2016.10.001
https://doi.org/10.3846/tede.2022.16617
https://doi.org/10.3109/03639045.2012.744417