Stata基础操作：绘图（六）人口金字塔及脸谱图

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

Stata 12 Graphics

Dawn Koffman

Office of  Population Research   Princeton University

人口金字塔条形图绘制

1. sysuse pop2000, clear
   
2. replace maletotal = -maletotal
   
3. #delimit ;
   
4. twoway bar maletotal agegrp, horizontal ||
   
5. bar femtotal agegrp, horizontal
   
6. title("US Male and Female Population by Age, Year 2000")
   
7. note("Source: US Census Bureau, Census 2000, Tables 1, 2 and 3", span);
   
8. #delimit cr

复制代码

人口金字塔（含坐标轴标签）：

1. sysuse pop2000, clear
   
2. replace maletotal = -maletotal
   
3. replace maletotal = maletotal / 1000000
   
4. replace femtotal = femtotal / 1000000
   
5. #delimit ;
   
6. twoway bar maletotal agegrp, horizontal ||
   
7. bar femtotal agegrp, horizontal
   
8. title("US Male and Female Population by Age, Year 2000")
   
9. note("Source: US Census Bureau, Census 2000, Tables 1, 2 and 3", span)
   
10. xtitle("Population in Millions") ytitle("Age Group Number")
    
11. xlabel( -12 "12" -10 "10" -8 "8" -6 "6" -4 "4" 4(2)12)
    
12. ylabel(1(1)17, angle(0))
    
13. legend(order(1 "Male" 2 "Female"));
    
14. #delimit cr

复制代码

人口金字塔（年龄分组标签）

1. sysuse pop2000, clear
   
2. replace maletotal = -maletotal
   
3. replace maletotal = maletotal / 1000000
   
4. replace femtotal = femtotal / 1000000
   
5. gen zero = 0
   
6. #delimit ;
   
7. twoway bar maletotal agegrp, horizontal bfc(gs7) blc(gs7) ||
   
8. bar femtotal agegrp, horizontal bfc(gs11) blc(gs11) ||
   
9. scatter agegrp zero, mlabel(agegrp) mlabcolor(black) msymbol(none)
   
10. title("US Male and Female Population by Age, Year 2000")
    
11. note("Source: US Census Bureau, Census 2000, Tables 1, 2 and 3", span)
    
12. xtitle("Population in Millions") ytitle("Age Group Number")
    
13. ytitle("") yscale(noline) ylabel(none)
    
14. xlabel( -12 "12" -10 "10" -8 "8" -6 "6" -4 "4" 4(2)12)
    
15. legend(off) text(15 -8 "Male") text(15 8 "Female");
    
16. #delimit cr

复制代码

注：关于人口金字塔的绘制可参阅   演示教程命令集及Pyramid命令修改

/*Chernoff 脸谱图*/

1. * Chernoff faces
   
2. * Few cases, each face is a row case.
   
3. 
   
4. use "http://www.princeton.edu/~otorres/chernoff.dta", clear
   
5. *ssc install chernoff    //chernoff移至gr0038下
   
6. findit gr0038    //安装
   
7. chernoff, hdark(gdppc) bdens(trade) nose(unemp) mcurv(polity2) order(gdppc) ilabel(country)

复制代码

/*圆环图*/

1. // Circle using function
   
2. #delimit ;
   
3. twoway (function y = sqrt(1 -(x)^2),
   
4. range(-1 1) lwidth(thick) lcolor(red))
   
5. (function y = -sqrt(1 - (x)^2),
   
6. range(-1 1) lwidth(thick) lcolor(red)), aspect(1);

复制代码

---

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享1 收藏20 回帖

关键词：Stata基础 Stata 基础操作 tata 脸谱图 University horizontal 金字塔 Office 坐标轴

好贴，谢谢分享，再接再厉。

顶起{:2_38:}

mark，学习了。

great~学习了

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

复制代码

1. //修改后
   
2. // Circle using function
   
3. #delimit ;
   
4. twoway (function y = sqrt(1 -(x)^2),
   
5. range(-1 1) lwidth(thick) lcolor(red))
   
6. (function y = -sqrt(1 - (x)^2),
   
7. range(-1 1) lwidth(thick) lcolor(red)), aspect(1)
   
8. legend(off) yla(, ang(h)) ytitle(, orient(horiz));

复制代码

关闭图例，调整y标签方向

好！！！！

好贴，再次mark 一下。

1. sysuse nlsw88, clear
   
2. #delimit ;
   
3. graph hbar wage, over(ind, sort(1)) over(collgrad)
   
4.   title("Average hourly wage, 1988, women aged 34-46", span)
   
5.   subtitle(" ")
   
6.   note("Source:  1988 data from NLS, U.S. Dept. of Labor, Bureau of Labor Statistics", span)
   
7.   ysize(7) ;
   
8. #delimit cr

复制代码