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
ivregress 2sls lnincome2
age agesqr edu train workexpr workexprsqr indus1-indus5 sex gtshare (gift juhuifrqc
=frdgov frdqs married1)
*之后,以下会给出相同的检验结果:
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
)