stata 计算基尼系数 igini

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

使用dasp 中的igini计算基尼系数，group可分类计算。但是是group中只能用一个变量。igini不能喝by连用。我用的是chns数据，就有省份又有年份还有城乡（urban=0、urban=1）。现在我想分年份分城乡分省份分别计算gini系数，该怎么做？我用最笨的办法把数据分割成不同文件，这个办法太笨了，有没有什么简单的办法，使用命令的

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享1 收藏52 回帖

关键词：Stata tata 基尼系数 CHNS数据 Group 基尼系数

ssc install ineqdeco
[quietly] ineqdeco var [fweight=...], by(group)标量里
r(gini)是总体基尼
r(gini_k)为每组的gini
还列出其它熵值GE(..)系列
GE(1)是Theil index of inequality，比Gini强，
它（总值Total Theil）可以分解为组内不平等(within)和跨组间不平等(between)
gen gini=0
gen gini_urban=0
gen gini_rural=0
foreach ind of numlist 1990/2011{
        foreahc i of numlist 1/20{
                qui ineqdeco income if year==`ind' & provid==`i', by(urban)
                replace gini=r(gini) if year==`ind' & provid==`i'
                replace gini_urban=r(gini_1) if year==`ind' & provid==`i'
                replace gini_rural=r(gini_2) if year==`ind' & provid==`i'
        }
}

要是有数据就好了，可以更具体些。

learning.

谢谢

monstersivle 发表于 2012-12-28 09:15
ssc install ineqdeco
[quietly] ineqdeco var [fweight=...], by(group)标量里
r(gini)是总体基尼

大侠，我按您的指点，写了一段小程序，但执行不了，一直显示无对应观察值，我收集了1999-2010年鉴的小学生均公用经费pbp，分中东西三个地区region3来考察各个地区gini系数变化的时序趋势，程序如下：
*没有成功执行
gen gini=0
gen gini_c=0
gen gini_e=0
gen gini_w=0
foreach ind of numlist 1999/2010{
        foreach i of numlist 1/12{
                qui ineqdeco pbp if year==`ind' & code==`i', by(region3)
                replace gini=r(gini) if year==`ind' & code==`i'
                replace gini_c=r(gini_1) if year==`ind' & code==`i'
                replace gini_e=r(gini_2) if year==`ind' & code==`i'
                replace gini_w=r(gini_3) if year==`ind' & code==`i'
        }
}
但是，显示“no observations”

monstersivle 发表于 2012-12-28 09:15
ssc install ineqdeco
[quietly] ineqdeco var [fweight=...], by(group)标量里
r(gini)是总体基尼

不好意思，code是省份的代码

monstersivle 发表于 2012-12-28 09:15
ssc install ineqdeco
[quietly] ineqdeco var [fweight=...], by(group)标量里
r(gini)是总体基尼

wealth        P141        PROVINCE
83000        6        11
232000        6        11
111000        6        11
20000        5        11
300000        6        12
560000        5        11
123000        6        12
上面分别是收入，职业类型和省份。
我想按照职业类型和省份计算gini系数，两个职业，两个省份，这样就可以分成四类，得到四个gini系数。
按照您给出这个命令，总是出现错误提示“没有观测值”，是为什么？

ranbo7856 发表于 2014-1-18 11:28
wealth        P141        PROVINCE
83000        6        11
232000        6        11

已解决。
原来是由于，没有5,12这样的观测值，命令没有任何问题。
由于我是省略，真实的数据没有这样情况，直接可出结果。

另外，之所以报错，是因为ineqdeco计算时必须有观测值。解决办法，首先生成一个分组变量。
gen code=0
foreach ind of numlist 1/9 11 .{
     foreach i of numlist 11 14 21 32 34 41 42 44 50 51 53 62 {
                foreach j of numlist 1/5 {
                                        replace code=`ind'`i'`j' if (P141==`ind' & PROVINCE==`i' & age==`j')
                          
       }
           }
}


*再把此分组变量变成连续的分组编号
egen code_new = group(code)
sum code_new

*查看最大值

gen gini=0
foreach ind of numlist 1/max值{
                         qui ineqdeco wealth if (code_new==`ind')
                replace gini=r(gini) if (code_new==`ind')

}



这个是笨办法，有可能ineqdeco有现成选项避免此情况。仅供参考，也希望知道的同学可以指导。

基尼系数是靠官员们说的，不用算的。