Сравните индексы концентрации для случайных портфелей

Этот пример показывает, как моделировать случайные портфели с различными дистрибутивами и сравнить их индексы концентрации. В целях рисунка используются логарифмически нормальное и распределение Weibull. Параметры распределения выбраны произвольно, чтобы получить подобную область значений значений для обоих случайных портфелей.

Сгенерируйте случайные портфели с различными дистрибутивами.

rng('default'); % for reproducibility
PLgn = lognrnd(1,1,1,300);
PWbl = wblrnd(2,0.5,1,300);

Отобразите самое большое моделируемое значение ссуды.

fprintf('\nLargest loan Lognormal: %g\n',max(PLgn));
Largest loan Lognormal: 97.3582
fprintf('Largest loan Weibull: %g\n',max(PWbl));
Largest loan Weibull: 91.5866

Постройте гистограммы портфеля.

figure;
histogram(PLgn,0:5:100)
hold on
histogram(PWbl,0:5:100)
hold off
title('Random Loan Histograms')
xlabel('Loan Amount')
ylabel('Frequency')
legend('Lognormal','Weibull')

Вычислите и отобразите меры по концентрации.

ciLgn = concentrationIndices(PLgn,'ID','Lognormal');
ciWbl = concentrationIndices(PWbl,'ID','Weibull');
disp([ciLgn;ciWbl])
        ID            CR          Deciles        Gini         HH          HK           HT          TE   
    ___________    ________    _____________    _______    ________    _________    _________    _______

    "Lognormal"    0.066363    [1x11 double]     0.5686    0.013298    0.0045765    0.0077267    0.66735
    "Weibull"      0.090152    [1x11 double]    0.72876    0.020197    0.0062594     0.012289     1.0944
ProportionLoans = 0:0.1:1;
figure;
area(ProportionLoans',[ciWbl.Deciles; ciLgn.Deciles-ciWbl.Deciles; ProportionLoans-ciLgn.Deciles]')
axis([0 1 0 1])
legend('Weibull','Lognormal','Diversified','Location','NorthWest')
title('Lorenz Curve (by Deciles)')
xlabel('Proportion of Loans')
ylabel('Proportion of Value')

Смотрите также

Связанные примеры

Больше о