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

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. // 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. // Different color bars for the over option
   
2. clear
   
3. 
   
4. input ///
   
5. id w pos mark
   
6. 1  1 1  69.55
   
7. 2  1 2  65.16
   
8. 3  1 3  64.91
   
9. 4  1 4  64.53
   
10. 5  1 5  63.70
    
11. 6  0 6  84.58
    
12. 7  0 7  84.51
    
13. 8  0 8  84.12
    
14. 9  0 9  83.34
    
15. 10 0 10 82.8
    
16. end
    
17. 
    
18. label define w 1 "Women" 0 "Men"
    
19. label value w w
    
20. 
    
21. gen g1=mark if _n<6
    
22. gen g2=mark if _n>5
    
23. 
    
24. list
    
25. 
    
26. graph bar (asis) g1 g2, ascat over(mark,sort(mark) descending ///
    
27. relabel(  ///
    
28. 5  "Barbora - CZE "    ///
    
29. 4  "Christina - GER "  ///
    
30. 3  "Linda - GER "      ///
    
31. 2  "Sunette - RSA"     ///
    
32. 1  "Huihui - CHN"      ///
    
33. 10 "Keshorn - TRI "    ///
    
34. 9  "Oleksandr - UKR "  ///
    
35. 8  "Antti - FIN "      ///
    
36. 7  "Vitezslav - CZE"   ///
    
37. 6  "Tero - FIN"        ///
    
38. ) label( angle(45)))   over(w) ylabel(0(10)90)                                 ///
    
39. ytitle("Metres") blabel(bar, position(inside) format(%9.1f) color(white))      ///
    
40. scheme(s2color)  legend(off) bargap(5) title("London Olympics 2012" "Javelin") ///
    
41. nofill  bar(1, color(dknavy)) bar(2, color(red))
    
42. 
    
43. exit

复制代码

1. // Adding average bar to group
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. summ price if for==0 & !missing(rep78)
   
6. local ave1=r(mean)
   
7. 
   
8. summ price if for==1 & !missing(rep78)
   
9. local ave2=r(mean)
   
10. 
    
11. set obs `=_N+3'
    
12. 
    
13. replace for=0 in -3
    
14. replace rep78=6 in -3
    
15. 
    
16. replace for=0 in -2
    
17. replace rep78=7 in -2
    
18. replace price=`ave1' in -2
    
19. 
    
20. replace for=1 in l
    
21. replace rep78=7 in l
    
22. replace price=`ave2' in l
    
23. 
    
24. label define kk 1 1 2 2 3 3 4 4 5 5 6 "-"  7 "AVE"
    
25. label value rep78 kk
    
26. 
    
27. graph bar (mean) price , over(rep78 ) over(for) showyvars ///
    
28.   asyvar legend(off)  scheme(color) bar(7, color(black))

复制代码

1. // Detail and Summary box plots together
   
2. 
   
3. sysuse auto, clear
   
4. generate order = _n
   
5. expand 3
   
6. bysort order : generate which = _n
   
7. drop if which == 1 & price > 5000
   
8. drop if which == 2 & price > 10000
   
9. label def which 1 "<= $5000" 2 "<= $10000" 3 "all"
   
10. label val which which
    
11. graph box mpg, over(which) over(foreign)
    
12. 
    
13. *(Nick Cox - Stata list,  December 06, 2013)

复制代码

1. sysuse nlsw88, clear
   
2. 
   
3. clonevar wagelog10=wage
   
4. replace wagelog10=log10(wagelog10)
   
5. 
   
6. mylabels 0(10)40 , myscale(log10(@)) local(labels)
   
7. 
   
8. graph hbox wagelog10, over(ind, sort(1)) nooutside ytitle("") ylabel(`labels') ///
   
9. title("Hourly wage, 1988, woman aged 34-46", span) subtitle(" ") ///
   
10. note("Source:1988 data from NLS, U.S. Dept. of Labor, "          ///
    
11. "Bureau of Labor Statistics",span)
    
12. exit

复制代码

1. // Labelling the over variable
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. levelsof for, local(a)
   
6. levelsof rep78, local(a1)
   
7. 
   
8. tempfile file
   
9. 
   
