根据陈强《高级计量经济学及stata应用》P163及论坛帖https://bbs.pinggu.org/thread-2317592-1-1.html
在stata中编写的代码：
quietly xtreg y1 lnm3 w2 w4 p r lniv lngdp lnpop,re
scalar theta=e(theta)
global yandxforhausman y1 lnm3 w2 w4 p r lniv lngdp lnpop
sort province
foreach x of varlist $yandxforhausman {
by province:egen mean`x'=mean(`x')
gen md`x'=`x'-mean`x'
gen red`x'=`x'-theta*mean`x'
}
quietly xtreg redy1 redlnm3 redw2 redw4 redp redr redlniv redlngdp redlnpop /// mdlnm3 mdw2 mdw4 mdp mdr mdlniv mdlngdp mdlnpop, vce(cluster province)
test mdlnm3 mdw2 mdw4 mdp mdr mdlniv mdlngdp mdlnpop
est table OLS FE_robust RE,b se
exit

但是出现以下问题：
. test mdlnm3 mdw2 mdw4 mdp mdr mdlniv mdlngdp mdlnpop
mdlnm3 not found
r(111);

去掉lnm3也不行
. test mdw2 mdw4 mdp mdr mdlniv mdlngdp mdlnpop
mdw2 not found
r(111);

但是数据里是有这些变量也可以跑reg的，请教大家