Stata基础操作：绘图（三）

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

Stata 12 Graphics

Dawn Koffman Office of  Population Research   Princeton University

矩阵关系图、循环语句图.rar (838 Bytes) 本附件包括： 矩阵关系图、循环语句图.do 2015-7-30 15:54:00 上传

输入数据，修改标签

1. clear  
   
2. input  str14   country    tvhome birth5years idealnum age1stbirth school agemarriage
   
3.             bangladesh     33.9     .6        2.3         17.7     4.5    15.2
   
4.             bolivia        74.8     .8        2.6         20.3     7.6    20.1
   
5.             colombia       93.4     .4        2.4         20.7     8.6    20.3
   
6.             dr             83       .5        3.2         19.9     8.6    18.3
   
7.             egypt          95.9     .7        2.9         21.3     7.3    19.7
   
8.             haiti          27.9     .8        3.2         20.7     4.3    19.4
   
9.             india          47.9     .6        2.4         19.2     4.3    17.1
   
10.             indonesia      74       .5        2.8         20.7     7.5    19.3
    
11.             morocco        66.6     .7        3.3         21.4     2.7    19.6
    
12.             nepal          50.7     .6        2.2         19.6     3.6    17.4
    
13.             pakistan       57.8     .9        4.1         20.5     2.8    18.4
    
14.             peru           78.1     .6        2.5         21.1     8.8    20.6
    
15. end

复制代码

1. label variable tvhome "TV at home (%)"
   
2. label variable school "years of school"
   
3. label variable birth5years "births in last 5 years"
   
4. label variable idealnum "ideal number of children"
   
5. label variable age1stbirth "age at first birth"
   
6. label variable school "years of school"
   
7. label variable agemarriage "age at first marriage"

复制代码

变量标签

矩阵关系图：

基本格式 graph matrix variable1 variable2 variable3, graph_options

Graph options:

1. /*        
   
2. half                 展示一半，分析对比较易
   
3. by(groupingvariable) 使用分组变量，注意先排序
   
4. jitter(#)            大数据集用     
   
5. 技巧：msymbol(p)     数据集较大时 point
   
6. msize(choice)        标记符号 [可选vtiny,tiny,small,medsmall,medium,medlarge,huge,vhuge,large,vlarge,ehuge]
   
7. */

复制代码

1. #delimit ;
   
2. local note  "Source: Most recent DHS std survey: bangladesh bolivia colombia dr egypt
   
3. haiti india indonesia morocco nepal pakistan peru, as of 5/2013";
   
4. graph matrix school tvhome idealnum agemarriage age1stbirth birth5years,
   
5. half note("`note'", size(vsmall));
   
6. #delimit cr

复制代码

关闭half,显示另一半

---

1. #delimit ;
   
2. local note  "Source: Most recent DHS std survey: bangladesh bolivia colombia dr egypt
   
3. haiti india indonesia morocco nepal pakistan peru, as of 5/2013";
   
4. graph matrix school tvhome idealnum agemarriage age1stbirth birth5years,
   
5. note("`note'", size(vsmall));
   
6. #delimit cr

复制代码

使用循环语句画图：

1. #delimit ;
   
2. local sample MARRIED;
   
3. local note  "Source: Most recent DHS standard survey, as of 5/2013"; foreach x in school tvhome {;
   
4. foreach y in idealnum agemarriage age1stbirth birth5years {;
   
5. twoway scatter `y' `x', mlabel(country) mlabsize(large) ylabel(, angle(0)) note("`note'") name("`x'_`y'", replace); graph export `x'_`y'_`sample'_mostrecent.emf, replace;
   
6. };
   
7. };
   
8. #delimit cr

复制代码

将多个直方图同时保存到一张图的方法：

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. }

复制代码

3.核密度图

1. clear all
   
2. sysuse auto
   
3. kdensity length, bw(20)

复制代码

叠加核密度图方法一：（https://bbs.pinggu.org/thread-3606957-1-1.html）感谢蓝色版主

1. sysuse auto,clear
   
