stata计算基尼系数-分地区分年份计算

***计算单个基尼系数（参考陈传波老师《stata十八讲》）
capture program drop gini //如果 gini 命令已存在，则删除，否则跳过该步
program gini //定义命令名为 gini
egen tinc=sum(1') //计算 gini 命令后第一个变量（收入支出等）的总和 egen tp=sum(2’) //计算 gini 命令后第二个变量（人口）的总和
gen m=1'/2’ //计算人均水平值 m
sort m //按人均水平值排序
gen gini=1-sum(2'/tp*(2*sum(1’/tinc)-`1’/tinc)) //计算基尼系数
dis gini[_N] //显示基尼系数
end //程序结束

gini 变量1(收入变量) 变量2(人数变量)

*分地区分年份计算基尼系数（参考陈传波老师《stata十八讲》后进一步拓展）
capture program drop gini
program gini
bysort cid10 year:egen tinc=sum(1') bysort cid10 year:egen tp=sum(2’)
gen m=1'/2’
sort cid10 year m
bysort cid10 year:gen 3'=1-sum(2’/tp(2sum(1'/tinc)-1’/tinc))
bysort cid10 year:replace 3'= 3’[_N] //显示基尼系数
end

gini 变量1(收入变量) 变量2(人数变量)