通过Stata程序⽣成⼀组随机数，然后将其排序再⽣成年份变量从1900——2000年，其中数据前55个数据是⼀个规律后56个数据则进⾏变换⽣成⽅式，以此模拟在1955年发⽣⼀次冲击并作为实验组，最后将数据按year排序保存到D:\爬⾍俱乐
部/temp,replace
拟合直线
clear
set obs 101
gen treat = cond(_n<=55,3*runiform(20,25),4*runiform(20,30))
sort treat
gen a=_n
gen year =_n+1899
sort year
save D:\爬⾍俱乐部/temp,replace
模拟⼀组数据为对照组⽽且对照组的数据趋势与实验组前55个数据趋势相同，以便对两组数据进⾏⽐较
clear
set obs 101
gen control=3*runiform(15,25)
sort control
gen a=_n
gen year =_n+1899
sort year
通过merge语句与实验组拼接
merge 1:1 year using D:\爬⾍俱乐部/temp
drop _merge
构造diff观察两组数据的差异
gen diff=treat-control
第⼆步：绘制图形
twoway (line control year)( line treatyear)( line diff year)(lfit diff year)
其中Fitted values是diff与year的拟合曲线
如果给图形添加参考线，我们可以使⽤line或者scatter绘图命令，在命令中加⼊yline() 和xline()
第三步：优化图形
针对上述例⼦在1955年加⼀个参考线，使⽤xline()命令
line control year,lwidth( medthick) yaxis(1) xline(1955) xaxis(1) lwidth(medium) ||
line treat year ,lwidth(medthick)||
line diff year,lwidth( medthick)||
lfit diff year,lwidth( medthick)
这样可以更加清楚地看见1955年曲线发⽣的变动
通过上述两个⽅法绘图我们可以轻松将四条曲线绘制在同⼀张图表中
如果您嫌背景⾊不好看，可以通过以下⽅法修改底⾊，⼤家可以⾃⾏尝试：
执⾏完 do ⽂档图⽚显⽰出来后，在图⽚上点右键，然后可以选择编辑图⽚（点 preferences），⾥⾯可以在 scheme ⾥调主题颜⾊，默认的是 s2color。不过调完之后，则需要在下次执⾏ do ⽂档的时候才会⽣效，⽽且是今后⼀直⽣效，要想⽤默认颜⾊需要⾃⼰再改回来。例如： 我们也可以看⼀下scheme的值换为s2manual-mono，图形变为
⼤家可以根据⾃⾝喜好选择Scheme的类型
如果想继续对图⽚进⾏编辑也可以右键Start Graph Editor
如果对图形背景颜⾊进⾏选择应对Color进⾏选择，图形边界修改则是Margin，如果想修改更多则选择more，如图：
对曲线进⾏修改只需双击此条曲线
如果读者想添加辅助线可以选择Reference lines
aaplot回归拟合曲线图的奇妙
强⼤的aaplot命令可以画线性和⼆次拟合以及两者的叠加。并给出了⽅程表达式、R^2（显⽰在图的顶部）、观测值的个数和RMSE（显⽰在图的底部）。 第⼀步：模拟数据
clear
set obs 100
gen a=_n
gen x=a*runiform(2,3)
gen y=x^2+x+10+runiform(1,40000)
第⼆步：绘制线性拟合直线图
aaplot y x, name(a1) //a1是所绘制图形的名称
由图可以得出线性拟合的的回归⽅程以及R^2，以及图下⽅的观测值和RMSE均⽅根误差
第三步：绘制⼆次拟合曲线
aaplot y x, quadratic qopts(lc(pink)) name(a2)
lc(pink)是定义⼆次拟合曲线为粉红⾊，⼤家可以根据⾃⼰喜好⾃⾏定义
根据图形我们可以通过判定系数R^2观测⼆次曲线拟合的更好
第四步：线性与⾮线性的⽐较
aaplot y x, both name(a3)
通过使⽤both我们将线性和⾮线性图组合到⼀起，且同时给出两次拟合的回归⽅程
为了更⽅便观察⼩编对命令进⾏了修饰
aaplot y x, both backdrop(lfitci y x,color(gs12)) name(a4)
对线性拟合曲线的背景进⾏定义，可以更加明显区分线性和⾮线性拟合的不同
如果⼤家想进⼀步了解aaplot，help⼀下你就知道！
以上数据均为⼩编随机⽣成的，希望通过使⽤上述命令可以对⼤家处理数据绘图有所帮助。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议