stata命令(两阶段最小二乘,弱工具变量检验)

webuse nlswork(
两阶最小二乘；弱工具变量检验；collapse；table高级用法；面板数据的设定
help ivregress（实例看怎么用）
ssc install ivreg2
ssc install ranktest
d
ivregress 2sls ln_wage i.race ttl_exp collgrad (union=c_city south not_smsa)内生变量：是否加入公会， 工具变量：后面3个
ivregress 2sls ln_wage i.race ttl_exp collgrad (union=c_city south not_smsa),robust
稳健的方差（怀特）

弱工具变量检验：
reg union c_city south not_smsa
help test
test (c_city=0) (south=0) (not_smsa)
test (c_city=0) (south=0) (not_smsa=0)

F大于10——联合显著

扩展命令reg2没有i.这个命令生成虚拟变量
使用的旧的语法格式：
xi: ivreg2 ln_wage i.race ttl_exp collgrad (union=c_city south not_smsa)
weak identification test的结果：197.042 大于5%那个临界值13.91就可以
xi: ivreg2 ln_wage i.race ttl_exp collgrad (union=c_city south not_smsa),robust

面板数据回归：
软件默认数据——截面数据
个人维度+时间维度组合起来可以确定唯一的观测值
xtset x——截面 t——时间
reg ln_wage L.ln_wage,robust
xtset idcode year——设置成面板数据
xtdes——看数据分布
过去的收入如何影响当前的工资：reg ln_wage L.ln_wage,robust L.ln_wage——上一期的工资
时间序列数据设定——tset（后面加时间变量的信息）
工资的增长率如何受过去工资的影响：reg D.ln_wage L.ln_wage,robust
不同行业人的平均工资：
tabstat ln_wage,by(ind_code)——只有mean
tabstat ln_wage,by(ind_code) stats(mean median sd)——有其他的

tab ind_code south,su(ln_wage)
tab ind_code south,su(ln_wage) nofreq mean——不同地区，不同行业的人的平均工资的列联表（只求均值，不求频数）
help table
table ind_code south,contents(mean ln_wage)
table ind_code south,contents(mean ln_wage) row col
table ind_code south union,contents(mean ln_wage) row col

preserve
画图：
collapse(mean) mln_wage=ln_wage (sd) sdln_wage=ln_wage,by(ind_code)
twoway connected mln_wage ind_code|| connected sdln_wage ind_code
twoway connected mln_wage ind_code|| connected sdln_wage ind_code,yaxis(2)
restore
reg ln_wage ttl_exp
help statsby
preserve
statsby coef=_b[ttl_exp],by(ind_code):reg ln_wage ttl_exp
list——看不同行业前面的系数
不同行业工作经验多一年，工资平均增加多少—— twoway connected coef ind_code

restore
help forvalues
doedit
\行拼接 逗号——列拼接
forvalues i=1/10{
display="`i'"
}
forvalues i=1/12{
reg ln_wage ttl_exp if ind_code=`i'
}


capture mat drop coef
forvalues i=1/12{
reg ln_wage ttl_exp if ind_code==`i' i左边的符号是数字1左边那个，右边的是单引号
mat coef=nullmat(coef)\_b[ttl_exp]————输出结果是行拼接的
}