Stata - 分年度-行业计算销售额前四名的行业占比

作者：连玉君 (知乎 | 简书 | 码云)

Stata连享会 精彩推文1  || 精彩推文2

空间计量专题-海报.jpg

1. 问题

在 Stata 中如何分年度计算一个行业当中排名前四位的公司的总收入占整个行业的收入比重？

Note： 如果是算一般的赫芬达尔指数，可以使用连玉君老师编写的 hhi5 命令，很便捷。安装命令的方法为：在 Stata 命令窗口中执行 ssc install hhi5, replace 即可。

本例的情况稍微特殊一些，需要自己写命令，但核心思路与 hhi5 是一样的。有兴趣的读者可以使用 adoedit hhi5.ado 打开 hhi5 命令的原始 ado 文件，以便知晓内部计算思路和过程。

2. 解决方法

我们先自行模拟生成一份 Panel data 数据。自己操作时，换用自己的数据即可。

2.1 自行模拟生成一份数据

1. *---

2. *-A: Data Generation: simulation

3. clear

4. set obs 50 // N=50 家

5. gen id = _n

6. egen industry = cut(id), at(1(10)55) // 行业

7. expand 4 // T=4 年

8. bysort id: gen year = _n + 2010

9. xtset id year

10. gen sale = int(100+rnormal()*20)

生成的数据长这样：

1. . gsort year industry -sale

2. . list in 1/20, sepby(year industry) noobs

3. +-----------------------------+

4. | id industry year sale |

5. |-----------------------------|

6. | 10 1 2011 137 |

7. | 1 1 2011 119 |

8. | 6 1 2011 102 |

9. | 4 1 2011 99 |

10. | 3 1 2011 97 |

11. | 2 1 2011 97 |

12. | 9 1 2011 96 |

13. | 7 1 2011 85 |

14. | 8 1 2011 84 |

15. | 5 1 2011 82 |

16. |-----------------------------|

17. | 11 11 2011 135 |

18. | 15 11 2011 134 |

19. | 20 11 2011 133 |

20. | 12 11 2011 119 |

21. | 13 11 2011 117 |

22. | 16 11 2011 104 |

23. | 17 11 2011 101 |

24. | 19 11 2011 95 |

25. | 14 11 2011 93 |

26. | 18 11 2011 88 |

27. +-----------------------------+

2.2 计算目标变量

1. *---

2. *-B: The percentage

3. gsort year industry -sale

4. bysort year industry: gen sale123 = _n //排序

5. list in 1/20, sepby(year industry) noobs

6. *-行业总收入

7. bysort year industry: egen sale_ind = total(sale)

8. *-Top 4

9. bysort year industry: egen sale_top4 = total(sale) if sale123<=4

10. *-Ratio

11. gen ratio4_temp = sale_top4/sale_ind

12. *-expand

13. bysort year industry: egen ratio4 = min(ratio4_temp)

最终的数据长这样：

1. . list in 1/20, sepby(year industry) noobs

2. +-------------------------------------------------------------------------------------+

3. | id industry year sale sale123 sale_ind sale_t~4 ratio4_~p ratio4 |

4. |-------------------------------------------------------------------------------------|

5. | 5 1 2011 130 1 1016 498 .49015748 .49015748 |

6. | 7 1 2011 125 2 1016 498 .49015748 .49015748 |

7. | 3 1 2011 124 3 1016 498 .49015748 .49015748 |

8. | 9 1 2011 119 4 1016 498 .49015748 .49015748 |

9. | 10 1 2011 106 5 1016 . . .49015748 |

10. | 4 1 2011 95 6 1016 . . .49015748 |

11. | 2 1 2011 86 7 1016 . . .49015748 |

12. | 1 1 2011 83 8 1016 . . .49015748 |

13. | 6 1 2011 74 9 1016 . . .49015748 |

14. | 8 1 2011 74 10 1016 . . .49015748 |

15. |-------------------------------------------------------------------------------------|

16. | 18 11 2011 130 1 984 493 .50101626 .50101626 |

17. | 15 11 2011 126 2 984 493 .50101626 .50101626 |

18. | 19 11 2011 124 3 984 493 .50101626 .50101626 |

19. | 17 11 2011 113 4 984 493 .50101626 .50101626 |

20. | 13 11 2011 110 5 984 . . .50101626 |

21. | 16 11 2011 92 6 984 . . .50101626 |

22. | 14 11 2011 91 7 984 . . .50101626 |

23. | 11 11 2011 86 8 984 . . .50101626 |

24. | 12 11 2011 68 9 984 . . .50101626 |

25. | 20 11 2011 44 10 984 . . .50101626 |

26. +-------------------------------------------------------------------------------------+

3. 其他解决方法

当然，还有其他解决方法，或许更为简洁。若有，请在下方留言或发邮件至 [e-mail]。

关于我们

• 【Stata 连享会(公众号：StataChina)】由中山大学连玉君老师团队创办，旨在定期与大家分享 Stata 应用的各种经验和技巧。

• 公众号推文同步发布于 【简书-Stata连享会】 和 【知乎-连玉君Stata专栏】。可以在简书和知乎中搜索关键词Stata或Stata连享会后关注我们。

• 推文中的相关数据和程序，以及 Markdown 格式原文 可以在 【Stata连享会-码云】 中获取。【Stata连享会-码云】 中还放置了诸多 Stata 资源和程序。如 Stata命令导航 ||  stata-fundamentals ||  Propensity-score-matching-in-stata || Stata-Training 等。
  

联系我们

• 欢迎赐稿： 欢迎将您的文章或笔记投稿至Stata连享会(公众号: StataChina)，我们会保留您的署名；录用稿件达五篇以上，即可免费获得 Stata 现场培训 (初级或高级选其一) 资格。

• 意见和资料： 欢迎您的宝贵意见，您也可以来信索取推文中提及的程序和数据。

• 招募英才： 欢迎加入我们的团队，一起学习 Stata。合作编辑或撰写稿件五篇以上，即可免费获得 Stata 现场培训 (初级或高级选其一) 资格。

• 联系邮件： [e-mail]

往期精彩推文

---

欢迎加入Stata连享会(公众号: StataChina)

空间计量专题-海报.jpg