stata怎么画分类图_Stata数据可视化: 十幅精美图形的绘制

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"