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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

本贴转自：http://www.surveydesign.com.au/tipsgraphs.html _该例讲述如何在图像中添加直角坐标系 极坐标系须先安装polar命令

扫码加我 拉你入群

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

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

分享0 收藏9 回帖

关键词：Stata tata 坐标系 Das SDA 平面 坐标系 如何 black local

谢谢分享~~~

1. sysuse nlsw88, clear       
   
2. 
   
3. generate g=max(grade,8)
   
4. label  var g "Years of education"       
   
5. rename  married m       
   
6. collapse wage, by(g m)       
   
7. generate  w=round(wage,.01)       
   
8. summarize wage, meanonly       
   
9. generate  c=round((wage-r(min))/(r(max)-r(min))*255)       
   
10. quiet levelsof c, loc(cs)       
    
11. local g               
    
12. foreach c of loc cs {       
    
13.    local c1=120+round(`c'/2)
    
14.    local c2=255-round(`c'/2)       
    
15.    local m mc("`c2' `c2' `c1'")       
    
16.    local g `g'||sc m g if c==`c', ms(S) msize(ehuge) `m'
    
17. }       
    
18. local g `g'||sc m g, ms(i) mlab(w) mlabp(0)       
    
19. local g `g' leg(off) yla(-.75 " " 0 "N" 1 "Y" 1.75 " ", notick)            
    
20. local g `g' xla(7 " " 8/18 19 " ", notick)               
    
21. 
    
22. twoway `g'  scheme(s1mono) ti(Wage heat map) xsize(8) ysize(5)

复制代码

1. clear               
   
2.                
   
3.   input input low   high  mean       
   
4.          1     1.6   2.2  1.9       
   
5.          2     1.5   2.3  1.9
   
6.          5     2.9   0.01 1.9
   
7.                  3     1.0   2.9  1.9
   
8.          4     0.5   3.2  1.9       
   
9.   end                                               
   
10.                                                                                
    
11.   label var input "Input"       
    
12.   label var low   "Low value"       
    
13.   label var high  "High value"       
    
14.   label var mean  "Mean"       
    
15.                                                                                
    
16.   label define inputlbl             ///       
    
17.   1 "Annual Donations in U.S."      ///
    
18.   2 "Seroprevalence"                ///
    
19.   3 "Exposures per Donation"            ///
    
20.   4 "Probability of Transmission"   ///       
    
21.   5 "Sensitivity of Test"                         
    
22.                                  
    
23.   label val input inputlbl       
    
24.                                                        
    
25.   gen sens = high-low       
    
26.   sort sens
    
27.                                                                                
    
28.   local ml1 = (2.2+1.9)/2       
    
29.   local mh1 = (1.6+1.9)/2       
    
30.   local ml2 = (2.3+1.9)/2       
    
31.   local mh2 = (1.5+1.9)/2       
    
32.   local ml3 = (2.9+1.9)/2       
    
33.   local mh3 = (1.0+1.9)/2       
    
34.   local ml4 = (3.2+1.9)/2       
    
35.   local mh4 = (0.5+1.9)/2       
    
36.   local ml5 = (2.9+1.9)/2       
    
37.   local mh5 = (0  +1.9)/2       
    
38.                
    
39.                
    
40. local format fcolor(white) lcolor(white)
    
41.                                                                                
    
42.   twoway                                                                   ///
    
43.   (rbar low mean input, horizontal)                                        ///
    
44.   (rbar mean high input, horizontal),                                      ///
    
45.   ylabel(, val angle(horizontal)) xsize(7.1) ysize(5.2) scheme(sj)         ///
    
46.   xtitle(" " "Predicted Cases of Transfusion-Transmitted T.   Cruzi")      ///
    
47.   title("Sensitivity of Results to Inputs") subtitle("Test Once Scenario") ///
    
48.   legend(label(1 "Low Input Values") label(2 "High Input Values") rows(2)) ///
    
49.   text(5 `ml5' "93.5%", box j(center)  margin(l+1 r+1 t+1 b+1)   `format') ///
    
50.   text(5 `mh5' "99.99%", box j(center) margin(l+1 r+1 t+1 b+1)   `format') ///
    
51.   text(4 `ml4' "5%", box j(center)     margin(l+1 r+1 t+1 b+1)   `format') ///
    
52.   text(4 `mh4' "0.8%", box j(center)   margin(l+1 r+1 t+1 b+1)   `format') ///
    
53.   text(3 `ml3' "3", box j(center)      margin(l+1 r+1 t+1 b+1)   `format') ///
    
54.   text(3 `mh3' "1", box j(center)      margin(l+1 r+1 t+1 b+1)   `format') ///
    
55.   text(2 1.4 "1:28,600", box           margin(l+1 r+1 t+1 b+1)   `format'  ///
    
56.   placement(9) j(right))                                                   ///
    
57.   text(2 2.4 "1:18,200", box           margin(l+1 r+1 t+1 b+1)   `format'  ///
    
58.   placement(3) j(left))                                                    ///
    
59.   text(1 1.5 "13.4 million", box       margin(l+1 r+1 t+1 b+1)   `format'  ///
    
60.   placement(9) j(right))                                                   ///
    
61.   text(1 2.4 "19 million", box         margin(l+1 r+1 t+1 b+1)   `format'  ///
    
62.   placement(3) j(left))       

复制代码

1. // Controlling the side of the graph that the axis is on
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. twoway ///
   
6. (histogram mpg, width(5) yscale(alt  axis(1)) ) ///
   
7. (line  weight mpg, yaxis(2) yscale(alt axis(2)) sort)

复制代码

1. // Grid on top of graph
   
2. 
   
3. ssc install freduse, replace  //gets program online
   
4. freduse MPRIME, clear
   
5. generate ym = mofd(daten)
   
6. tsset ym, monthly             //sets for time series
   
7. twoway                                                                    ///
   
8. function y=20.705,range(119 130) recast(area) color(gs12) base(4.7025) || ///
   
9. function y=20.705,range(166 182) recast(area) color(gs12) base(4.7025) || ///
   
10. function y=20.705,range(240 274) recast(area) color(gs12) base(4.7025) || ///
    
11. function y= 5    ,range(119 274) lstyle(grid)                          || ///
    
12. function y=10    ,range(119 274) lstyle(grid)                          || ///
    
13. function y=15    ,range(119 274) lstyle(grid)                          || ///
    
14. function y=20    ,range(119 274) lstyle(grid)                          || ///
    
15. function y=ym(1980,11), range(4.7025 20.705) horizontal lstyle(grid)   || ///
    
16. tsline MPRIME if tin(1970m1,1990m1), xlabel(,format(%tm)) lstyle(p1)      ///
    
17. legend(order(5 1 "Recession")) tlabel(,grid) scheme(s2color)
    
18. 
    
19. *(Maarten Buis - Stata list,  30 October 2011)

复制代码

1. // stacked bar graph
   
2. 
   
3. sysuse auto, clear
   
4. twoway__histogram_gen mpg if foreign, gen(freq_foreign x) ///
   
5. freq start(0) width(2)       
   
6. drop if missing(x)
   
7. keep freq_foreign x
   
8. save a1, replace
   
9. 
   
10. sysuse auto, clear
    
11. twoway__histogram_gen mpg if !foreign, gen(freq_domestic x) ///
    
12. freq start(0) width(2)       
    
13. drop if missing(x)
    
14. keep freq_domestic x mpg
    
15. 
    
16. merge 1:1         x using a1
    
17. 
    
18. mvencode freq_domestic  freq_foreign , mv(.=0)
    
19.        
    
20. egen freq = rowtotal(freq_domestic  freq_foreign )
    
21. 
    
22. twoway                                                                        ///
    
23. (bar freq_domestic x if freq < ., barw(2) bcolor(blue) xsc(r(10 .)) )         ///
    
24. rbar freq_domestic freq x if freq < ., barw(2) bcolor(red)                    ///
    
25. legend(order(1 "Domestic" 2 "Foreign") pos(11) col(1) ring(0))                ///
    
26. xtitle("`: var label mpg'") ytitle(Frequency) xla(10(5)45, alt) yla(, ang(h)) ///
    
27. xsc(r(10 .))

复制代码

1. // Shade Part of a Graph
   
2. 
   
3. sysuse auto, clear
   
4. sort mpg
   
5. generate x=_n
   
6. gen max=ceil(_N/10)*10
   
7. 
   
8. summarize mpg, detail
   
9. 
   
10. twoway ///
    
11. (area max x   if inrange(x, `=r(p1)', `=r(p25)'),   ///
    
12. lcolor(bg) fcolor(green)  fintensity(inten20))      ///
    
13. (area  max x   if inrange(x, `=r(p75)', `=r(p99)'), ///
    
14. lcolor(bg) fcolor(green) fintensity(inten20))       ///
    
15. (line x mpg),legend(off) plotregion(margin(zero))   ///
    
16. xtitle("MPG")                                       ///
    
17. xlabel(10 `=r(p1)' `=r(p25)' `=r(p75)' `=r(p99)' )  ///
    
18. yscale(range(0 `=_N'))

复制代码

1. sysuse auto, clear
   
2. 
   
3. graph bar  mpg, asyvars  over(rep78) ///
   
4. bar(1, bcolor(red*0.4))              ///
   
5. bar(2, bcolor(red*0.2))              ///
   
6. bar(3, bcolor(blue*0.2))             ///
   
7. bar(4, bcolor(blue*0.4))             ///
   
8. bar(5, bcolor(blue*0.6))             ///
   
9. legend(row(1) title("REP78"))        ///
   
10. title("Mean MPG for rep78")

复制代码

1. // box plot
   
2. 
   
3. sysuse auto, clear       
   
4.                                                
   
5. statsby "su mpg" mean=r(mean) max=r(max) min=r(min) sd=r(sd), by(rep78)
   
6.        
   
7. // here k = 1; tune to choice
   
8. gen meanpsd = mean + sd
   
9. gen meanmsd = mean - sd
   
10.                                                
    
11. // use -line- or -connected- not -scatter-
    
12. twoway scatter mean rep78 ||                            ///
    
13. rbar mean max rep78, bcolor(none) barw(0.2) ||          ///
    
14. rbar mean min rep78, bcolor(none) barw(0.2) ||          ///
    
15. rbar meanmsd meanpsd rep78 , barw(0.1) bcolor(red) || , ///
    
16. legend(off) ytitle(Miles per gallon) yla(, ang(h))      ///
    
17. note(red bars: mean +/- sd; empty bars show range)

复制代码