clear

capture log close
log using unbiased.txt, text replace
di "unbiased.txt"
version 7

drop _all

set obs 20

gen x = invnorm(uniform())
gen u = invnorm(uniform())
gen y = 1 + 2*x + u


regress y x


capture program drop regwrs
program define regwrs
if "`1'" == "?" {
	global S_1 "b_const se_const b_slope se_slope"
	exit
	}
	drop _all
	set obs 20
	gen x = invnorm(uniform())
	gen u = invnorm(uniform())
	gen y = 1 + 2*x + u
	reg y x
	post `1' (_b[_cons]) (_se[_cons]) (_b[x]) (_se[x])
end

simul regwrs, reps(1000)
summarize

log close