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

1. //heatmap-Graphs a heatmap.
   
2. 
   
3. /*Requires: cbarplot
   
4. to download this program type the following on the Stata
   
5. command line (if not already loaded):*/
   
6. ssc install hmap
   
7. 
   
8. //To run the example: copy the following into a do file and run
   
9. 
   
10. 
    
11.    // producing the data
    
12.        clear all
    
13.        set obs 64
    
14.        generate n=int(uniform()*10)
    
15.        generate x=1+int((_n-1)/8)
    
16.        generate y=1+mod((_n-1),8)
    
17.        label define xlab 1 "one" 2 "two" 7 "seven" 8 "eight"
    
18.        label define ylab 3 "three" 4 "four" 5 "five" 6 "six"
    
19.        label value x xlab
    
20.        label value y ylab
    
21.        table y x [fw=n]
    
22.    
    
23.    hmap x y n

复制代码

1. //mm_plot
   
2. /*mm_plot is a user written Mata graph command that allow you
   
3. to plot a Mata matrix
   
4. Requires: moremata
   
5. to download this program type the following on the Stata command
   
6. line (if not already loaded):.*/  
   
7. 
   
8. ssc install moremata
   
9. 
   
10. //To run the example: copy the following into a do file and run
    
11. 
    
12. 
    
13. set more off
    
14. mata:
    
15. 
    
16. for(i=0;i<=2*pi()*1000; i++) {
    
17. 
    
18.   a=(i/1000)*cos(i/1000)
    
19.   a1=(i/1000)*sin(i/1000)
    
20. 
    
21. if(i==0) {
    
22.   za=a
    
23.   za1=a1
    
24. }
    
25. else {
    
26.   za=za\a
    
27.   za1=za1\a1
    
28. }
    
29.   c=za,za1
    
30. }
    
31. 
    
32. mata: mm_plot(c)
    
33. 
    
34. end
    

复制代码

1. //Fractal
   
2. /*Generates fractals
   
3. 
   
4. Requires: fractal
   
5. to download this program type the following on the Stata
   
6. command line (if not already loaded):*/
   
7. ssc install fractal
   
8. 
   
9. clear
   
10. 
    
11. fractal  ,hs(0,.2,.3,.40,.5,.618,.75,.85,1.0)   ///
    
12. vs(0,.25,.2,.35,.25,.516,.366,.42,.2) hr(0 100) ///
    
13. vr(0 200) iter(3) savegraph
    
14. 
    
15. graph combine _frctl1.gph _frctl2.gph _frctl3.gph

复制代码

1. // riskplot
   
2. /*The riskplot is a graphical aid to the investigation
   
3. of the contributions of risk factors on outcomes of interest.
   
4. (see SJ Vol 10 No 1)
   
5. 
   
6. Requires: riskplot
   
7. to download this program type the following on the Stata command
   
8. line (if not already loaded):*/
   
9. findit gr0044
   
10. /*(The all option ensures that the data (data_riskplot.dta)
    
11. will also be downloaded)
    
12. 
    
13. To run the example: copy the following into a do file and run*/
    
14. 
    
15. 
    
16. use d:\StataMP\ado\personal\data\data_riskplot.dta, clear  //自行调整  
    
17. set scheme sj
    
18. 
    
19. riskplot depr1995 sex sclass if Idep91==1 [pw=wg], path obs        /*
    
20.    */ thick(20) c(. red)  title(Risk plot for subjects with mild   /*
    
21.    */ or severe depression at baseline) subtitle((results using    /*
    
22.    */ sampling weights), margin(b+5)) scale(0.9) ytitle(depression /*
    
23.    */ score 1995) saving(riskplotWG, replace)

复制代码

1. //zmap
   
2. 
   
3. /*zmap graphs (or maps) binned values of a variable z with respect to two
   
4. variables x and y treated as Cartesian coordinates.
   
5. 
   
6. Requires: zmap
   
7. to download this program type the following on the Stata
   
8. command line (if not already loaded):*/
   
9. 
   
10. 
    
11. //To run the example: copy the following into a do file and run
    
12. 
    
13. ssc install zmap
    
14. 
    
15. //修改color为mcolor
    
16. //修改nlswork.dta下载地址
    
17. //修正mxla为xla
    
18. //色彩暂未修正
    
19. use http://www.stata-press.com/data/r11/nlswork.dta, clear
    
20. 
    
21. clear all
    
22. use nlswork.dta,clear
    
23. egen mean = mean(ln_wage), by(age grade)
    
24. egen tag = tag(age grade)
    
25. label variable mean "mean ln wage"
    
26. summarize ln_wage if !missing(age, grade), detail
    
27. 
    
28. 
    