10. foreach i of local a {     //for
    
11. foreach i1 of local a1 {   //rep78
    
12. 
    
13. sysuse auto, clear
    
14. 
    
15. count if for==`i' & rep78==`i1'  //making sure there are obs to regress
    
16. 
    
17. if r(N)>=1 {
    
18.    regress mpg weight if for==`i' & rep78==`i1'
    
19.    parmest, norestore
    
20.    egen axis= axis(estimate), label(parm)
    
21.    gen for=`i'
    
22.    gen rep78=`i1'
    
23.    capture append using `file'
    
24.    save `file', replace
    
25.    local kk=`i1'+1
    
26. }
    
27. }
    
28. }  //end loop
    
29. 
    
30. summarize rep78
    
31. foreach i of local a {  //for
    
32.   set obs `=_N+1'
    
33.   replace for=`i' if for==.
    
34.   local max1=`r(max)'
    
35.   replace rep78=`max1'  +1 if rep78==.
    
36. }
    
37. 
    
38. gsort for -rep78
    
39. 
    
40. drop if parm=="_cons"
    
41. gen order=_N-_n
    
42. 
    
43. label define jj 9 "foreign=0       "  ///
    
44.                 3 "foreign=1       "  ///
    
45.                                 0 "rep78_3"        ///
    
46.                                 1 "rep78_4"        ///
    
47.                                 2 "rep78_5"        ///
    
48.                                 4 "rep78_1"        ///
    
49.                                 5 "rep78_2"        ///
    
50.                 6 "rep78_3"        ///
    
51.                                 7 "rep78_4"        ///
    
52.                                 8 "rep78_5"
    
53. 
    
54. label value order  jj
    
55. 
    
56. twoway (scatter order est , msymbol(none))  ///
    
57. (scatter order est )                        ///
    
58. (rcap min95 max95 order, horizontal ),      ///
    
59. ylabel(0/9, valuelabel angle(0)) legend(off) ytitle("rep78 over foreign") ///
    
60. xtitle("Coefficients and confidence interval")
    
61. 
    
62. exit

复制代码

1. // Including a distribution on top of a histogram.
   
2. 
   
3. //The normal option for graph histogram uses the dataset's mean
   
4. //and standard deviation.
   
5. //This method allows you to specify these.  
   
6. 
   
7. sysuse auto, clear
   
8. histogram mpg , addplot(function normalden(x,26,10), ra(10 42))
   
9. 
   
10. *(Nick Cox - Stata list,  Fri, 03 June 2011)

复制代码

1. // Double graph - if there is too much detail for one
   
2. // graph then 2 graphs can be used
   
3. 
   
4. clear
   
5. set more off
   
6. use "http://www.stata-press.com/data/r11/nlswork2.dta"
   
7. describe
   
8. 
   
9. replace year=1900+year
   
10. 
    
11. gen age11= yofd(date(c(current_date)),"DMY") -(1900+birth_yr)
    
12. 
    
13. recode age11 (56/58=1) (59/61=2) (62/64=3) (65/67=4) (68/70=5)
    
14. 
    
15. collapse (count) union  (count) idcode ,by(year age11)
    
16. gen per=(union/idcode)*100
    
17. 
    
18. twoway                                                                   ///
    
19. (line  per year   if age11==1 , lpattern(dash)  connect(l))              ///
    
20. (line  per year   if age11==2 , lpattern(dash_dot)  connect(l))          ///
    
21. (line   per year   if age11==3 , lpattern(dot) connect(l)) ,             ///
    
22. name(a1,replace) legend(  bmargin(l=40)     stack    region(style(none)) ///
    
23. label( 1 "56-58")                                                        ///
    
24. label( 2 "59-61")                                                        ///
    
25. label( 3 "62-64")                                                        ///
    
26. rows(1)   ) xtitle("")                                                   ///
    
27. xlabel(1968 "            1968"  1978 "1978             ", noticks)
    
28. 
    
29. twoway                                                                      ///
    
30. (line  per year   if age11==4 ,  yaxis(1)  lpattern(solid)  connect(l))     ///
    
31. (line  per year   if age11==5 ,  yaxis(1)  lpattern(longdash)  connect(l)), ///
    
32. name(a2,replace) legend(  bmargin(r=40)  stack    region(style(none))       ///
    
33. label(1 "65-67")                                                            ///
    
34. label( 2 "68-70")                                                           ///
    
35. rows(1)   ) xtitle("")                                                      ///
    
36. xlabel(1968 "            1968"  1978 "1978             ", noticks) ///
    
37. yscale(alt) xline(1968)
    
38. 
    
39. graph combine a1 a2 , imargins(zero) ycommon

复制代码

1. sysuse auto, clear
   
2. su mpg
   
3. local width=1
   
4. histogram mpg ,width(`width') freq ///
   
5. addplot(function y =`width'*`=_N'*normalden(x,`r(mean)',`r(sd)'), ra(10 45))
   
6. 
   
7. 
   
8. sysuse auto, clear
   
9. su mpg
   
10. local width=2   //   <-----------根据需要改变
    
11. histogram mpg ,width(`width') freq ///
    
12. addplot(function y =`width'*`=_N'*normalden(x,`r(mean)',`r(sd)'), ra(10 45))

复制代码

注：原do-file修改后