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

1. // Displaying part of a histogram
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. twoway__histogram_gen mpg , gen(h y, replace) freq bin(10)
   
6. 
   
7. twoway ///
   
8. (histogram y [fw=h] if inrange(y,19,35),discrete color(green)lcolor(red)freq) ///
   
9. (histogram y [fw=h],color(none) lcolor(red) discrete freq name(a1, replace))
   
10. 
    
11. twoway ///
    
12. (histogram y [fw=h] if inrange(y,19,35),discrete color(green) lcolor(red)freq ///
    
13. name(a2, replace))
    
14. 
    
15. graph combine a1 a2 ,xcommon cols(1)

复制代码

1. // Changing bar labels
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. graph bar  mpg  gear, showyvars ///
   
6. yvaroptions(relabel(1 "MPG" 2"GEAR RATIO")) legend(off)
   
7. 
   
8. exit

复制代码

1. // Changing bar labels
   
2. //It is sometimes more convenient to use value lables rather than a graph
   
3. //label options to change graph bar labels. In the example it also allows the
   
4. //legend to be spread over the width of the graph.
   
5. 
   
6. clear all
   
7. sysuse auto
   
8. 
   
9. lab def origin 0 "Europe de l`=char(146)'Ouest" ///
   
10. 1 "Asie de l`=char(146)'Est", modify
    
11. 
    
12. graph hbar mpg trunk turn, over(foreign) ///
    
13. legend(row(1) span) stack name(two,replace)
    
14. 
    
15. exit
    
16. 
    
17. *( Stata list,  25 Jun 2011)

复制代码

1. // Different color bars for the over option - Using twoway graph
   
2. clear all
   
3. 
   
4. input ///
   
5. id w pos mark
   
6. 1  0 6  84.58
   
7. 2  0 7  84.51
   
8. 3  0 8  84.12
   
9. 4  0 9  83.34
   
10. 5  0 10 82.8
    
11. 6  1 1  69.55
    
12. 7  1 2  65.16
    
13. 8  1 3  64.91
    
14. 9  1 4  64.53
    
15. 10 1 5  63.70
    
16. end
    
17. 
    
18. label define w 1 "Women" 0 "Men"
    
19. label value w w
    
20. 
    
21. label define kk        ///
    
22. 1  "Keshorn - TRI "    ///
    
23. 2  "Oleksandr - UKR "  ///
    
24. 3  "Antti - FIN "      ///
    
25. 4  "Vitezslav - CZE"   ///
    
26. 5  "Tero - FIN"        ///
    
27. 6  " "                 ///
    
28. 7  "Barbora - CZE "    ///
    
29. 8  "Christina - GER "  ///
    
30. 9  "Linda - GER "      ///
    
31. 10  "Sunette - RSA"    ///
    
32. 11 "Huihui - CHN"      
    
33. 
    
34. label value id  kk
    
35. 
    
36. replace id=id+1 if id>5
    
37. 
    
