STATA画图命令（二）

曲线标绘图和连线标绘图

[twoway] line varlist [if] [in] [,options] //绘制曲线标绘图

quietly regress mpg weight  //quietly要求不显示回归结果，简写为qui

predict hat  //得到预测值hat

predict stdf, stdf  //得到预测的标准差

gen lo = hat - 1.96*stdf  //构建置信度为95%的双侧检验t值置信区间下限

gen hi= hat + 1.96*stdf  //构建置信度为95%的双侧检验t值置信区间上限

scatter mpg weight || line hat lo hi weight, pstyle (p2 p3 p3) sort  

//要求以hat lo hi 为y轴变量，以weight为x轴变量绘制散点图；

sort要求对weight进行排序，需要强调的是sort必须加上，否则会得到无意义图形


gen diff = le_wm - le_bm  //生成一个新变量diff，其数值是女性的预期寿命减去男性预期寿命

label var diff “Difference”  //为变量设置标签

line le_bm year, yaxis(1 2) xaxis(1 2)  //男性预期寿命随时间变化的标绘图，并绘制两个坐标轴

  || line le_wm year  //女性预期寿命和时间变化曲线标绘图

  || line diff year  //性别预期寿命差值和时间曲线标绘图

  || lfit diff year  //性别预期寿命差值和时间一次拟合线

  ||,

  ylabel(0(5)20,axis(2) gmin angle(horizontal))  //设置右边的坐标轴在0-20之间，步长为5，并在最小值处标注，表示为水平而不是默认的垂直

  ylabel(0 20(10)80 ,axis(2)  gmax angle(horizontal))  //设置左边的坐标轴在0-80之间，其中20-80的步长为10，并在最大值处标注，表示为水平而不是默认的垂直

  ytitle(“”,axis(2))  //省略右边y轴的标题

  xlabel(1980,axis(2)) xtitle(“”,axis(2))  //上侧x轴标注1980位置，且无标题

  ylabel(,axis(2) grid)  //y轴标注网格线

  ytitle(“LEBY”)  //y轴坐标轴标题

  title(“WBLE”)  //总标题设置

  subtitle(“USA, 1900-1999”)  //副标题设置

  note("ASDFGH""ARTHHU")  //下方注释，两行文字

legend(label(1 “WM”) label(2 “BM”))  //第一条图例文字说明改为WM，第二条改为BM

legend(col(1) pos(3))  //图例拍成一列，同时位于图形的三点钟方向

twoway connected varlist [if] [in] [weight] [,scatter_options]

连线标绘图与散点图本质上是一致的，连线标绘图可以与默认情况下的散点图相互转换

拟合标绘图

(1)绘制一次拟合图

twoway lfit yvar xvar [if] [in] [weight] [,options]

[,options]：range(##)计算预期值的x区域；

n(#)计算预期值的点的数目；

atobs在xvar赋值处计算yvar的预期值；

cline改变连线的外形；

axis选择需要的坐标轴

scatter mpg weight || lfit mpg weight

等价于：

reg mpg weight

predict fitted

scatter mpg weight || line fitted weight

*根据分类变量绘制复合图形

scatter mpg weight || lfit mpg weight ||,by(foreign,total row(1))  //根据分类变量foreign绘制所有一次拟合图，并排成一行

(2)绘制二次拟合图

twoway qfit yvar xvar [if] [in] [weight] [,options]  //命令选项同一次拟合图

scatter mpg weight || qfit mpg weight

等价于：

gen tempvar = weight^2

reg mpg weight tempvar

predict fitted

scatter mpg weight || line fitted weight

*根据分类变量绘制复合图形

scatter mpg weight || qfit mpg weight ||,by(foreign,total row(1))  //根据分类变量foreign绘制所有二次拟合图，并排成一行

(3)绘制lowess拟合图

twoway lowess yvar xvar [if] [in] [,options]

[,options]：bwidth设定波段宽度，默认是bwidth(0.8)

mean用移动平均数法修匀，默认是最小二乘法

noweight未加权的修匀，默认为局部加权函数

logit将修匀值做logit变换

adjust调整修匀值的函数

twoway scatter mpg weight

||lfit mpg weight

||lowess mpg weight

*根据分类变量绘制复合图形

scatter mpg weight || lfit mpg weight || lowess mpg weight ||,by(foreign,total row(1))  //根据分类变量foreign绘制所有一次拟合图及修匀拟合图，并排成一行