stata如何计算空间基尼系数？

潜水的猫

18808

收藏 2013-10-09

请教各位，想按照下图文献中的公式计算空间基尼系数rg，不知道应该如何写code，

我有尝试用这个帖子里说的算，但是出现了负值（应该是0到1之间）：https://bbs.pinggu.org/thread-494345-3-1.html
sort province income
by province: egen t=sum(_n*income)
by province: egen s=sum(income)
by province: g gini=2*t/(_N*s)-1-1/_N

请高手指教！

扫码加我拉你入群

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

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

全部回复

夏目贵志

2013-10-10 02:05:52

发一些数据上来，不然这个问题不好回答

扫码加我拉你入群

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

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

潜水的猫

2013-10-10 14:46:22

夏目贵志发表于 2013-10-10 02:05
发一些数据上来，不然这个问题不好回答

注意力完全被娘口三三吸引。。。

回来。。。。

数据大概长这样（简化为有2个地区、2年的数据）：n是地区行业数量：3（共有11、12、13三个行业）；
i,j属于｛1，2，3｝；
x应该是行业｛11，12，13｝各自在每年(year)每地区(prvcid)三个行业总值中的份额，也就是根据下表中的prod来计算的；
u应该是在每年(year)每省(prvcid)x的均值。

year	indid	prvcid	prod
1999	11	110000	12356
1999	11	120000	56423
1999	12	110000	5698
1999	12	120000	7842
1999	13	110000	99632
1999	13	120000	10254
2000	11	110000	16897
2000	11	120000	66812
2000	12	110000	10235
2000	12	120000	5996
2000	13	110000	78985
2000	13	120000	33000

然后不会在stata中写|xi-xj|。。。和同学讨论了，rg的公式和底下参考的帖子给出的这个是不同的表达式，所以不能直接用。。。

谢谢猫老师！

扫码加我拉你入群

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

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

夏目贵志

2013-10-10 23:45:00

猫老师很开心。托我送你这个礼物：

复制代码

你给的示例数据算出来的结果是

复制代码

扫码加我拉你入群

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

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

夏目贵志

2013-10-10 23:50:16

猫老师还说，你的这个公式有问题，按照这个公式，rg是不会等于1的。假定x1=1，x2=0，n=2。求和的结果是2，按照公式rg=(2*n^2*u)^(-1)*2=2/(8*u)。这里如果要rg=1，那么u必须等于1/4。显然在只有两个行业，总产值是1的情况下，u无论如何也不该是1/4的。

猫老师要喝酒吃包子了。。。

扫码加我拉你入群

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

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

潜水的猫

2013-10-11 20:44:44

夏目贵志发表于 2013-10-10 23:50
猫老师还说，你的这个公式有问题，按照这个公式，rg是不会等于1的。假定x1=1，x2=0，n=2。求和的结果是2，按 ...

= =!!!猫老师在STATA界也是kamisama一般的存在么。。。

不过猫老师的例子中，如果x1=1,x2=0，其实不是n=2，而是n=1只有一个行业的情况了，求和的结果。。。我有点晕了...

扫码加我拉你入群

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

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

点击查看更多内容…

潜水的猫

2013-10-11 20:48:32

另外想请教猫老师一下，实际中我的数据里indid里面有11A这种形式，造成indid变成了string，有没有办法把它变成数值的形式。。=================
后来我回到txt源文件把indid替换掉了。。。。

扫码加我拉你入群

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

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

夏目贵志

2013-10-11 22:04:49

潜水的猫发表于 2013-10-11 20:48
另外想请教猫老师一下，实际中我的数据里indid里面有11A这种形式，造成indid变成了string，有没有办法把它变 ...

替换掉也是可以的。或者就用encode命令。help encode就能看到用法了。

扫码加我拉你入群

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

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

夏目贵志

2013-10-11 22:12:32

潜水的猫发表于 2013-10-11 20:44
= =!!!猫老师在STATA界也是kamisama一般的存在么。。。

不过猫老师的例子中，如果x1=1,x2=0，其实不是 ...

如果某个地区产出是完全集中的，那么就表示除了一个行业之外，其他所有行业的产出都为0。即地区工业总产出等于这一个行业的产出。所以我给的那个例子里，n=2是没问题的。我说求和的结果rg公式里除了(1/(2n^2*u))的部分，即double summation的结果。

一般情况下，如果某地区有n个行业，x_1, ..., x_n，且除x_1外的所有x都为0，求和的结果就是x_1*n*(n-1)，将x_1 normalize成1，得n*(n-1)。带入：rg=(n*(n-1))/(2n^2*u)=(n-1)/(2*n*u)，可见u的取值无法使rg=1。

不过无所谓啦。等解决你的问题就好。