Histogram + normal 图

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

刚刚看看 Joseph Coveney 之建议的图，还蛮有趣的，与大家方享： 其程序为：

1. clear *
   
2. 
   
3. sysuse auto, clear
   
4. quietly summarize mpg if foreign
   
5. local var1mean: display %6.2f r(mean)
   
6. local var1sd: display %6.2f r(sd)
   
7. local var1min: display %6.2f r(min)
   
8. local var1max: display %6.2f r(max)
   
9. quietly summarize mpg if !foreign
   
10. local var2mean : display %6.2f r(mean)
    
11. local var2sd: display %6.2f r(sd)
    
12. local var2min: display %6.2f r(min)
    
13. local var2max: display %6.2f r(max)
    
14. 
    
15. graph twoway ///
    
16.   (histogram mpg if !foreign, start(10) width(2) frequency bfcolor(brown%20) blcolor(brown)) ///
    
17.   (histogram mpg if foreign, frequency start(10) width(2) barw(1.8) bfcolor(navy%20) blcolor(navy)) ///
    
18.   (function y=normalden(x, `var1mean', `var1sd')  * `var1max' , range(`var1min' `var1max') lcolor(navy)) ///
    
19.   (function y=normalden(x, `var2mean', `var2sd') * `var2max' , range(`var2min' `var2max') lcolor(brown)) ///
    
20.    , xlabel(10(5)45, labsize(small)) ylabel(, angle(horizontal) labsize(vsmall) format(%3.0f)) ///
    
21.    legend(order(1 "Domestic" 2 "Foreign") pos(2) ring(0) col(1)) ///
    
22.    note(" " "MPG Domestic:  {it:M} = `var2mean'  {it:SD} = `var2sd'" ///
    
23.    "MPG Foreign:     {it:M} = `var1mean'  {it:SD} = `var1sd'", size(vsmall)) ///
    
24.    scheme(s1color) name(HistPlusNormal1)
    
25. 
    
26. graph twoway ///
    
27.   (histogram mpg if !foreign, start(10) width(2) fraction bfcolor(brown%20) blcolor(brown)) ///
    
28.   (histogram mpg if foreign, fraction start(10) width(2) barw(1.8) bfcolor(navy%20) blcolor(navy)) ///
    
29.   (function y=normalden(x, `var1mean', `var1sd') , range(`var1min' `var1max') lcolor(navy)) ///
    
30.   (function y=normalden(x, `var2mean', `var2sd') , range(`var2min' `var2max') lcolor(brown)) ///
    
31.    , xlabel(10(5)45, labsize(small)) ylabel(, angle(horizontal) labsize(vsmall) format(%3.2f)) ///
    
32.    legend(order(1 "Domestic" 2 "Foreign") pos(2) ring(0) col(1)) ///
    
33.    note(" " "MPG Domestic:  {it:M} = `var2mean'  {it:SD} = `var2sd'" ///
    
34.    "MPG Foreign:     {it:M} = `var1mean'  {it:SD} = `var1sd'", size(vsmall)) ///
    
35.    scheme(s1color) name(HistPlusNormal2)
    
36.   
    
37. graph combine HistPlusNormal1 HistPlusNormal2, xcommon

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏9 回帖

关键词：histogram normal gram Norm hist

支持下黄老师传的新命令

不是看了这个code，再查了twoway hist的说明，还真不知道原来hist之所以能附加一个正态分布图是通过twoway hist才实现的

谢谢黄老师上传的新命令

静心品茶 发表于 2018-7-10 14:19
谢谢黄老师上传的新命令

我也喜欢这个图！

黃河泉 发表于 2018-7-1 16:03
刚刚看看 Joseph Coveney 之建议的图，还蛮有趣的，与大家方享：其程序为：

谢谢黄老师！

谢谢黄老师！