自用Stata绘图（2）-经管之家官网！

自用Stata绘图（2）

1. /*1.对数坐标散点图*/
   
2. sysuse lifeexp, clear
   
3. gen gnp000 = gnppc/1000
   
4. label var gnp000 "GNP per capita, thousands of dollars"
   
5. scatter lexp gnp000, xsca(log) ///
   
6. xlabel(.5 2.5 10(10)40, grid)

1. /*2.泡泡图加权散点图*/
   
2. sysuse census, clear
   
3. gen drate = divorce / pop18p
   
4. label var drate "Divorce rate"
   
5. scatter drate medage [w=pop18p] if state!="Nevada", msymbol(Oh) ///
   
6. note("State data excluding Nevada" ///
   
7. "Area of symbol proportional to state's population aged 18+")

1. /*3.带置信区间的直线回归图*/
   
2. sysuse auto, clear
   
3. quietly regress mpg weight
   
4. predict hat
   
5. predict stf, stdf
   
6. gen lo = hat - 1.96*stf
   
7. gen hi = hat + 1.96*stf
   
8. scatter mpg weight || line hat lo hi weight, pstyle(p2 p3 p3) sort

1. /*4.期望寿命折线图*/
   
2. sysuse uslifeexp, clear
   
3. gen diff = le_wm - le_bm
   
4. label var diff "Difference"
   
5. #delimit ;
   
6. line le_wm year, yaxis(1 2) xaxis(1 2) ||
   
7. line le_bm year || line diff year ||
   
8. lfit diff year ||, ylabel(0(5)20, axis(2) gmin angle(horizontal))
   
9. ylabel(0 20(10)80, gmax angle(horizontal))
   
10. ytitle("", axis(2)) xlabel(1918, axis(2)) xtitle("", axis(2))
    
11. ytitle("Life expectancy at birth (years)") title("White and black life expectancy")
    
12. subtitle("USA, 1900-1999")
    
13. note("Source: National Vital Statistics, Vol 50, No. 6" "(1918 dip caused by 1918 Influenza Pandemic)")
    
14. legend(label(1 "White males") label(2 "Black males")) legend(col(1) pos(3)) ;
    
15. #delimit cr

1. /*5.带置信区间的散点图*/
   
2. sysuse auto, clear
   
3. quietly regress mpg weight
   
4. predict hat
   
5. predict s, stdf
   
6. gen low = hat - 1.96*s
   
7. gen hi = hat + 1.96*s
   
8. #delimit ;
   
9. twoway rarea low hi weight,
   
10. sort bcolor(gs14) || scatter mpg weight ;
    
11. #delimit cr

1. /*6.折线、散点复合图*/
   
2. sysuse sp500, clear
   
3. #delimit ;
   
4. twoway line close date,
   
5. yaxis(1) || bar change date, yaxis(2) || in 1/52,
   
6. ysca(axis(1) r(1000 1400)) ylab(1200(50)1400,
   
7. axis(1)) ysca(axis(2) r(-50 300)) ylab(-50 0 50, axis(2))
   
8. ytick(-50(25)50, axis(2) grid) legend(off)
   
9. title("S&P 500") subtitle("January - March 2001")
   
10. note("Source: Yahoo!Finance and Commodity Systems, Inc.")
    
11. yline(1150, axis(1) lstyle(foreground)) ;
    
12. #delimit cr

1. /*7.人口金字塔*/
   
2. sysuse pop2000, clear
   
3. replace maletotal = -maletotal/1e+6
   
4. replace femtotal = femtotal/1e+6
   
5. gen zero = 0
   
6. #delimit ;
   
7. twoway bar maletotal agegrp, horizontal xvarlab(Males) ||
   
8. bar femtotal agegrp, horizontal xvarlab(Females) ||
   
9. sc agegrp zero , mlabel(agegrp) mlabcolor(black) msymbol(i) || ,
   
10. xtitle("Population in millions") ytitle("") plotregion(style(none))
    
