[Stata] 如何添加直角坐标系及极坐标系及SDAS_绘图实例

1. sysuse sp500,clear
   
2. twoway parea high date, pattern(pattern8) || ///
   
3.         parea low date, pattern(pattern6) scale(0.5)
   
4. 
   
5. sysuse auto,clear
   
6.         twoway parea l w, sort pattern(pattern6)
   
7. 
   
8. twoway parea d w, sort pattern(pattern9) ||  ///
   
9.         parea l w if w<=3600, sort pattern(pattern6) || ///
   
10.         parea l w if w>=3600, sort pattern(pattern4) ,legend(rows(1))

复制代码

1. /*
   
2. 
   
3.    Sergiy Radyakin 12. Aug, 2007
   
4. 
   
5. Demonstration of -parea- capabilities.
   
6. Draw 10 patterned rectangles in a 2 by 5 formation
   
7. 
   
8. */
   
9. 
   
10. #delimit ;
    
11. 
    
12. clear;
    
13. set obs 60;
    
14. gen x=. ;
    
15. gen y=. ;
    
16. gen p=. ;
    
17. 
    
18. forvalues i=1/10 { ;
    
19.   local sh_y=0 ;
    
20.   if `i'>5 local sh_y=1 ;
    
21.   quietly replace x=`i'-.2-`sh_y'*5  in `=`i'*6-5' ;
    
22.   quietly replace x=`i'-1 -`sh_y'*5  in `=`i'*6-4' ;
    
23.   quietly replace x=`i'-1 -`sh_y'*5  in `=`i'*6-3' ;
    
24.   quietly replace x=`i'-.2-`sh_y'*5  in `=`i'*6-2' ;
    
25.   quietly replace x=`i'-.2-`sh_y'*5  in `=`i'*6-1' ;
    
26. 
    
27.   quietly replace y=`sh_y'*0.55      in `=`i'*6-5' ;
    
28.   quietly replace y=`sh_y'*0.55      in `=`i'*6-4' ;
    
29.   quietly replace y=`sh_y'*0.55+.5   in `=`i'*6-3' ;
    
30.   quietly replace y=`sh_y'*0.55+.5   in `=`i'*6-2' ;
    
31.   quietly replace y=`sh_y'*0.55      in `=`i'*6-1' ;  
    
32. 
    
33.   quietly replace p=`i' in `=`i'*6-5'/`=`i'*6-1'   ;
    
34.   local labl `"`labl' label(`i' "pattern`i'")"'    ;
    
35. } ;
    
36. 
    
37. local GRAPHCMD="twoway " ;
    
38. forvalues i=1/10 { ;
    
39.   local GRAPHCMD `"`GRAPHCMD' parea y x if p==`i',lc(black) fc(black)
    
40.                            pattern(pattern`i') fi(100) nodropb `=cond(`i'<10,"||","")'"' ;
    
41. };
    
42. 
    
43. local GRAPHCMD `"`GRAPHCMD' graphregion(color(white) margin(zero))
    
44.                             plotregion(style(none) color(white))
    
45.                             xlabel(none) xscale(off r(0 4.5))
    
46.                             ylabel(none) yscale(off r(0 1.05))
    
47.                             legend(cols(5) `labl') scale(0.8)
    
48.                             title("twoway parea y x,...") subtitle("Sergiy Radyakin") "'  ;
    
49. 
    
50. `GRAPHCMD' ;

复制代码

好贴，收藏了，感谢楼主的辛勤劳动！