38. set obs `=_N+1'
    
39. replace id=6 in l
    
40. 
    
41. twoway (bar mark id if id==1, base(0) ylabel(0(10)90) ///
    
42. xtitle("Men                                  Women")  ///
    
43. xlabel(1(1)11,valuelabel angle(45))                   ///
    
44. barwidth(.4)   color(cyan))                           ///
    
45. (bar mark id if id==2,  barwidth(.4)   color(cyan))   ///
    
46. (bar mark id if id==3,  barwidth(.4)   color(cyan))   ///
    
47. (bar mark id if id==4,  barwidth(.4)   color(cyan))   ///
    
48. (bar mark id if id==5,  barwidth(.4)   color(cyan))   ///
    
49. (bar mark id if id==7,  barwidth(.4)   color(red))    ///
    
50. (bar mark id if id==8,  barwidth(.4)   color(red))    ///
    
51. (bar mark id if id==9,  barwidth(.4)   color(red))    ///
    
52. (bar mark id if id==10, barwidth(.4)   color(red))    ///
    
53. (bar mark id if id==11, barwidth(.4)   color(red))    ///
    
54. (scatter mark id,  mlabel(mark) msymbol(none)         ///
    
55. mlabposition(12)), legend(off)                        ///
    
56. title("London Olympics 2012" "Javelin")               ///
    
57. plotregion(margin(l=5 b=0))
    
58. 
    
59. exit

复制代码

1. // Difference from Average by levels of groups
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. graph drop _all
   
6. 
   
7. bysort for rep78 : egen mpg1=mean(mpg)
   
8. bysort  rep78: egen mpg2=mean(mpg)
   
9. 
   
10. gen mpg3=mpg1-mpg2
    
11. 
    
12. replace mpg2=0 if rep78==.
    
13. 
    
14. twoway ///
    
15. (bar mpg3 rep78 if for==0 & rep78==1 ,fcolor(green) ///
    
16. lcolor(green) barwidth(.5)horiz  yscale(noline))  ///
    
17. (bar  mpg3 rep78 if for==0 & rep78==2 ,fcolor(blue) ///
    
18. lcolor(blue) barwidth(.5)horiz  yscale(noline))   ///
    
19. (bar  mpg3 rep78 if for==0 & rep78==3 ,fcolor(red)  ///
    
20. lcolor(red) barwidth(.5)horiz  yscale(noline))    ///
    
21. (bar  mpg3 rep78 if for==0 & rep78==4 ,fcolor(pink) ///
    
22. lcolor(pink) barwidth(.5)horiz  yscale(noline))   ///
    
23. (bar  mpg3 rep78 if for==0 & rep78==5 ,fcolor(cyan) ///
    
24. lcolor(cyan) barwidth(.5)horiz  yscale(noline)) , ///
    
25. xtitle(MPG) xline(0)  ///
    
26. title("Foreign made cars") legend(off) name(k1, replace)
    
27. 
    
28. twoway ///
    
29. (bar  mpg3 rep78 if for==1 & rep78==1 ,fcolor(green) ///
    
30. lcolor(green) barwidth(.5)horiz  yscale(noline)) ///
    
31. (bar  mpg3 rep78 if for==1 & rep78==2 ,fcolor(blue) ///
    
32. lcolor(blue) barwidth(.5)horiz  yscale(noline))   ///
    
33. (bar  mpg3 rep78 if for==1 & rep78==3 ,fcolor(red)  ///
    
34. lcolor(red) barwidth(.5)horiz  yscale(noline))    ///
    
35. (bar  mpg3 rep78 if for==1 & rep78==4 ,fcolor(pink) ///
    
36. lcolor(pink) barwidth(.5)horiz  yscale(noline))    ///
    
37. (bar  mpg3 rep78 if for==1 & rep78==5 ,fcolor(cyan) ///
    
38. lcolor(cyan) barwidth(.5)horiz  yscale(noline)) , ///
    
39. xtitle(MPG) xline(0) yscale(off) ///
    
40. title("Domestically made cars") legend(off) name(k2, replace)
    
41. 
    
42. graph combine k1 k2, row(1)  imargin(10 0 0 0) ycommon ///
    
43. title("Difference from average of mpg for " ///
    
44. "Foreign or Domestically produced cars" " " )

复制代码

1. *修改后，添加第一行
   
2. 
   
3. sysuse lifeexp.dta, clear
   
4. replace lexp = 35 if country == "Haiti"
   
5. egen median = median(lexp), by(region)
   
6. egen upq = pctile(lexp), p(75) by(region)
   
7. egen loq = pctile(lexp), p(25) by(region)
   
8. egen iqr = iqr(lexp), by(region)
   
9. egen upper = max(min(lexp, upq + 1.5 * iqr)), by(region)
   
10. egen lower = min(max(lexp, loq - 1.5 * iqr)), by(region)
    
11. twoway rbar med upq region, ///
    
12. pstyle(p1) blc(gs15) bfc(gs8) barw(0.35) || ///
    
13. rbar med loq region, pstyle(p1) blc(gs15) bfc(gs8) barw(0.35) || ///
    
14. rspike upq upper region, pstyle(p1) || ///
    
15. rspike loq lower region, pstyle(p1) || ///
    
16. rcap upper upper region, pstyle(p1) msize(*2) || ///
    
17. rcap lower lower region, pstyle(p1) msize(*2) || ///
    
18. scatter lexp region if !inrange(lexp, lower, upper) & lexp > 50, ///
    
19. ms(Oh)  mla(country) mlabcolor(gs8) xscale(off)  legend(off) ///
    
20. yla(, ang(h)) ytitle(Life expectancy (years)) xtitle("") ///
    
21. name(gr1,replace)
    
22. 
    
23. scatter lexp region if !inrange(lexp, lower, upper) & lexp < 50, ///
    
24. ylabel(35)   fysize(20) ytitle("") ms(Oh) mlabcolor(blue) ///
    
25. mcolor(blue)  mla(country)  xlabel(1 `" "Europe and" "Central Asia" "' ///
    
