program hausmantest,eclass
tempname b bfe bre
xtreg fatal beertax spircons unrate perinck, fe r
matrix `bfe'=e(b)
xtreg fatal beertax spircons unrate perinck, re r
matrix `bre'=e(b)
matrix `b'=`bfe'-`bre'
ereturn post`b'
end
bootstrap _b,reps(400) nodots nowarn:hausmantest
结果如下