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

1. //Predictive margins with 95% confidence interval - using margins
   
2. 
   
3. /*Requires: parmest
   
4. to download this program type the following on the Stata command
   
5. line (if not already loaded):*/
   
6. ssc install parmest
   
7. 
   
8. sysuse auto, clear
   
9. gen kprice=price/1000
   
10. regress  mpg c.kprice##c.kprice i.foreign i.rep78 turn trunk headroom
    
11. 
    
12. margins, at(kprice=(3/15)) vsquish post
    
13. 
    
14. return list
    
15. ereturn list
    
16. 
    
17. parmest, norestore  //saves parameters
    
18. l
    
19. gen kprice=3 + _n - 1
    
20. l
    
21. label var kprice "Price(\$1000)"
    
22. 
    
23. twoway ///
    
24. (rarea max95 min95 kprice, pstyle(ci)) ///
    
25. (line estimate kprice), ///
    
26. ytitle(Miles per gallon) ///
    
27. title(Predictive margins with 95% confidence interval) ///
    
28. xlabel(3(3)15 10.4) legend(off)

复制代码

1. //Paired Plot 2
   
2. //Looking at the effect of treatment on the weight of anorexic patients.
   
3. 
   
4. //Making up the data; based on the data found at:
   
5. //http://users.polisci.wisc.edu/kritzer/Teaching/ps551/anorexia.txt
   
6. 
   
7. clear
   
8. set obs 76
   
9. gen treatment=mod(_n,3)+1
   
10. 
    
11. gen w_before1=rnormal(82, 4.8)
    
12. replace w_before1=rnormal(81, 5.7)  if treatment==2
    
13. replace w_before1=rnormal(83, 5)    if treatment==3
    
14. 
    
15. gen w_after1=rnormal(85.6, 8.3)
    
16. replace w_after1=rnormal(81.1, 4.7) if treatment==2
    
17. replace w_after1=rnormal(90, 5.4)   if treatment==3
    
18. //end of making up the data
    
19. 
    
20. label var w_before " Weight before treatment, lb"
    
21. label var w_after  " Weight after treatment, lb"
    
22. label define treat ///
    
23.         1  "Cognitive behavioural" ///
    
24.         2  "Control"               ///
    
25.         3  "Family therapy"
    
26.                
    
27. label values treatment treat
    
28.        
    
29. gen byte one = 1
    
30. gen byte two = 2       
    
31.        
    
32. twoway ///
    
33. pcspike w_before one w_after two if w_before <= w_after, ///
    
34. lcolor(gs12) ||                                          ///
    
35. pcspike w_before one w_after two if w_before > w_after,  ///
    
36. lw(*1.2) lcolor(gs2) xla(1 " before" 2 "after ")         ///
    
37. xtitle("") yla(, nogrid ang(h)) ytitle("weight, lb")     ///
    
38. by(treatment, row(1) note("") legend(off) noixtick)       

复制代码

1. //Paired Plot 3
   
2. //Looking at the effect of treatment on the weight of anorexic patients.
   
3. 
   
4. //Making up the data; based on the data found at:
   
5. //http://users.polisci.wisc.edu/kritzer/Teaching/ps551/anorexia.txt
   
6. 
   
7. clear
   
8. set obs 76
   
9. gen treatment=mod(_n,3)+1
   
10. 
    
11. gen w_before1=rnormal(82, 4.8)
    
12. replace w_before1=rnormal(81, 5.7)  if treatment==2
    
13. replace w_before1=rnormal(83, 5)    if treatment==3
    
14. 
    
15. gen w_after1=rnormal(85.6, 8.3)
    
16. replace w_after1=rnormal(81.1, 4.7) if treatment==2
    
17. replace w_after1=rnormal(90, 5.4)   if treatment==3
    
18. //end of making up the data
    
19. 
    
20. label var w_before " Weight before treatment, lb"
    
21. label var w_after  " Weight after treatment, lb"
    
22. label define treat ///
    
23.         1  "Cognitive behavioural" ///
    
24.         2  "Control"               ///
    
25.         3  "Family therapy"
    
26.                
    
27. label values treatment treat
    
28.        
    
29. 
    
30. bysort treatment (w_before w_after) : gen order1 = _n - _N/2
    
31.        
    
32. twoway                                                   ///
    
33. pcarrow w_before order1 w_after order1, pstyle(p1)       ///
    
34. || scatter w_before order1, pstyle(p1) ms(o)             ///
    
35. xla(none) xtitle("") yla(, ang(h)) ytitle("weight, lb")  ///
    
36. by(treatment, row(1) note("") legend(off))
    
37.        

复制代码

1. // Adding labels - standard error bar chart
   
2. 
   
3. webuse assembly, clear
   
4. 
   
5. generate upper=mean+2*std
   
6. 
   
7. serrbar mean std date, scale(2) yline(195)              ///
   
8. addplot( scatter  up date, msymbol(none) mlabel(upper))

复制代码

1. // Odds ratio graph - selecting order of bars
   
2. 
   
3. ** requires 2 user written graphs: parmest and egenmore
   
4. ** if not already installed:
   
5. 
   
6. ssc install parmest
   
7. ssc install egenmore
   
8. 
   
9. sysuse nlsw88, clear
   
10. logit union married never_married collgrad south
    
11. 
    
12. parmest, norestore eform
    
13. 
    
14. tab estimate, gen(kk)
    
15. 
    
16. //order of bars
    
17. input a
    
18. 1
    
19. 2
    
20. 3
    
21. 4
    
22. 5
    
