tomhanks 发表于 2009-12-6 00:24
.ivreg lnincome2 age agesqr edu train workexpr workexprsqr indus1-indus5 sex gtshare (gift juhuifrqc=frdgov frdqs married1)
.predict e,residual
.reg e age agesqr edu train workexpr workexprsqr indus1-indus5 sex gtshare frdgov frdqs married1
.scalar sargan=(e(N)-16)*e(r2)
.display sargan
.68169887
.display chi2tail(3-2,sargan)
.40900234
*之后,以下会给出相同的检验结果:
estat over
loc n=wordcount("`e(insts)'")-e(df_m)
predict e,r
reg e age agesqr edu train workexpr workexprsqr indus1-indus5 sex gtshare frdgov frdqs married1
sca sargan=e(N)*e(r2)
di "Sargan (score) chi2(`n')="sargan
di "p="chi2tail(`n',sargan)


雷达卡

京公网安备 11010802022788号







