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

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. // Stacked bar graph - percentage
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. gen a=_n
   
6. collapse a, by(for rep78)
   
7. drop if missing(for, rep78)
   
8. 
   
9. bysort for rep78: gen sum = sum(a)
   
10. replace a=0 if a==.
    
11. 
    
12. levelsof rep78, local(aa)
    
13. 
    
14. generate csum=.
    
15. 
    
16. foreach i of local aa {
    
17.   summ a if `i'==rep78
    
18. replace csum=a/`r(sum)' if `i'==rep78
    
19. }
    
20. 
    
21. graph bar (sum)csum ,  over(for) over(rep78)        ///
    
22. asyvars stack legend(size(vsmall)) per  blabel(bar, ///
    
23. position(center) format(%3.1f))
    
24. 
    
25. 
    
26. *(Also see the user written program catplot on:
    
27. http://www.survey-design.com.au/Usergraphs.html
    
28. )  

复制代码

1. // Editing a graph with the gr_edit command
   
2. 
   
3. //The graph with the y axis line
   
4. 
   
5. sysuse auto, clear
   
6. graph hbar   weight price, yscale(off)
   
7. 
   
8. //The graph with the y axis line removed
   
9. 
   
10. sysuse auto, clear
    
11. graph hbar   weight price, ysc(off)
    
12. gr_edit .varaxis.style.editstyle linestyle(color(none)) editcopy
    
13. 
    
14. //The above can also be produced with:
    
15. 
    
16. graph hbar weight price ,showyvar  yvaroptions(axis(lc(none)))
    
17. 
    
18. From the Statalist:
    
19. Kieran McCaul
    
20. 11/6/2010
    
21. 
    
22. David Airey
    
23. 11/6/2010

复制代码

1. // Breaking up text automatically in a long note
   
2. 
   
3. webuse auto, clear
   
4. set more off
   
5. 
   
6. //input the text to place in the graph - can be on one line if you wish
   
7. 
   
8. local z1 "The development of Linux is one of the most prominent"
   
9. local z2 " examples of free  and  open source software"
   
10. local z3 "local z3 collaboration; typically all the underlying source code"
    
11. local z4 "can beused, freely  modified, and redistributed, both commercially"
    
12. local z5 "and  non-commercially, by anyone under licenses such as the GNU"
    
13. local z6 "General Public License."
    
14. local z7 "Typically Linux is packaged in a format known as a Linux"
    
15. local z8 "distribution for desktop and server use."
    
16. 
    
17. local total "`z1' `z2' `z3' `z4' `z5' `z6' `z7' `z8'"
    
18. 
    
19. local size=40   //adjust for length of line
    
20. 
    
21. local t1
    
22. 
    
23. forvalues i=1/12 {                //increase the 12 if more lines required
    
24. local a`i' : piece `i' `size' of "`total'"   //extended macro function
    
25. 
    
26. local t1  `"`t1' `=char(34)+ "`a`i''"+char(34)' "'
    
27. }
    
28. 
    
29. graph bar price,bargap(-30) over(foreign)         ///
    
30. note( `t1' ,   margin(large) justification(left))
    
31. 
    
32. exit

复制代码

1. // Breaking up long labels
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. #delimit
   
6. label define origin 0 `"group0 group0 group0 group0 group0
   
7. group0 group0 group0 group0 group0"', modify;
   
8. 
   
9. label define origin 1 "group1 group1 group1 group1 group1
   
10. group1 group1 group1 group1 group1 ", modify;
    
11. #delimit cr
    
12. 
    
13. 
    
14. local relabels
    
15. local relabels1
    
16. 
    
17. levelsof for, local(groups)
    
18.   
    
19. local s_len=20  //Change to required length
    
20. 
    
21. foreach g of local groups {
    
22. 
    
23.   local label : label origin `g'
    
24.   local len : length local label
    
25. 
    
26.   if `len'>`s_len' {
    
27.         forvalues i=1/  `=`len'/`s_len'+1' {
    
28.             local p1 : piece `i' `s_len' of `"`label'"', nobreak
    
29.             local relabels `"`relabels' `=char(34)+ "`p1'" +char(34)' "'
    
30.         }
    
31. 
    
32.         local relabels1  `relabels1'   `=`g'+1'    `"`relabels'"'
    
33. local relabels
    
34.        
    
35.    }
    
36. }
    
37. 
    