2.      kdensity price  if foreign==0 ,  addplot(   ///
   
3.      kdensity price  if foreign==1   )               ///
   
4.      legend(label(1 "国内") label(2 "国外") col(2))

复制代码

叠加核密度图方法二：

1. sysuse auto,clear
   
2. 
   
3.        tw     (kdensity price  if foreign==0)          ///
   
4.                 (kdensity price  if foreign==1) ,       ///  
   
5.                  legend(label(1 "国内") label(2 "国外") col(2))

复制代码

面积图：注：只有在线图下的面积部分具有理论及分析意义时使用。

1. sysuse auto, clear        
   
2. summarize price, mean
   
3. local mean = r(mean)
   
4. kdensity price, gen(x h)
   
5. line h x, || area h x if  x < `mean'

复制代码

1. sysuse gnp96, clear
   
2. #delimit ;
   
3. 
   
4. twoway area d.gnp96 date, xlabel(36(8)164, angle(90))
   
5.   ylabel(-100(50)200, angle(0))
   
6.   ytitle("Billions of 1996 Dollars")
   
7.   xtitle("")
   
8.   subtitle("Change in U.S. GNP", position(11))
   
9.   note("Source: U.S. Department of Commerce, Bureau of Economic Analysis")
   
10.   ;
    
11. #delimit cr

复制代码

//左右同时浮现且保持标题简洁的操作技巧：
//添加xtitle("")

1. sysuse auto,clear
   
2. sort foreign
   
3. set scheme lean1
   
4. histogram rep78, discrete by(foreign) percent   ///
   
5. addlabels xlabel(1 "1" 2 "2" 3 "3" 4 "4" 5 "5") ///
   
6. gap(25) title("Repair Record in 1978") xtitle("")

复制代码

在下端回避Repair Record in 1978

扫码加我 拉你入群

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

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

分享1 收藏21 回帖

关键词：Stata基础 Stata 基础操作 tata Bangladesh University country school 标签

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

楼主，这里第一句话的时候申明str14的意义在哪里

还有这个核密度图里面第一行clear all与正常的clear的差别大不大？

玄一无相 发表于 2015-10-2 21:02
楼主，这里第一句话的时候申明str14的意义在哪里

您好，对于字符型变量在录入时往往需要先定义变量长度，如这里的country：变量值最长为14个字符，其余数值型变量按照默认或者profile.do里面预设的存储格式录入。

玄一无相 发表于 2015-10-2 21:14
还有这个核密度图里面第一行clear all与正常的clear的差别大不大？

您好，在绘图时，clear all在clear清除内存数据和标签的基础上，同时会将现有绘制的图形窗口关闭。

1. //twoway function_UCLA Graph library example
   
2. * 将原normden修正为normalden
   
3. 
   
4. #delimit ;
   
5. twoway
   
6.   function y=normalden(x), range(-4 -1.96) bcolor(gs12) recast(area) ||
   
7.   function y=normalden(x), range(1.96 4) bcolor(gs12) recast(area) ||
   
8.   function y=normalden(x), range(-4 4) clstyle(foreground) ||,
   
9.   plotregion(style(none))
   
10.   yscale(off) xscale(noline)
    
11.   legend(off)
    
12.   xlabel(-4 "-4 sd" -3 "-3 sd" -2 "-2 sd" -1 "-1 sd" 0 "mean"
    
13.   1 "1 sd"   2 "2 sd"   3 "3 sd"   4 "4 sd", grid gmin gmax)
    
14.   xtitle("")
    
15.   ;
    
16. #delimit cr

复制代码

1. sysuse lifeexp, clear
   
2. gen lgnppc = ln(gnppc)
   
3. gr matrix popgr lexp lgnp safe

复制代码

1. sysuse lifeexp, clear
   
2. generate lgnppc = ln(gnppc)
   
3. graph matrix popgr lgnp safe lexp, half

复制代码

1. sysuse lifeexp, clear
   
2. generate lgnppc = ln(gnppc)
   
3. graph matrix popgr lgnp safe lexp, maxes(ylab(#4, grid) xlab(#4, grid))

复制代码