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

1. /*Producing a similar plot to the density plot at:
   
2. http://www.survey-design.com.au/Usergraphs.html
   
3. using Stata 12's new contour plot
   
4. 
   
5. Requires: spgrid, spkde  & mylabels
   
6. to download these programs type the following
   
7. on the Stata command line (if not already loaded):
   
8. 耗时约3分钟
   
9. */
   
10. 
    
11. ssc install spgrid
    
12. ssc install spkde
    
13. ssc install mylabels
    
14. 
    
15. sysuse "auto.dta", clear
    
16. set more off
    
17. 
    
18. summarize price mpg
    
19. clonevar x = mpg
    
20. clonevar y = price
    
21. replace x = (x-0) / (50-0)
    
22. replace y = (y-0) / (20000-0)
    
23. mylabels 0(10)50, myscale((@-0) / (50-0)) local(XLAB)
    
24. mylabels 0(5000)20000, myscale((@-0) / (20000-0)) local(YLAB)
    
25. keep x y
    
26. save "xy.dta", replace
    
27. 
    
28. * 2. Generate a 100x100 grid
    
29. 
    
30. spgrid, shape(hexagonal) xdim(100)    ///
    
31. xrange(0 1) yrange(0 1)               ///
    
32. dots replace                          ///
    
33. cells("2D-GridCells.dta")             ///
    
34. points("2D-GridPoints.dta")
    
35. 
    
36. * 3. Estimate the bivariate probability density function
    
37. 
    
38. spkde using "2D-GridPoints.dta",    ///
    
39. xcoord(x) ycoord(y)                 ///
    
40. bandwidth(fbw) fbw(0.1) dots        ///
    
41. saving("2D-Kde.dta", replace)
    
42. 
    
43. use "2D-Kde.dta", clear
    
44. 
    
45. merge 1:1 _n using xy.dta
    
46. 
    
47. twoway (contour p  spgrid_ycoord spgrid_xcoord  if p!=0 , ///
    
48. levels(15))                                               ///
    
49. (scatter y x, mcolor(black) msize(small) )                ///
    
