stata 计算基尼系数 igini

43474

收藏 2012-12-24

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

扫码加我拉你入群

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

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

全部回复

monstersivle

2012-12-28 09:15:30

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'
}
}

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

扫码加我拉你入群

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

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

peyzf

2013-1-1 07:18:22

learning.

扫码加我拉你入群

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

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

holywill

2013-1-1 11:45:01

谢谢

扫码加我拉你入群

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

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

HM_CUHK

2013-9-20 22:11:13

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”

扫码加我拉你入群

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

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

HM_CUHK

2013-9-20 22:12:25

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

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

扫码加我拉你入群

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

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

点击查看更多内容…

ranbo7856

2014-1-18 11:28:59

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 14:24:34

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

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

扫码加我拉你入群

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

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

ranbo7856

2014-1-18 14:26:04

另外，之所以报错，是因为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有现成选项避免此情况。仅供参考，也希望知道的同学可以指导。

扫码加我拉你入群

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

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

doublethree

2014-1-24 05:22:54

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

扫码加我拉你入群

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

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

wangyong123aza

2014-4-25 15:00:56

想咨询楼主，stata中igini这种方法谁最先提出来？

扫码加我拉你入群

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

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

zhangyawei1027

2014-6-10 15:36:29

大侠，我是用CHNS数据2000、2004、2006、2009年的数据，用9个省的数据，分年份分城乡分身份算gini系数，我写的命令是下面，但是为什么总是出现没有观察值？急急急
gen gini=0
gen gini_urban=0
gen gini_rural=0
foreach ind of numlist 2000 2004 2006 2009{
      foreach i of numlist 1/9{
            qui ineqdeco C8 if wave==`ind' & T1==`i', by(urban)
            replace gini=r(gini) if wave==`ind' & T1==`i'
            replace gini_urban=r(gini_1) if wave==`ind' & T1==`i'
            replace gini_rural=r(gini_2) if wave==`ind' &T1==`i'
      }
}

扫码加我拉你入群

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

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

lucy_xue

2014-7-22 06:20:54

. ineqdeco avincome, by(famtype) w
怎么看不见基尼斯数？

急用！谢谢！！

扫码加我拉你入群

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

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

tambei0610

2014-10-27 12:06:27

求救！
我的是分产业，分年度的销售额数据。想计算各年度和各产业的基尼系数。
我用了DASP，分组的话只能算出年度或产业的基尼系数，算不出来一起的。
按照这里各位大侠的指点，写了一小段程序，也是运行不出来
这是我的部分数据，大侠帮看看改怎么算
KSIC2 CODE year SALE
C10 2009 4.2006E+13
C10 2010 4.53942E+13
C10 2011 5.26003E+13
C10 2012 5.77099E+13
C10 2013 5.93993E+13
C11 2009 6.16927E+12
C11 2010 6.09334E+12
C11 2011 7.91286E+12
C11 2012 9.52084E+12
C11 2013 9.7327E+12
C12 2009 3.22536E+12
C12 2010 2.96616E+12
C12 2011 2.94959E+12
C12 2012 3.11312E+12
C12 2013 2.99509E+12