В этом примере показано, как выполнить расчеты ожидаемой потери кредита (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 только требуются крайние ФУНТЫ. Однако сначала можно визуализировать пожизненные ФУНТЫ.
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); данные = соединение (LoanData (IndCompany, :), MultipleScenarios (IndScenario, :)); LifetimePD (: ii) = predictLifetime (pdModel, данные); end график (годы, LifetimePD) xticks (годы) сетка on xlabel'Year') yLabel 'Lifetime PD') заголовок'Lifetime PD By Scenario') легенда (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); данные = соединение (LoanData (IndCompany, :), MultipleScenarios (IndScenario, :)); MarginalPD (: ii) = predictLifetime (pdModel, данные,'ProbabilityType','marginal'); end DiscTimes = годы лет (1) +1; DiscFactors = 1. / (1+EffRate).^DiscTimes; ProbScenario = ScenarioProbabilities.Probability; ECL_t_s = (MarginalPD*LGD*EAD).*DiscFactors; % 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 | predict | predictLifetime | modelDiscrimination | modelAccuracy | Logistic | Probit | Cox