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

clear input str14 country school1 school2 agemarriage1 agemarriage2 bangladesh 3.3 4.5 14.8 15.2 bolivia 6.9 7.6 19.8 20.1 colombia 7.9 8.6 20.3 20.3 dr 7.9 8.6 18.3 18.3 egypt 5.5 7.3 18.9 19.7 haiti 3.1 4.3 19.6 19.4 india 3.6 4.3 16.9 17.1 indonesia 5.9 7.5 18.1 19.3 morocco 1.6 2.7 18.7 19.6 nepal 2.4 3.6 16.9 17.4 pakistan 1.5 2.8 17.9 18.4 peru 8.1 8.8 20.2 20.6 end set more off keep if inlist(country,"bangladesh","bolivia","colombia","dr","egypt","haiti","india","indonesia") collapse (mean) school2, by(school1 country) sort country school1 generate date1=_n tsset date1  //set to time series so lag "l." and              // forward "f." operators can be used                          // see help tsvarlist                          twoway pcarrow school1 l.school1 f.school1 school1, mlab(date1) /// ytitle("This month Deg. C")                            /// xtitle("Previous month Deg. C")                        /// title("Temperature Cycle Melbourne - Deg. C")   /// subtitle("Years 1940 1950 1960 1970 2010") 复制代码

一个粗糙的例子

1. insheet using "IDCJAC0002_087031_Data1.csv", clear  
   
2. //or whatever data you have
   
3. set more off
   
4. keep if inlist(year,1940, 1950, 1960, 1970,2010)
   
5. collapse (mean) MeanmaximumtemperatureC, by(month year)
   
6. sort year month
   
7. generate date1=_n
   
8. 
   
9. tsset date1  //set to time series so lag "l." and
   
10.              // forward "f." operators can be used
    
11.                          // see help tsvarlist
    
12.                         
    
13. twoway pcarrow MeanmaximumtemperatureC l.MeanmaximumtemperatureC f.MeanmaximumtemperatureC MeanmaximumtemperatureC, mlab(date1) ///
    
14. ytitle("This month Deg. C")                            ///
    
15. xtitle("Previous month Deg. C")                        ///
    
16. title("Temperature Cycle Melbourne - Deg. C")   ///
    
17. subtitle("Years 1940 1950 1960 1970 2010")

复制代码

IDCJAC0002_087031.zip (8.16 KB) 本附件包括：

• IDCJAC0002_087031_Data1.csv
• IDCJAC0002_087031_Data12.csv
• IDCJAC0002_087031_Note.txt

2015-9-1 15:27:27 上传

数据csv

                                                 对元数据及do-file变量进行了修改

1. webuse nlswork, clear
   
2. 
   
3. local ytitle "`: var label ln_wage'"
   
4. 
   
5. set seed 2803
   
6. generate rnd = runiform()
   
7. bysort idcode (rnd): generate byte select = _n == 1
   
8. statsby mean=r(mean) ub=r(ub) lb=r(lb) N=r(N) if select, ///
   
9. by(grade) clear: ci ln_wage
   
10. 
    
11. levelsof grade, local(levels)
    
12. 
    
13. twoway rcap ub lb grade || ///
    
14. scatter mean grade, ytitle(`ytitle') legend(off) ///
    
15. subtitle("95% confidence intervals for mean", place(w)) xlabel(`levels')

复制代码

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. // Stacked bar graph - reshaping the data
   
2. 
   
3. //creating the data
   
4. 
   
5. clear
   
6. set obs 1000
   
7. egen id = seq(), block(10)
   
8. egen time = seq(), to(10)
   
9. format time %tq
   
10. gen time2=string(time, "%tq")
    
11. gen sector = ceil(3 * runiform())
    
12. bysort id: replace sector=sector[1]
    
13. label define sector 1 "SECTOR 1" 2 "SECTOR 2" 3 "SECTOR 3"
    
14. label values sector sector
    
15. gen sales=exp(rnormal())*1000
    
16. gen sales2=2*sales
    
17. 
    
18. //reshaping the data
    
19. gen obs = _n
    
20. rename sales sales1
    
21. reshape long sales, i(obs) j(which)
    
22. 
    
23. graph bar (sum) sales ,  over(sector) over(which) ///
    
24. over(time2, label(labsize(vsmall))) asyvars stack legend(size(vsmall))

复制代码

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 bar  estimate axis, base(1) horizontal barw(.5) ///
    
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("")
    
21. 
    
22. *(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. //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. //sequence plot  序列图
   
2. *修改数据下载地址
   
3. ssc install sq
   
4. use http://fmwww.bc.edu/repec/bocode/y/youthemp.dta,clear
   
5. reshape long st, i(id) j(order)
   
6. sqset st id order
   
7. sqindexplot, scheme(vg_palec)

复制代码

1. //3.asciiplot  二维码图
   
2. 
   
3. ssc install asciiplot
   
4. 
   
5. asciiplot       

复制代码

1. //Venn diagram using venndiag 维恩图之一
   
2. 
   
3. ssc install venndiag
   
4. 
   
5. clear
   
6. input hayfever        eczema        asthma        freq
   
7. 1        0        0        31088
   
8. 1        1        0        9863
   
9. 0        1        0        43522
   
10. 0        1        1        9258
    
11. 0        0        1        35299
    
12. 1        0        1        11024
    
13. 1        1        1        6200
    
14. 0        0        0        345262
    
15. end
    
16. 
    
17. list
    
18. 
    
19. expand freq
    
20. 
    
21. venndiag asthma eczema hayfever

复制代码