求如何画此Lorenz Curve

做了有关收入的调查问卷，想做一个rural和urban收入对比的Lorenz Curve看收入差距，求大神帮忙看看。CSV不能上传所以打包上传了。(100个rural和400个urban和8组收入区间)，

见图

1. keep annualincome area
   
2. 
   
3. // 将分类数据编码为数值
   
4. gen income=10000
   
5. replace income=20000 if annualincome=="10000-30000"
   
6. replace income=40000 if annualincome=="30000-50000"
   
7. replace income=65000 if annualincome=="50000-80000"
   
8. replace income=100000 if annualincome=="80000-120000"
   
9. replace income=135000 if annualincome=="120000-150000"
   
10. replace income=165000 if annualincome=="150000-180000"
    
11. replace income=200000 if annualincome=="200000 or above"
    
12. 
    
13. 
    
14. // 计算累计人口，累计总收入
    
15. sort area income
    
16. 
    
17. by area income: gen num=_N
    
18. 
    
19. collapse (sum) sum_income = income (count) num = income, by(area annualincome)
    
20. 
    
21. by area: egen total_income = total(sum_income)
    
22. by area: egen total_people = total(num)
    
23. by area: gen cum_share_people = 100*sum(num)/total_people
    
24. by area: gen cum_share_income = 100*sum(sum_income)/total_income
    
25. 
    
26. keep area cum_share_people cum_share_income
    
27. 
    
28. 
    
29. // 画图前补上（0,0）这个点
    
30. input
    
31. "Rural" 0 0
    
32. "Urban" 0 0
    
33. end
    
34. 
    
35. sort area cum_share_people
    
36. 
    
37. 
    
38. // 画线图
    
39. twoway  line  cum_share_income cum_share_people if area=="Rural"   ///
    
40. ||  line  cum_share_income cum_share_people if area=="Urban", legend(label(1 "Rural") label(2 "Urban"))
    
41. 
    
42. // 如果希望线条平滑
    
43. twoway  qfit  cum_share_income cum_share_people if area=="Rural"   ///
    
44. ||  qfit  cum_share_income cum_share_people if area=="Urban", legend(label(1 "Rural") label(2 "Urban"))
    
45. 
    
46. 
    

复制代码

deem 发表于 2017-8-21 05:58
见图

哇，大神能顺便帮忙解释下此图么

shangwuqiwen 发表于 2017-8-21 07:27
哇，大神能顺便帮忙解释下此图么

每个（x,y）点表示，财富最低的x%人群拥有y%的财富，用下面的数据举例，在农村地区收入最低的29%的人仅占有9.446%的财富。曲线越向右下弯曲就说明不平等越严重。上图显示农村地区曲线在下方，所以农村地区不平等更加严重。

所用数据:

1. area        cum_share_people        cum_share_income
   
2. Rural        0        0
   
3. Rural        29        9.446254
   
4. Rural        66        33.55049
   
5. Rural        67        37.94788
   
6. Rural        69        48.69707
   
7. Rural        88        73.45277
   
8. Rural        99        96.74267
   
9. Rural        100        100
   
10. Urban        0        0
    
11. Urban        9        1.329394
    
12. Urban        31.5        7.976367
    
13. Urban        38        20.93796
    
14. Urban        40.5        27.03102
    
15. Urban        48        49.18759
    
16. Urban        68        61.00443
    
17. Urban        84        76.36633
    
18. Urban        100        100

复制代码

deem 发表于 2017-8-21 08:05
每个（x,y）点表示，财富最低的x%人群拥有y%的财富，用下面的数据举例，在农村地区收入最低的29%的人仅占 ...

谢谢解释，非常清楚。还有个问题问下哈，如果根据分布来看
这个是总结的图表，rural的收入普遍偏低，而urban的相对平均分散，是不是也说明rural不平等，而urban相对平等些？还是我一直以为的是rural的收入都低所以差距不大所以相对平等，而urban的收入差距大所以不平等....我被图表和自己的逻辑弄混了，帮我分析下，谢谢。

shangwuqiwen 发表于 2017-8-21 08:26
谢谢解释，非常清楚。还有个问题问下哈，如果根据分布来看
这个是总结的图表，rural的收入普遍偏低，而u ...

和绝对收入无关，就是看图就行了。一般讲的基尼系数就是曲线和45度线围成区域的面积。

deem 发表于 2017-8-21 08:41
和绝对收入无关，就是看图就行了。一般讲的基尼系数就是曲线和45度线围成区域的面积。

哦，好的，谢谢！

deem 发表于 2017-8-21 08:05
每个（x,y）点表示，财富最低的x%人群拥有y%的财富，用下面的数据举例，在农村地区收入最低的29%的人仅占 ...

大神 想问下是否会用DASP插件画洛伦兹曲线 如何把纵轴设置为各省人口比重 横轴设置为各省面积比重呢