23. 
    
24. egen axis= axis( a estimate), label(parm)
    
25. 
    
26. twoway bar  estimate axis, base(1) horizontal barw(.5) ///
    
27. xline(1) ylab(, valuelabel angle(0)) ||                ///
    
28. rcap min95 max95 axis, horizontal                      ///
    
29. legend(order(1 "point estimate"                        ///
    
30. 2 "95% conf. int.") pos(6))                            ///
    
31. xtitle("odds ratio") ytitle("") scheme(color)
    
32. 
    
33. *(Maarten L. Buis - Statalist,  14/09/2010)

复制代码

1. // Odds ratio graph
   
2. 
   
3. ** requires 2 user written graphs: parmest and egenmore
   
4. ** if not already installed:
   
5. 
   
6. ssc install parmest
   
7. ssc install egenmore
   
8. 
   
9. sysuse nlsw88, clear
   
10. logit union married never_married collgrad south
    
11. 
    
12. parmest, norestore eform
    
13. egen axis= axis(estimate), label(parm)
    
14. 
    
15. twoway scatter   axis estimate, ///
    
16.             xline(1) ylab(, valuelabel angle(0)) ||    ///
    
17.        rcap min95 max95 axis, horizontal               ///
    
18.             legend(order(1 "point estimate"            ///
    
19.                          2 "95% conf. int.") pos(6))   ///
    
20.             xtitle("odds ratio") ytitle("")

复制代码

1. // Table on the side of a graph
   
2. 
   
3. sysuse auto, clear
   
4. twoway (scatter mpg rep78) , graphregion( margin(r+30) )
   
5. 
   
6. gr_edit AddTextBox added_text editor `=82+8' `=101'
   
7. gr_edit added_text_new = 1
   
8. gr_edit added_text_rec = 1
   
9. gr_edit added_text[1].text = {}
   
10. gr_edit added_text[1].text.Arrpush "Mean mpg by rep78"
    
11. 
    
12. gr_edit AddTextBox added_text editor `=82+2' `=112'
    
13. gr_edit added_text_new = 2
    
14. gr_edit added_text_rec = 2
    
15. gr_edit added_text[2].text = {}
    
16. gr_edit added_text[2].text.Arrpush "rep78  mpg"
    
17. 
    
18. collapse mpg, by(rep78)
    
19. 
    
20. local z=2
    
21. 
    
22. forvalues i=0/4 {  //rows
    
23.   local ++z
    
24.   gr_edit AddTextBox added_text editor `=80-(`i'*4)' `=114'
    
25.   gr_edit added_text_new = `z'
    
26.   gr_edit added_text_rec = `z'
    
27.   gr_edit added_text[`z'].text = {}
    
28.   gr_edit added_text[`z'].text.Arrpush ///
    
29.   "`=rep78[`z'-2]'      `=string(mpg[`z'-2],"%8.2f")' "
    
30. }
    
31. exit

复制代码

1. // Table on the side of a graph
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. forval n = 1/20 {
   
6.         loc sidetable `" `sidetable'  " `=make[`n']'  "  "'
   
7.         }
   
8. 
   
9. forval n = 1/20 {
   
10.         loc sidecap `" `sidecap'  "      =price[`n']'"  "'
    
11.         }
    
12.        
    
13. sysuse auto, clear
    
14. twoway (scatter mpg rep78) ,        ///
    
15. note("{bf: MAKE    }"               ///
    
16. "-----------------"    `sidetable'  ///
    
17. , size(medsmall) color(green)       ///
    
18. position(2) margin(small)           ///
    
19. justification(left))                ///
    
20. caption(" {bf:   Price }"           ///
    
21. "---------"    `sidecap'            ///
    
22. , size(medsmall) color(midgreen)    ///
    
23. position(2) margin(vsmall)          ///
    
24. justification(right))
    
25. exit
    
26. 
    
27. //Based statalist Jun 9, 2011
    
28. //Eric Booth

复制代码

1. // Putting Data set results in a graph
   
2. 
   
3. clear all
   
4. sysuse auto
   
5. 
   
6. quiet: summarize mpg
   
7. 
   
8. twoway ///
   
9. scatter mpg rep78, ///
   
10. text(40 3 "The mean mpg is:`=string(`=r(mean)',"%8.2f")'", placement(right)) ///
    
11. text(38 3 "The sd mpg is: `=string(`=r(sd)',"%8.2f")'" ,   placement(right))
    
12. 
    
13. exit

复制代码

1. sysuse auto, clear
   
2. 
   
3. egen mpg1=cut(mpg), group(5) icode
   
4. 
   
5. levelsof mpg1, local(kk)
   
6. 
   
7. gen mean=.
   
8. gen ub=.
   
9. gen lb=.
   
10. gen loc=.
    
11. local z=1
    
12. 
    
13. foreach i of local kk {
    
14. ci for if mpg1==`i', wilson binomial
    
15. 
    
16.   replace mean=r(mean) in `z'
    
17.   replace ub=r(ub) in `z'
    
18.   replace lb=r(lb) in `z'
    
19.   replace loc=`i' in `z'
    
20.   local ++z
    
21. }
    
22. 
    
23. label define a 0 "very poor"  1 poor  2  OK 3 good  4 vgood
    
24. label value  loc a
    
25. 
    
26. twoway (bar mean loc) (rcap ub lb loc), ///
    
27. xlab(, valuelabel angle(45))            ///
    
28. xtitle("Fuel Consumption")              ///
    
29. ytitle(mean) legend(off)

复制代码