matlab画基尼系数,Dagum基尼系数分解的MATLAB程序代码(更新)

很多同学在初次使用Dagum基尼系数法做研究时苦于找不到计算软件,手动计算又耗时费力。去年写论文时,参考Dagum(1997)论文中提供的算例,自己编写了一个MATLAB版的计算程序,在这里分享给大家。转载和使用请注明出处!

一、Dagum基尼系数分解简介

1、总体基尼系数

mathormula=%5Cmathrm%7BG%7D_%7B%5Cmathrm%7Bjj%7D%7D%3D%5Cfrac%7B1%7D%7B2%20%5Cbar%7By%7D_%7Bj%7D%7D%5Cleft(%5Csum_%7Bi%3D1%7D%5E%7Bn_%7Bj%7D%7D%20%5Csum_%7Br%3D1%7D%5E%7Bn_%7Bj%7D%7D%5Cleft%7Cy_%7Bj%20i%7D-y_%7Bj%20r%7D%5Cright%7C%5Cright)%20%2F%20n_%7Bj%7D%5E%7B2%7D

3、子群内差异对总体基尼系数的贡献

mathormula=%5Cmathrm%7BG%7D_%7B%5Cmathrm%7Bjh%7D%7D%3D%5Cleft(%5Csum_%7Bi%3D1%7D%5E%7Bn_%7Bj%7D%7D%20%5Csum_%7Br%3D1%7D%5E%7Bn_%7Bh%7D%7D%5Cleft%7Cy_%7Bj%20i%7D-y_%7Bh%20r%7D%5Cright%7C%5Cright)%20%2F%20n_%7Bj%7D%20n_%7Bh%7D%5Cleft(%5Cbar%7By%7D_%7Bj%7D%2B%5Cbar%7By%7D_%7Bh%7D%5Cright)

5、子群间净值差距对总体基尼系数的贡献度

mathormula=%5Cmathrm%7BG%7D_%7B%5Cmathrm%7Bt%7D%7D%3D%5Csum_%7Bj%3D1%7D%5E%7Bk%7D%20%5Csum_%7Bh%3D1%7D%5E%7Bj-1%7D%20%5Cmathrm%7BG%7D_%7B%5Cmathrm%7Bjh%7D%7D%5Cleft(P_%7Bj%7D%20S_%7Bh%7D%2BP_%7Bh%7D%20S_%7Bj%7D%5Cright)%5Cleft(1-D_%7Bj%20h%7D%5Cright)

其中:

mathormula=d_%7Bj%20h%7D%3D%5Cint_%7B0%7D%5E%7B%5Cinfty%7D%20d%20F_%7Bj%7D(y)%20%5Cint_%7B0%7D%5E%7By%7D(y-x)%20d%20F_%7Bh%7D(x)

bca653fd62eetm_campaign=shakespeare

输出结果1.png

注1:G_sub(1)为子群A基尼系数,G_sub(2)为子群B基尼系数。按变量输入顺序逐一对应。

注2:G_jh_1_2为A和B的子群间基尼系数,G_jh_1_3为A和C的子群间基尼系数。同样按输入变量顺序逐一对应。

注3:上述计算结果仅为单一年份的结果,要计算不同年份的基尼系数分解结果,需根据每年的子群数据逐一使用该程序计算。

注4:该程序支持3个以上的子群划分。例如存在A、B、C、D四个子群,且有B的均值>A的均值>D的均值>C的均值,则输入Dagum_Gini(B,A,D,C)即可。

四、特殊情形

1、子群的分布完全不重合:

例如,有A、B、C三个子群,子群C的所有元素均大于子群B的元素,子群B的所有元素均大于子群A的元素。此时,子群间的相对影响

mathormula=G_%7Bt%7D%3D0

测试代码:

A=[1:0.5:10]’;

B=[11:0.8:20]’;

C=[21:0.6:30]’;

%C的均值>B的均值>A的均值

Dagum_Gini(C,B,A);

输出结果:

mathormula=%5Coverline%20A%3D%5Coverline%20B%3D%5Coverline%20C,则,子群间相对影响

mathormula=G_%7Bnb%7D%3D0

测试代码:

a=magic(10);%生成魔方矩阵(每一列均值都相等)

A=a(:,1);

B=a(:,2);

C=a(:,3);

%A的均值=B的均值=C的均值

Dagum_Gini(A,B,C);

输出结果

bca653fd62eetm_campaign=shakespeare

输出结果3.png

相关资源:经典刻录软件尼禄Nero6.6.1.4中文安装版及序列号-其它工具类资源…

来源:董德生

声明:本站部分文章及图片转载于互联网,内容版权归原作者所有,如本站任何资料有侵权请您尽早请联系jinwei@zod.com.cn进行处理,非常感谢!

上一篇 2021年2月15日
下一篇 2021年2月15日

相关推荐