如何绘制散点图和箱线图结合的图。

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

原先之问题在此 https://bbs.pinggu.org/thread-6252218-1-1.html，而我上美国 Stata forum 寻求解答之答案于此 https://www.statalist.org/forums ... -in-the-same-figure，所有 credit 归于回答者。

其之 code 为：

1. // ericmelse
   
2. * Load example file:
   
3. sysuse nlsw88, clear
   
4. 
   
5. * The distribution underlying the box plot for inspection:
   
6. histogram wage, width(0.1) start(0) lc(navy%50) horizontal yla(, ang(h))
   
7. histogram hours, width(0.1) start(0) lc(navy%50) vertical
   
8. 
   
9. * Graph the box plots for inspection:
   
10. graph box wage, marker(1, ms(oh) mc(navy%50)) yla(, ang(h))
    
11. graph box hours, marker(1, ms(oh) mc(navy%50)) yla(, ang(h)) hor
    
12. 
    
13. * Compute required data for Boxplot of Wage:
    
14. gen maxX = 84    // boxplot x-scale position
    
15. * The median and quartiles are easiest:
    
16. egen med_w = median(wage) //
    
17. egen upq_w = pctile(wage), p(75) //
    
18. egen loq_w = pctile(wage), p(25) //
    
19. egen mean_w = mean(wage) //
    
20. * We could now get the IQR by subtraction, upq - loq, which would be more efficient,
    
21. * but we will mention that it has its own egen function.
    
22. egen iqr_w = iqr(wage) //,
    
23. * For upper limits of whiskers:
    
24. egen upper_w = max(min(wage, upq_w + 1.5 * iqr_w)) //,
    
25. * For lower limits of whiskers:
    
26. egen lower_w = min(max(wage, loq_w - 1.5 * iqr_w)) //,
    
27. 
    
28. * Compute required data for Boxplot of Wage:
    
29. gen maxY = 43 // boxplot y-scale position
    
30. * The median and quartiles are easiest:
    
31. egen med_h = median(hours) //
    
32. egen upq_h = pctile(hours), p(75) //
    
33. egen loq_h = pctile(hours), p(25) //
    
34. egen mean_h = mean(hours) //
    
35. * We could now get the IQR by subtraction, upq - loq, which would be more efficient,
    
36. * but we will mention that it has its own egen function.
    
37. egen iqr_h = iqr(hours) //,
    
38. * For upper limits of whiskers:
    
39. egen upper_h = max(min(hours, upq_h + 1.5 * iqr_h)) //,
    
40. * For lower limits of whiskers:
    
41. egen lower_h = min(max(hours, loq_h - 1.5 * iqr_h)) //,
    
42. 
    
43. * Boxplots and scatterplot in the same figure:
    
44. twoway /// Box plot of Wage
    
45. (rbar med_w upq_w maxX, blc(navy) bfc(white) barw(1.7) ) ///
    
46. (rbar med_w loq_w maxX, blc(navy) bfc(white) barw(1.7) ) ///
    
47. (rspike upq_w upper_w maxX , blc(navy) ) ///
    
48. (rspike loq_w lower_w maxX , blc(navy) ) ///
    
49. (rcap upper_w upper_w maxX , blc(navy) msize(*1) ) ///
    
50. (rcap lower_w lower_w maxX , blc(navy) msize(*1) ) ///
    
51. (scatter mean_w maxX, pstyle(p1) ms(dh) msize(*1) mls(*.8) mlc(gs9) ) /// Diamond for mean
    
52. (scatter wage maxX if !inrange(wage, lower_w, upper_w), ms(oh) mc(navy%20) legend(off)) /// Outliers
    
53. (rbar med_h upq_h maxY, blc(navy) bfc(white) barw(1.3) hor) /// Box plot of Hours
    
54. (rbar med_h loq_h maxY, blc(navy) bfc(white) barw(1.3) hor) ///
    
55. (rspike upq_h upper_h maxY , blc(navy) hor) ///
    
56. (rspike loq_h lower_h maxY , blc(navy) hor) ///
    
57. (rcap upper_h upper_h maxY , blc(navy) msize(*1) hor) ///
    
