stata测量健康不平等的方法，集中指数估算、健康不平等的分解

在文献中，健康不平等的测算包括健康集中指数、健康不平等的分解、基尼系数、广义熵值数的泰尔指数阿特金森指数、Shapley和RIF-I-OLS（只是记录可能的方法，并不一定适用于每个人的文章，按自己要求采用方法才是重要的）。
在学习这些方法过程中，记录下来。

第一：健康集中指数：如何算出集中指数和集中曲线

help conindex
***哦豁，没有这个命令，安装

ssc install conindex
**我可真棒，安装成功。

help conindex
***看看这个语句的命令是咋写的，原理都是公式，放过它，放过自己。

conindex health, rankvar(age) truezero
**先随便写一个，找到感觉

conindex health, rankvar(hhincome) truezero compare(hukou) 
///测算城乡之间的健康不均等：health为健康状况，hhincome为家庭收入，hukou为户籍, truezero是广义集中度指数

****进一步演化：wagstaff和erreygers获得归一化浓度指数
gen B_health = invnormal(uniform())>0
conindex B_health, rankvar(hhincome) truezero bounded limits(0 1) wagstaff rankvar((hukou)
conindex B_health, rankvar(hhincome) truezero bounded limits(0 1) erreygers rankvar((hukou)

conindex health, rankvar(hhincome) truezero compare(hukou) graph ytitle(“健康”) xtitle(“收入情况”)
我发现我用graph语句也画不出图像，出错：command lorenz is unrecognized
这样，我就继续help conindex看看graph语句，发现因为 lorenz这个命令我没有安装，因此

ssc install lorenz
conindex health, rankvar(hhincome) truezero compare(hukou)  graph ytitle("健康") xtitle("收入情况")
****这样集中指数估算和画图都完成啦，我们可太棒啦，一会会就学会了新的方法~~——~~

第二，学习两个健康不平等的分解方法
1.shapley分解方法

ssc install shapley2

 reg health pension age hukou gender exerice yrs_edu num_hhmember ln_ave_hhincome_1  per_GDP 
***health 是健康，pension是养老金，age hukou gender yrs_edu exerice是个体特征，num_hhmember ln_ave_hhincome_1 是家庭特征，per_GDP 是省份层面的变量

shapley2,stat(r2) group(pension age hukou gender exerice yrs_edu num_hhmember ln_ave_hhincome_1  per_GDP )
too few variables specified
***出现问题，解决问题。
help shapley2
drop _mysample _est_myreg
reg health pension age hukou gender yrs_edu exerice num_hhmember ln_ave_hhincome_1  per_GDP 
shapley2, stat(r2)
/*
 shapley2, stat(r2)
Factor     | Shapley value |  Per cent 
           |  (estimate)   | (estimate)
-----------+---------------+-----------+
insure1    |  0.00055      |    2.25 % |
age        |  0.00004      |    0.17 % |
hukou      |  0.00497      |   20.19 % |
gender     |  0.00200      |    8.13 % |
yrs_edu    |  0.00210      |    8.51 % |
exerice    |  0.00078      |    3.17 % |
num_hhme~r |  0.00439      |   17.84 % |
ln_ave_h~1 |  0.00965      |   39.17 % |
per_GDP    |  0.00014      |    0.57 % |
-----------+---------------+-----------+
TOTAL      |  0.02463      |  100.00 % |
-----------+---------------+-----------+
*/
**可以看出贡献的百分比，哇塞，我们太棒啦，又解锁了新技能。
**group是分组，可以按需要多尝试

2.Oaxaca分解的stata实现及解读

ssc install oaxaca
help oaxaca 
oaxaca health_self_report insure1 age  gender yrs_edu exerice num_hhmember, by(hukou)

***学会三个方法，并掌握上述健康不平等测算，一篇论文实证部分可以算是完成了。后续有用到其他的，再更新。方法都不难，多尝试！