29. zmap mean grade age if tag, breaks(.993 1.166 1.361 1.641 ///
    
30. 1.964 2.275 2.456) ms(S ..) ysc(on) xsc(on) legend(on pos(3) ///
    
31. col(1)) yla(0/18, ang(h)) ytitle(`: var label grade') ///
    
32. xla(15(5)45) note("") mcolor(blue blue*0.5 orange*0.3 orange)   

复制代码

1. //Statplot
   
2. /*Statplot is an alternative to graph hplot, graph hbar and graph dot.
   
3. It allows grouping on the axis and other options.  
   
4. 
   
5. Requires: statplot
   
6. to download this program type the following on the Stata command
   
7. line (if not already loaded):*/
   
8. ssc install statplot
   
9. 
   
10. //To run the example: copy the following into a do file and run
    
11. 
    
12. 
    
13. webuse systolic, clear
    
14. 
    
15. sysuse census, clear
    
16. statplot marriage divorce, over(region) s(sum) xpose varnames

复制代码

1. //tddens
   
2. /*Bivariate kernel density graphs over a grid
   
3. 
   
4. Requires: tddens
   
5. to download this program type the following on the Stata command
   
6. line (if not already loaded):*/
   
7. ssc install tddens
   
8. 
   
9. //To run the example: copy the following into a do file and run
   
10. 
    
11. 
    
12. sysuse auto, clear
    
13. 
    
14. tddens price mpg, s b

复制代码

1. //Deviation plots
   
2. /*Plots each data point relative to some appropriate mean with a marker
   
3. symbol and a spike connecting marker and mean.
   
4. 
   
5. Requires: devnplot
   
6. to download this program type the following on the Stata command
   
7. line (if not already loaded):*/
   
8. 
   
9. ssc install devnplot
   
10. 
    
11. //To run the example: copy the following into a do file and run
    
12. 
    
13. 
    
14. webuse systolic, clear
    
15. 
    
16. version 9: anova systolic drug disease drug*disease
    
17. predict predict
    
18. predict residual, residual
    
19. 
    
20. devnplot systolic drug disease, level(predict) superplines

复制代码

1. //ellip
   
2. /*ellip calculates a confidence ellipse from the elliptically
   
3. distributed variables yvar and xvar, and then graphs the
   
4. confidence ellipse using graph twoway line.
   
5. 
   
6. Requires: ellip
   
7. to download this program type the following on the Stata command
   
8. line (if not already loaded):*/
   
9. ssc install ellip
   
10. 
    
11. //To run the example: copy the following into a do file and run
    
12. 
    
13. sysuse auto, clear
    
14. 
    
15. ellip mpg weight, by(foreign, total legend(off)) ///
    
16. total tlabel(Total as a by-group) plot(scatter mpg weight)

复制代码

1. //Density Plot
   
2. 
   
3. /*Requires: spgrid, spkde, spmap & mylabels
   
4. to download these programs type the following
   
5. on the Stata command line (if not already loaded):*/
   
6. 
   
7. ssc install spgrid
   
8. ssc install spkde
   
9. ssc install spmap
   
10. ssc install mylabels
    
11. 
    
12. *------------------------------------------------
    
13. 
    
14. * 1. Normalize variables in the range [0,1]       
    
15. 
    
16. sysuse "auto.dta", clear
    
17. set more off
    
18. 
    
19. summarize price mpg
    
20. clonevar x = mpg
    
21. clonevar y = price
    
22. replace x = (x-0) / (50-0)
    
23. replace y = (y-0) / (20000-0)
    
24. mylabels 0(10)50, myscale((@-0) / (50-0)) local(XLAB)
    
25. mylabels 0(5000)20000, myscale((@-0) / (20000-0)) local(YLAB)
    
26. 
    
27. keep x y
    
28. save "xy.dta", replace
    
29. 
    
30. * 2. Generate a 100x100 grid
    
31. 
    
32. spgrid, shape(hexagonal) xdim(100) ///
    
33. xrange(0 1) yrange(0 1)  ///
    
34. dots replace   ///
    
35. cells("2D-GridCells.dta")  ///
    
36. points("2D-GridPoints.dta") ///
    
37. 
    
38. * 3. Estimate the bivariate probability density function
    
39. 
    
40. spkde using "2D-GridPoints.dta", ///
    
41. xcoord(x) ycoord(y)                 ///
    
42. bandwidth(fbw) fbw(0.1) dots        ///
    
43. saving("2D-Kde.dta", replace)
    
44. 
    
45. * 4. Draw the density plot
    
46. 
    
47. use "2D-Kde.dta", clear
    
48. recode lambda (.=0)
    
49. spmap lambda using "2D-GridCells.dta",  ///
    
50. id(spgrid_id) clnum(20) fcolor(Rainbow) ///
    
51. ocolor(none ..) legend(off)             ///
    
52. point(data("xy.dta") x(x) y(y))         ///
    
53. freestyle aspectratio(1)                ///
    
54. xtitle(" " "Mileage (mpg)")             ///
    
55. xlab( `XLAB')                           ///
    
56. ytitle(" " "Price ({c S|}US)")          ///
    
57. ylab(`YLAB', angle(0))
    
58. 
    
59. exit

复制代码