|
capture program drop reschi2
program reschi2,rclass
version 14.0
syntax[,obs(integer 1)]
drop _all
set obs `obs'
tempvar x y e
gen double `x' = rchi2(1)
gen `e' = rnormal(0,3)
gen `y' = 3+2*`x'+rchi2(1)-3+`e'
reg `y' `x'
return scalar b = _b[`x']
return scalar se = _se[`x']
return scalar t = (return(b)-0)/return(se)
return scalar re = abs(return(t)) > invttail(2,0.025)
end
simulate b = r(b) se = r(se) t = r(t) r = r(re),reps(1000) seed(101): reschi2,obs(100)
|