26. 2 "North America" 3 "South America", noticks) ///
    
27. ytitle(Life expectancy (years), color(white))  yla(, ang(h))  ///
    
28.   xtitle("")  name(gr2,replace)
    
29. 
    
30. graph combine gr1 gr2, col(1) xcommon imargin(zero)

复制代码

1. // Combining graphs and including a common legend
   
2. //this uses the user written command "grc1leg" which must be downloaded first
   
3. 
   
4. net install grc1leg,from( http://www.stata.com/users/vwiggins/)
   
5. 
   
6. clear        all
   
7. 
   
8. input quest str25 q      a1 a2 a3 a4 a5 a6
   
9. 1 "Question 1"           0  2  37 45 12 4
   
10. 1 "Benchmark Q1"         2  5  25 47 17 4
    
11. 2 "Question 2"           1  37  2 40 17 3
    
12. 2 "Benchmark Q2"         2  5  25 47  4 17
    
13. 3 "Question 3"           1  2  40 37 17 3
    
14. 3 "Benchmark Q3"         2  5  25 47 17 4
    
15. 4 "Question 4"           1  2  37  17 3 40
    
16. 4 "Benchmark Q4"         2  5  47 25 17 4
    
17. end       
    
18. 
    
19. graph hbar a1-a6 if quest==1, percent over(quest, gap(1)) over(q, gap(10))   ///
    
20. stack legend(off) yscale(off) yline(20 40 60 80 ,lwidth(0.25) lcolor(black)) ///
    
21. saving(a1, replace)yscale(off) plotregion(margin( b+3 t+3))  outergap(30)
    
22. 
    
23. graph hbar a1-a6 if quest==2, percent over(quest, gap(1)) over(q, gap(10))   ///
    
24. stack legend(off) yscale(off) yline(20 40 60 80 ,lwidth(0.25)lcolor(black))  ///
    
25. saving(a1a, replace)yscale(off) plotregion(margin( b+3 t+3))  outergap(30)
    
26. 
    
27. graph hbar a1-a6 if quest==3, percent over(quest, gap(1)) over(q, gap(10))    ///
    
28. stack legend(off) yscale(off) yline(20 40 60 80 ,lwidth(0.25)  lcolor(black)) ///
    
29. saving(a1b, replace)yscale(off) plotregion(margin( b+3 t+3))  outergap(30)
    
30. 
    
31. graph hbar a1-a6 if quest==4, percent over(quest, gap(1)) over(q, gap(20))    ///
    
32. stack legend(span rows(1) label(1 "Missing") label(2 "Never")                 ///
    
33. label(3 "Rarely") label(4 "Occasionly  ") label(5 "Mostly")                   ///
    
34. label(6 "Always    ") size(small))                                            ///
    
35. yline(20 40 60 80 ,lwidth(0.25 ) lcolor(black)) saving(a2, replace)
    
36. 
    
37. grc1leg a1.gph a1a.gph a1b.gph a2.gph, cols(1) ///
    
38. imargin(0 0 0 0) ycommon xcommon legendfrom(a2.gph)
    
39. 
    
40. exit
    
41. 
    
42. *(Also see the user written program catplot on:
    
43. http://www.survey-design.com.au/Usergraphs.html)

复制代码

1. // histogram - add x bar labels for the centre of each bar
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. twoway__histogram_gen mpg  , bin(6) generate(h x , replace )
   
6. 
   
7. count if !missing(x)
   
8. 
   
9. forvalues i=1/`=r(N)' {
   
10. 
    
11.   local label `label'  `=string(x[`i'],"%8.2f")'
    
