请教各位，如何计算各省以及各省分城乡的基尼系数，数据结构为prov urban income，其中prov为省行政编号11，12，……，65，urban为0（农村）和1（城镇），income就是收入的连续变量。
我尝试过用：https://bbs.pinggu.org/thread-494345-3-1.html
sort prov income
by prov: egen t=sum(_n*income)
by prov: egen s=sum(income)
g gini=2*t/(_N*s)-1-1/_N
但结果都是负值。
也尝试过：https://bbs.pinggu.org/thread-2161878-1-1.html
foreach i of numlist 11 12 13 14 15 21 22 23 31 32 33 34 35 36 37 41 42 43 44 45 46 50 51 52 53 54 61 62 63 64 65 {
quietly ineqdeco income if prov==`i',by(urban)
replace gini_prov=r(gini) if prov==`i' //各省的基尼系数
replace gini_prov_rural=r(gini_1) if prov==`i' //各省的农村基尼系数
replace gini_prov_urban=r(gini_2) if prov==`i' //各省的城镇基尼系数
}
但运行花了将近1个小时，而且只有gini_prov 和gini_prov_rural 是有数据的，gini_prov_urban 是缺失的。

求大神指点迷津！