В этом примере показано, как выполнить расчеты ECL с помощью моделируемых данных ссуды, данных макро-сценария и существующей модели вероятности дефолта (PD) во время жизни.
Загрузите данные ссуды, готовые к предсказанию, данные макроперехода и соответствующие вероятности сценария.
load DataPredictLifetime.mat
disp(LoanData) ID ScoreGroup YOB Year
____ _____________ ___ ____
1304 "Medium Risk" 4 2020
1304 "Medium Risk" 5 2021
1304 "Medium Risk" 6 2022
1304 "Medium Risk" 7 2023
1304 "Medium Risk" 8 2024
1304 "Medium Risk" 9 2025
1304 "Medium Risk" 10 2026
2067 "Low Risk" 7 2020
2067 "Low Risk" 8 2021
2067 "Low Risk" 9 2022
2067 "Low Risk" 10 2023
disp(head(MultipleScenarios,10))
ScenarioID Year GDP Market
__________ ____ ____ ______
"Severe" 2020 -0.9 -5.5
"Severe" 2021 -0.5 -6.5
"Severe" 2022 0.2 -1
"Severe" 2023 0.8 1.5
"Severe" 2024 1.4 4
"Severe" 2025 1.8 6.5
"Severe" 2026 1.8 6.5
"Severe" 2027 1.8 6.5
"Adverse" 2020 0.1 -0.5
"Adverse" 2021 0.2 -2.5
disp(ScenarioProbabilities)
Probability
___________
Severe 0.1
Adverse 0.2
Baseline 0.3
Favorable 0.2
Excellent 0.2
load LifetimeChampionModel.mat
disp(pdModel) Probit with properties:
ModelID: "Champion"
Description: "A sample model used as champion model for illustration purposes."
Model: [1x1 classreg.regr.CompactGeneralizedLinearModel]
IDVar: "ID"
AgeVar: "YOB"
LoanVars: "ScoreGroup"
MacroVars: ["GDP" "Market"]
ResponseVar: "Default"
Для расчетов ECL требуются только предельные PD. Однако сначала можно визуализировать PD в течение жизни.
CompanyIDChoice ="1304"; CompanyID = str2double (CompanyIDChoice); IndCompany = LoanData.ID = = CompanyID; Годы = LoanData.Year (IndCompany); NumYears = длина (годы); ScenarioID = уникальный (MultipleScenarios.ScenarioID,'stable'); NumScenarios = длина (ScenarioID); LifetimePD = нули (NumYears, NumScenarios); for ii = 1: NumScenarios IndScenario = MultipleScenarios.ScenarioID = = ScenarioID (ii); данные = joint (LoanData (IndCompany,:), MultipleScenarios (IndScenario,:)); LifetimePD (:, ii) = predictLifetime (pdModel, данные); end график (годы, LifetimePD) xticks (Годы) сетка on xlabel ('Year') ylabel ('Lifetime PD') заголовок ('Lifetime PD By Scenario') legend (ScenarioID,'Location','best')

Строго говоря, расчет ECL требуется пожизненная модель PD, пожизненная модель LGD и пожизненная модель EAD, плюс сценарии, вероятности сценария и эффективная процентная ставка.
Для простоты этот пример принимает постоянные LGD и EAD модели и заданная процентная ставка.
LGD = 0.55; EAD = 100000; EffRate = 0.045; CompanyIDChoice ="1304"; CompanyID = str2double (CompanyIDChoice); IndCompany = LoanData.ID = = CompanyID; Годы = LoanData.Year (IndCompany); NumYears = длина (годы); ScenarioID = уникальный (MultipleScenarios.ScenarioID,'stable'); NumScenarios = длина (ScenarioID); MarginalPD = нули (NumYears, NumScenarios); for ii = 1: NumScenarios IndScenario = MultipleScenarios.ScenarioID = = ScenarioID (ii); данные = joint (LoanData (IndCompany,:), MultipleScenarios (IndScenario,:)); MarginalPD (:, ii) = predictLifetime (pdModel, данные,'ProbabilityType','marginal'); end DiscTimes = Years-Years (1) + 1; DiscFactors = 1 ./( 1 + EffRate). ^ DiscTimes; ProbScenario = СценарийВероятности.Вероятность; ECL_t_s = (MarginalPD * LGD * EAD). * Факторы диска;% ECL by year and scenario ECL_s = сумма (ECL_t_s );% ECL total by scenario ECL = ECL_s*ProbScenario; % ECL weighted average over all scenarios % Arrange ECL data for display in table format % Append ECL total per scenario and scenario probabilities ECL_Disp = array2table ([ECL _ t _ s; ECL_s; ProbScenario ']); ECL_Disp.Properties.VariableNames = ScenarioID; ECL_Disp.Properties.RowNames = [strcat ("ECL ", строка (Годы ));"ECL total"; "Probability"]; disp (ECL_Disp)
Severe Adverse Baseline Favorable Excellent
______ _______ ________ _________ _________
ECL 2020 595.58 507.16 430.44 364.11 306.97
ECL 2021 394.24 349.95 310.02 274.11 241.9
ECL 2022 235.53 215.4 196.75 179.5 163.57
ECL 2023 143.05 135.23 127.75 120.59 113.77
ECL 2024 85.219 83.517 81.816 80.118 78.429
ECL 2025 51.346 51.514 51.665 51.798 51.917
ECL 2026 33.162 33.271 33.368 33.454 33.531
ECL total 1538.1 1376 1231.8 1103.7 990.08
Probability 0.1 0.2 0.3 0.2 0.2
fprintf('Lifetime ECL for company %s is: %g\n',CompanyIDChoice,ECL)Lifetime ECL for company 1304 is: 1217.32
fitLifetimePDModel | Logistic | modelAccuracy | modelDiscrimination | predict | predictLifetime | Probit