50. ,xlab(`XLAB', nogrid) xtitle("Mileage (mpg)")             ///
    
51. ylab(`YLAB', nogrid)                                      ///
    
52. ytitle("Price $US") plotregion(color(blue))

复制代码

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

复制代码

1. // Different color bars
   
2. 
   
3. clear
   
4. 
   
5. input ///
   
6. id w pos mark
   
7. 1  1 1  69.55
   
8. 2  1 2  65.16
   
9. 3  1 3  64.91
   
10. 4  1 4  64.53
    
11. 5  1 5  63.70
    
12. 6  0 6  84.58
    
13. 7  0 7  84.51
    
14. 8  0 8  84.12
    
15. 9  0 9  83.34
    
16. 10 0 10 82.8
    
17. end
    
18. 
    
19. label define kk        ///
    
20. 1  "Barbora - CZE "    ///
    
21. 2  "Christina - GER "  ///
    
22. 3  "Linda - GER "      ///
    
23. 4  "Sunette - RSA"     ///
    
24. 5  "Huihui - CHN"      ///
    
25. 6  "Keshorn - TRI "    ///
    
26. 7  "Oleksandr - UKR "  ///
    
27. 8  "Antti - FIN "      ///
    
28. 9  "Vitezslav - CZE"   ///
    
29. 10 "Tero - FIN"
    
30. 
    
31. label value pos kk
    
32. 
    
33. label define w 1 "Women" 0 "Men"
    
34. label value w w
    
35. 
    
36. graph bar (asis) mark,  over(id ) over(pos,label(angle(45))) over(w) ///
    
37. ylabel( 0(10)90, angle(45))                                         ///
    
38. blabel(bar, position(inside) format(%9.1f) color(black))            ///
    
39. scheme(s2color)  legend(off)                                        ///
    
40. bargap(5) title("London Olympics 2012" "Javelin")                   ///
    
41. nofill ytitle("Metres")                                             ///
    
42. bar(1, color(gold))      ///
    
43. bar(6, color(gold))      ///
    
44. bar(2, color(gs13))      ///
    
45. bar(7, color(gs13))      ///
    
46. bar(3, color(sienna))    ///
    
47. bar(8, color(sienna))    ///
    
48. bar(4, color(cyan))      ///
    
49. bar(9, color(cyan))      ///
    
50. bar(5, color(pink))      ///
    
51. bar(10, color(pink))
    
52. 
    
53. 
    
54. exit

复制代码

1. // Stacking some levels of a variable and not others
   
2. // Has the years 71 and 77 stacked and year 88 normaly displayed
   
3. 
   
4. webuse union, clear
   
5. keep idcode year age
   
6. recode year (71 77=1) (88=2) (*=.), gen(which)
   
7. graph bar age, over(year) over(which) asyvars stack
   
8. 
   
9. exit
   
10. 
    
11. Statalist
    
12. Nick Cox 22/12/11

复制代码

1. // bar graph - categorical variables - proportions
   
2. 
   
3. sysuse auto, clear
   
4. //creating some data
   
5. gen mpg1=mpg<30
   
6. recode mpg (1/25=1) (26/30=2) (31/50=3), gen(mpg2)
   
7. save kk, replace
   
8. list
   
9. 
   
10. capture erase kk1.dta
    
11. local z=1
    
12. 
    
13. foreach i of varlist for mpg1 mpg2 {
    
14.   use kk, clear
    
15.   keep `i'
    
16.   contract `i'
    
17.   egen total=total(_freq)
    
18.   replace _freq=_freq/total
    
19.   generate type=`z'
    
20.   capture append using kk1
    
21.   save kk1, replace
    
22.   local ++z
    
23. }
    
24. 
    
25. drop mp* for
    
26. bysort type : gen t1=_n
    
27. 
    
28. label define kk 1 cat1 2 cat2 3 cat3
    
29. label value type kk
    
30. 
    
31. graph bar (mean) _freq , over(t1,sort(_freq) ) over(type, ) showyvars ///
    
32. asyvar legend(off) nofill ytitle("Proportion") scheme(s2color)

复制代码

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

复制代码

1. // Controlling date labels
   
2. clear
   
3. 
   
4. input str40 time       x
   
5. "5jan2011 15:00:30"    1
   
6. "1feb2011 17:30:32"    2
   
7. "1march2011 18:10:30"  6
   
8. "1april2011 05:05:01"  2
   
9. "1march2011 15:10:30"  12
   
10. end
    
11. list
    
12. 
    
13. generate double time1=clock(time,"DMY hms")
    
14. 
    
15. local t1=tc(1jan2011 12:00:0)
    
16. local t2=1000*60*60*24*10
    
17. local t3=tc(1may2011 12:0:0)
    
18. 
    
19. scatter x time1, xlab(`t1'(`t2')`t3' ,                ///
    
20. format(%tcddmonccyy_hh:mm) angle(45)) xtitle("Time")
    
21. 
    
22. *(Nick Cox - Stata list,  Monday, 23 Jan 2011)

复制代码

1. //2nd Row of axis labels
   
2. 
   
3. webuse turksales, clear
   
4. line sales t,                              ///
   
5.     xtick(119.5(4)159.5, tlength(*1.5))    ///
   
6.     xmlabel(120/159,                       ///
   
7.         noticks format(%tqq))              ///
   
8.     xlabel(121.5(4)157.5,                  ///
   
9.         noticks format(%tqCY) labgap(2.5)) ///
   
10.     xtitle("")
    
11. 
    
12. *(Maarten Buis - Stata list,  Thu, 7 Nov 2013)

复制代码

1. // Adding frequency count bar labels in bold font
   
2. 
   
3. // Adding this using the Stata editor commands. This is produced
   
4. // using a program.
   
5. 
   
6. program g2
   
7. syntax varname, over(varname)
   
8. 
   
9. levelsof `over', local(kk)
   
10. 
    
11. graph bar  `varlist'  , over(`over')                 ///
    
12. blabel(bar, position(base) gap(*22.5) format(%9.1f)  ///
    
13. size(large) color(black))
    
14. 
    
15. local x=1
    
16. foreach i of local kk {
    
17.   summ `varlist' if `over'==`i'
    
18.   local a1="N= "+string(r(N),"%8.0f")
    
19. 
    
20. gr_edit plotregion1.barlabels[`x'].text = {}
    
21. gr_edit plotregion1.barlabels[`x']._set_orientation vertical
    
22. gr_edit plotregion1.barlabels[`x'].text.Arrpush `"{bf:  `a1'   }"'
    
23. local ++x
    
24. }
    
25. 
    
26. end
    
27. 
    
28. sysuse auto, clear
    
29. egen pr2=cut(price), group(10)
    
30. 
    
31. g2 weight , over(pr2)

复制代码

1. // Coloring bars in order of magnitude. All bars not covered by the first and
   
2. // second places are considered (colored) as 3rd place)
   
3. 
   
4. clear
   
5. input country gdp
   
6. 1        100       
   
7. 2        110       
   
8. 3        240       
   
9. 4        50       
   
10. 5        10       
    
11. end
    
12. 
    
13. list
    
14. 
    
15. gsort -gdp
    
16. 
    
17. gen gdpred=gdp if _n==1
    
18. gen gdpgreen=gdp if _n==2
    
19. gen gdpblue=gdp if _n>2
    
20. 
    
21. list
    
22. 
    
23. graph hbar gdpred gdpgreen gdpblue, over(country, sort(gdp)          ///
    
24.       descending) nofill  legend(label(1 "Best") label(2 "2nd best") ///
    
25.           label(3 "rest") rows(1)) bar(1, color(cranberry))              ///
    
26.           bar(2, color(dkgreen)) bar(3, color(dknavy))
    
27. 
    
28. *(Sergiy Radyakin - Stata list,  Sat, 19 Feb 2011)

复制代码