11. ysca(noline) ylabel(none) xsca(noline titlegap(-3.5))
    
12. xlabel(-12 "12" -10 "10" -8 "8" -6 "6" -4 "4" 4(2)12 ,
    
13. tlength(0) grid gmin gmax) legend(label(1 Males) label(2 Females))
    
14. legend(order(1 2))
    
15. title("US Male and Female Population by Age, 2000")
    
16. note("Source: U.S. Census Bureau, Census 2000, Tables 1, 2 and 3") ;
    
17. scheme(s1mono)
    
18. #delimit cr

1. /*8.折线穗式图*/
   
2. sysuse sp500, clear
   
3. #delimit ;
   
4. twoway line close date,
   
5. yaxis(1) || spike change date, yaxis(2) ||
   
6. , yscale(axis(1) r(700 1400))
   
7. ylabel(1000(100)1400, axis(1))
   
8. yscale(axis(2) r(-50 300))
   
9. ylabel(-50 0 50, axis(2))
   
10. ytick(-50(25)50, axis(2) grid)
    
11. legend(off)
    
12. title("S&P 500")
    
13. subtitle("January - December 2001")
    
14. note("Source: Yahoo!Finance and Commodity Systems, Inc.")
    
15. yline(950, axis(1) lstyle(foreground)) ;
    
16. #delimit cr

1. /*9合并图像*/
   
2. set scheme s1color
   
3. u cumbrian_cirques, clear
   
4. scatter width length , name(g1_1, replace)
   
5. scatter width length , ms(Oh) name(g1_2, replace)
   
6. graph combine g1_1 g1_2, saving(g1, asis replace)

1. /*10.抖动对比*/
   
2. set scheme s1color
   
3. u cumbrian_cirques, clear
   
4. scatter maxgrad grade, ms(Oh) yla(40(10)80, format(%4.0g)) xla(, valuelabel) xsc(r(0.8 5.2)) name(g2_1, replace)
   
5. scatter maxgrad grade, ms(Oh) yla(40(10)80, format(%4.0g)) xla(, valuelabel) xsc(r(0.8 5.2)) jitter(1) name(g2_2, replace)
   
6. graph combine g2_1 g2_2, saving(g2, replace asis)

1. /*11.点图*/
   
2. set scheme s1color
   
3. u cumbrian_cirques, clear
   
4. 
   
5. dotplot maxgrad, over(grade) saving(g3, replace asis) ms(Oh) yla(, format(%4.0g))

1. /*12.散点图*/
   
2. set scheme s1color
   
3. u cumbrian_cirques, clear
   
4. 
   
5. separate length, by(grade) veryshortlabel
   
6. scatter length? width, xsc(log) ysc(log) ms(O ..) mcolor(gs1 gs4 gs7 gs10 gs13) ///

1. /*13.循环作图*/
   
2. set scheme s1color
   
3. u cumbrian_cirques, clear
   
4. 
   
5. forval i = 1/5 {
   
6. scatter length width if grade != `i', xsc(log) ysc(log) xla(200 500 1000 2000) ///
   
7. yla(200 500 1000 2000) ms(oh) mc(blue*0.2) || scatter length width if grade ///
   
8. == `i', xsc(log) ysc(log) xla(200 500 1000 2000) yla(200 500 1000 2000) ///
   
9. ms(O) mc(blue) yla(, ang(h)) subtitle("`: label (grade) `i''", place(w)) ///
   
10. legend(off) name(g11_`i', replace)
    
11. }
    
12. 
    
13. scatter length width, xsc(log) ysc(log) xla(200 500 1000 2000) yla(200 500 1000 2000) ///
    
14. ms(O) mc(blue) yla(, ang(h)) name(g11_6, replace) subtitle("all", place(w))
    
15. 
    
16. graph combine g11_1 g11_2 g11_3 g11_4 g11_5 g11_6, saving(g11, asis replace) imargin(small)