如何用stata画下面这样的图：
      数据集中有变量ROA，盈余管理程度DA
第一步：先按ROA的取值范围，设置30个子区间，画出ROA在各个子区间的频数分布直方图
第二步：我想在上面的直方图的基础上，添加一个折线图，表示各个ROA区间上，盈余管理程度DA的均值（中位数），




希望各位大神指点，怎么在stata中实现。感激！感激！

附：我希望达到类似下图的效果：


图片来源：吴联生、薄仙慧、王亚平，2007：《避免亏损的盈余管理程度：上市公司与非上市公司的比较》，《会计研究》第2期。