58. (rcap lower_h lower_h maxY , blc(navy) msize(*1) hor) ///
    
59. (scatter maxY mean_h , pstyle(p1) ms(dh) msize(*1) mls(*.8) mlc(gs9) ) /// Diamond for mean
    
60. (scatter maxY hours if !inrange(hours, lower_h, upper_h), ms(oh) mc(navy%20) legend(off)) /// Outliers
    
61. (scatter wage hours, msymbol(oh) /// Scatter plot of Wage & Hours
    
62. legend(off) xscale(range(-1 80)) ///
    
63. mcolor(navy%50)) ///
    
64. , ytitle(Hourly wage) /// General plot controls
    
65. xtitle(Usual hours worked)
    
66. 
    
67. * The above final code can be edited to modify the plot using the regular twoway options etc.

复制代码

扫码加我 拉你入群

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

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

分享0 收藏16 回帖

关键词：散点图 箱线图 distribution scatterplot INSPECTION

谢谢老师分享！

谢谢老师分享！

谢谢老师分享！

另外也有彩色的 boxplots 如下：
其程序如下（接连上面的）：

1. // color: red/blue
   
2. * Boxplots and scatterplot in the same figure, colors adjusted
   
3. twoway /// Box plot of Wage (vertical)
   
4. (rbar med_w upq_w maxX, blc(blue) bfc(white) barw(3) ) ///
   
5. (rbar med_w loq_w maxX, blc(blue) bfc(white) barw(3) ) ///
   
6. (rspike upq_w upper_w maxX , blc(blue) ) ///
   
7. (rspike loq_w lower_w maxX , blc(blue) ) ///
   
8. (rcap upper_w upper_w maxX , blc(blue) msize(*1) ) ///
   
9. (rcap lower_w lower_w maxX , blc(blue) msize(*1) ) ///
   
10. (scatter mean_w maxX, pstyle(p1) ms(dh) msize(*1) mls(*.8) mlc(navy%20) ) /// Diamond for mean (horizontal)
    
11. (scatter wage maxX if !inrange(wage, lower_w, upper_w), ms(oh) mc(blue%20) legend(off)) /// Outliers
    
12. (rbar med_h upq_h maxY, blc(red) bfc(white) barw(1.5) hor) /// Box plot of Hours
    
13. (rbar med_h loq_h maxY, blc(red) bfc(white) barw(1.5) hor) ///
    
14. (rspike upq_h upper_h maxY , blc(red) hor) ///
    
15. (rspike loq_h lower_h maxY , blc(red) hor) ///
    
16. (rcap upper_h upper_h maxY , blc(red) msize(*1) hor) ///
    
17. (rcap lower_h lower_h maxY , blc(red) msize(*1) hor) ///
    
18. (scatter maxY mean_h , pstyle(p1) ms(dh) msize(*1) mls(*.8) mlc(red%20) ) /// Diamond for mean
    
19. (scatter maxY hours if !inrange(hours, lower_h, upper_h), ms(oh) mc(red%20) legend(off)) /// Outliers
    
20. (scatter wage hours, msymbol(oh) /// Scatter plot of Wage & Hours
    
21.            legend(off) xscale(range(-1 80)) ///
    
22.            mcolor(navy%50)) ///
    
23. , xsize(9) ysize(9) graphregion(margin(1 -3 1 -3)) /// General plot controls from here
    
24. yscale(noextend) xscale(noextend) /// Disable axis extension
    
25. ytitle(Hourly wage, margin(0 2 0 0)) /// specified margins; left, right, bottom, top
    
26. xtitle(Usual hours worked, margin(0 0 0 2)) // last code line

复制代码

这个程序有用，多谢黄老师分享

yuqr1986 发表于 2018-3-5 09:04
这个程序有用，多谢黄老师分享

老实说，还蛮 fancy 的。

最后画图的程序看不懂头都大了
老师厉害！

阿璇2016 发表于 2018-3-5 16:34
最后画图的程序看不懂头都大了
老师厉害！

这是别人写的，不是我，我也不是很懂（还没花时间去了解）。借花献佛而已。

厉害厉害，黄老师出品 必为精品