38. graph hbar mpg, over(for, relabel(`relabels1' ))
    
39. exit

复制代码

1. *<-------------修改后
   
2. 
   
3. ssc install splitvallabels
   
4. sysuse auto, clear
   
5.    
   
6. label define rep78 1 "Very low rating in 1978" ///
   
7. 2 "Mediocre - Low Rating in 1978"              ///
   
8. 3 "Medium Rating in 1978"                      ///
   
9. 4 "Not bad rating in 1978, considering"        ///
   
10. 5 "Excellent in 1978--the best of the year!"
    
11.         
    
12. label values rep78 rep78
    
13.         
    
14. splitvallabels rep78 if rep78>=3, length(11) recode
    
15.         
    
16. graph bar mpg if rep78>=3, over(rep78, relabel(`r(relabel)') ) scheme(s1)
    
17. 
    
18. exit

复制代码

1. //Combining many graphs on a page
   
2. sysuse auto, clear
   
3. graph drop _all
   
4. drop make
   
5. foreach i of varlist _all{
   
6.         capture confirm numeric variable `i' if _rc==0
   
7.         histogram `i', name(`i')
   
8.         local z "`z' `i'"
   
9.         graph combine `z'
   
10. }

复制代码

*修改后的do-file

1. //Changing the order of twoway graphs where the by() option is used
   
2. 
   
3. *Requires: egenmore
   
4. *to download this program type the following on the Stata command
   
5. *line (if not already loaded):
   
6. ssc install egenmore
   
7. 
   
8. //the following is an example of order that the graphs are in when
   
9. //determined by Stata
   
10. // ie. lowest rep78 value to highest
    
11. sysuse auto,clear
    
12. scatter mpg weight , by(rep78) name(a2, replace)
    
13. 
    
14. //If you require the order to be the mean mpg of the levels of rep78,
    
15. //the following can be done:  
    
16.   
    
17. egen mean1 = mean(mpg) , by(rep78)
    
18. 
    
19. egen axis = axis1(mean1 rep78), label(rep78) reverse  // egenmore command
    
20. 
    
21. scatter mpg weight , by(axis) name(a1, replace)

复制代码

1. //Using the compact option in the by() option
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. forvalues i=10(10)40 {
   
6.   local m`i'=string(`i'*0.4251437,"%6.2f")
   
7. }
   
8. 
   
9. twoway (scatter mpg weight, yaxis(1 2) xlabel(,angle(45)) ///
   
10. yscale( axis(2, alt))                                     ///
    
11. ylabel(                                                   ///
    
12. 10 "`m10'"                                                ///
    
13. 20 "`m20'"                                                ///
    
14. 30 "`m30'"                                                ///
    
15. 40 "`m40'"                                                ///
    
16. ,axis(2))) ,                                              ///
    
17. by(for, title( "km/L",  orientation(vertical)            ///
    
18. size(med) position(3)) xrescale compact)  

复制代码

1. //Adding titles to the left hand side of the graph
   
2. 
   
3. 
   
4. sysuse auto, clear
   
5. 
   
6. fillin rep78 foreign
   
7. twoway scatter price mpg, by(foreign rep78, cols(5) compact)
   
8. 
   
9. egen group = group(foreign rep78)
   
10. 
    
11. label define group 1 "Poor" ///
    
12. 2 "Fair"         ///
    
13. 3 "Average"      ///
    
14. 4 "Good"         ///
    
15. 5 "Excellent"    ///
    
16. 6" "             ///
    
17. 7" "             ///
    
18. 8" "             ///
    
19. 9" "             ///
    
20. 10 " "  
    
21. 
    
22. 
    
23. label value group group ///
    
24. 
    
25. 
    
26. twoway scatter price mpg, ///
    
27. by(group, cols(5)         ///
    
28. r1title("Car type",       ///
    
29. orientation(rvertical)    ///
    
30. size(medsmall))           ///
    
31. t1title("Repair status",  ///
    
32. size(medsmall))           ///
    
33. note("") compact)  
    
34. 
    
35. 
    
36. gr_edit plotregion1.r1title[5].style.editstyle drawbox(yes) editcopy
    
37. 
    
38. gr_edit plotregion1.r1title[5].as_textbox.setstyle, style(yes)
    
39. gr_edit plotregion1.r1title[10].as_textbox.setstyle, style(yes)
    
40. 
    
41. gr_edit plotregion1.r1title[5].text = {}
    
42. gr_edit plotregion1.r1title[5].text.Arrpush Domestic
    
43. 
    
44. gr_edit plotregion1.r1title[10].text = {}
    
45. gr_edit plotregion1.r1title[10].text.Arrpush Foreign
    
46. 
    
47. gr_edit plotregion1.r1title[5]._set_orientation rvertical
    
48.                                        

复制代码