12. 
    
13. }
    
14. 
    
15. histogram mpg , bin(6) xlabel(`label')

复制代码

1. // Changing bar labels using the Stata graph editor
   
2. 
   
3. // The type of bar label is limited in Stata graphs but sometimes
   
4. // additional information is required. This can be added with the
   
5. // Stata graphics editor commands.
   
6. // Below the number of observation for each region is added.
   
7. 
   
8. sysuse citytemp, clear
   
9. 
   
10. #delimit ;
    
11. graph bar tempjuly tempjan, over(region) bargap(-30)
    
12.   legend( label(1 "July") label(2 "January") )
    
13.   ytitle("Degrees Fahrenheit")
    
14.   title("Average July and January temperatures")
    
15.   subtitle("by regions of the United States")
    
16.   note("Source:  U.S. Census Bureau, U.S. Dept. of Commerce")
    
17. blabel(bar, position(inside) format(%9.1f) color(white)) ;
    
18. #delimit cr
    
19. 
    
20. generate No=_n
    
21. collapse (sum) N, by(region)
    
22. list
    
23. 
    
24. local i1=1
    
25. 
    
26. forvalues i=1/`=_N' {
    
27. 
    
28. gr_edit plotregion1.barlabels[`i1'].text = {}
    
29. gr_edit plotregion1.barlabels[`i1'].text.Arrpush N=`=No[`i']'
    
30. local ++i1
    
31. 
    
32. gr_edit plotregion1.barlabels[`i1'].text = {}
    
33. gr_edit plotregion1.barlabels[`i1'].text.Arrpush N=`=No[`i']'
    
34. local ++i1
    
35. 
    
36. }
    
37. 
    
38. exit

复制代码

1. // Adding bar labels
   
2. 
   
3. // The type of bar label is limited in Stata graphs but sometimes
   
4. // additional information is required. This can be added with the
   
5. // by adding another twoway plot.
   
6. 
   
7. sysuse census, clear
   
8. graph bar (sum) marri div, over(region) stack
   
9. 
   
10. more
    
11. 
    
12. collapse (sum) marri div, by(region)
    
13. gen ratio = div / (marri + div)
    
14. gen sratio = string(ratio, "%4.3f")
    
15. gen marriPdiv = marria + div
    
16. 
    
17. twoway bar marria region, barw(0.6) base(0) ///
    
18. || rbar marria marriPdiv region, barw(0.6) xla(1/4, noticks valuelabel) ///
    
19. || scatter marriPdiv region, ms(none) mla(sratio) mlabpos(12)  ///
    
20. legend(order(1 "marriages" 2 "divorces")) ///
    
21. yla(0 4e5 "400"  8e5 "800" 12e5 "1200" , ang(h)) ytitle(thousands) ///
    
22. ysc(r(. 14e5))
    
23. 
    
24. *(Nick Cox - Stata list,  Thu, 17 Feb 2011)

复制代码

---