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

1. //strippplot
   
2. 
   
3. ssc install stripplot
   
4. 
   
5. sysuse bplong, clear
   
6. 
   
7. egen group = group(age sex), label
   
8. 
   
9. stripplot bp*, bar over(when) by(group, compact col(1) note("")) ///
   
10. ysc(reverse) subtitle(, pos(9) ring(1) nobexpand bcolor(none) ///
    
11. placement(e)) ytitle("") xtitle(Blood pressure (mm Hg))
    
12. 
    
13. //From striplot help:Nick cox

复制代码

1. //strippplot
   
2. 
   
3. ssc install stripplot
   
4. 
   
5. 
   
6. sysuse auto, clear
   
7. gen pipe = "|"
   
8. stripplot price, over(rep78) box(barw(0.3)) ms(none) mla(pipe) boffset(0.3)
   
9. 
   
10. //From the Statalist:Nick cox 15/10/2010

复制代码

lexis.rar (796 Bytes) 本附件包括：

• lexis.do

2015-9-9 09:55:01 上传

1. //sixplot
   
2. 
   
3. /*Displays six diagnostic and descriptive graphs for a single variable
   
4. 
   
5. Requires: sixplot
   
6. to download this program type the following on the Stata command line
   
7. (if not already loaded):*/
   
8. ssc install sixplot
   
9. 
   
10. sysuse uslifeexp.dta
    
11. 
    
12. sixplot le_male

复制代码

netplot-do.rar (505 Bytes) 本附件包括：

• netplot-do.do

2015-9-9 10:15:09 上传

1. //添加surf3数据集来源，findit surf3下载
   
2. //修正surface为surface7
   
3. 
   
4. use surf3,clear
   
5. set more off
   
6. surface7 x y z, orient(xzy) saving(s3d ,replace)

复制代码

z注：Z-axis刻度反向，待修改。

1. /*binscatter_非参估计_其中添加mcol(green)
   
2. binscatter generates binned scatterplots, and is optimized for
   
3. speed in large datasets.
   
4. Binned scatterplots provide a non-parametric way of visualizing the
   
5. relationship between two variables. With a large number of observations,
   
6. a scatterplot that plots every data point would become too crowded to
   
7. interpret visually.  binscatter groups the x-axis variable into
   
8. equal-sized bins, computes the mean of the x-axis and y-axis variables
   
9. within each bin, then creates a scatterplot of these data points.*/
   
10. ssc install binscatter
    
11. clear all
    
12. 
    
13. sysuse nlsw88, clear
    
14. 
    
15. keep if inrange(age,35,44) & inrange(race,1,2)
    
16. 
    
17. scatter wage tenure ,  mcol(green) title("Graph produced by Scatter") name(g1)
    
18.         
    
19. /* The scatter was too crowded to be easily interpetable. The binscatter
    
20.    is cleaner, try a quadratic fit.*/
    
21. 
    
22. binscatter wage tenure, line(qfit) xscale(range(0,25)) ///
    
23. xlabel(0(5)25) ylabel(0(10)40) yscale(range(0,40))  ///
    
24. title("Graph produced by Binscatter") name(g2)
    
25. 
    
26. graph combine g1 g2

复制代码

1. ssc install contour
   
2. 
   
3. clear
   
4. set obs 40
   
5. gen a = 0.20 if [_n] == 1
   
6. replace a = a[_n-1] - 0.02 if [_n] != 1
   
7. gen b = a
   
8. fillin a b
   
9. gen pred = a*-1067.54 + b*-89.1
   
10. gen _90 = pred >= 90
    
11. 
    
12. twoway (contour _90 a b, levels(2) scolor(mint) ecolor(magenta) heatmap)
    
13. 
    
14. twoway (contour _90 a b, scolor(mint) ecolor(magenta) levels(2) heatmap ///
    
15.     zlabel(.25 "0" .75 "1"))

复制代码

好贴，谢谢分享。

xddlovejiao1314 发表于 2015-9-23 13:53
好贴，谢谢分享。

感谢师兄哈