Note: 助教招聘信息请进入「课程主页」查看。因果推断-内生性 专题 ⌚ 2020.11.12-15 主讲:王存同 (中央财经大学);司继春(上海对外经贸大学)空间计量 专题 ⌚ 2020.12.10-13 主讲:杨海生 (中山大学);范巧 (兰州大学)
编译:张晓明 (中国人民大学);连玉君(中山大学)
Stata 连享会: 知乎 | 简书 | 码云
下载 dofile: 连享会-十幅经典图形绘制.do (可右击另存,亦可在线查看)
1. 散点图:附加密度函数图和拟合线(Scatter Plot with Regression Results)
附加密度函数和拟合曲线的散点图能够更好地看清样本的概率分布,同时能够表示出参数估计的置信区间。
*-文件夹设定
cd D:\
mkdir myfigs
cd D:\myfigs // 后文输出的图形自从存储于此处
*-----------
*- F1 散点图:附加密度函数和拟合曲线图绘制
*-----------
sysuse auto, clear
version 9.2 //绘制这幅图的时候,需要加这条命令,后续版本的stata绘图命令的语法有所改变
keep if foreign
sort weight
gen weight2 = weight^2
regress mpg weight weight2
predict fit
predict se , stdp
#delimit ;
twoway
scatter mpg weight , pstyle(p3) ms(o) ||
fn weight[3] - 1000 * normden(x, `=fit[3]' , `=se[3]') ,
range(`=fit[3] -5' `=fit[3] +5') horiz pstyle(p1) ||
fn `=fit[3]' , range(`=weight[3]' `=weight[3]-1000*normden(0, se[3])')
pstyle(p1) ||
fn weight[17] - 1000 * normden(x, `=fit[17]', `=se[17]') ,
range(`=fit[17]-5' `=fit[17]+5') horiz pstyle(p1) ||
fn `=fit[17]', range(`=weight[17]' `=weight[17]-1000*normden(0, se[17])')
pstyle(p1) ||
fn weight[21] - 1000 * normden(x, `=fit[21]' , `=se[21]') ,
range(`=fit[21] -7' `=fit[21] +7') horiz pstyle(p1) ||
fn `=fit[21]', range(`=weight[21]' `=weight[21]-1000*normden(0, se[21])')
pstyle(p1) ||
line fit weight, clwidth(*2) legend(off)
ytitle(Miles per gallon) xtitle(Weight)
title("Scatter with Regression Line and Confidence Interval Densities"