program define report /* level dof truth coef stderr */
	version 3.0
	local lvl = `1'
	local dof = `2'
	local truth = `3'
	local b "`4'"
	local se "`5'"
	di _n "report for `b' `se' true value `3'
	di "dof = " `dof'
	quietly {
		capture drop lbt
		capture drop ubt
		capture drop lbz
		capture drop ubz
		local t = invt(`dof',`lvl'/100)
		gen lbt = `b'-`t'*`se'
		gen ubt = `b'+`t'*`se'
		quietly count if `b'!=.
		local nc=_result(1)
		quietly count if `b'!=. & `truth'>=lbt & `truth'<=ubt
		local nir = _result(1)
		noi di "`nir'/`nc' = " `nir'/`nc' " in `lvl'% t-range"
		local z = invnorm(1-(1-`lvl'/100)/2)
		gen lbz = `b'-`z'*`se'
		gen ubz = `b'+`z'*`se'
		quietly count if `b'!=. & `truth'>=lbz & `truth'<=ubz
		local nir = _result(1)
		noi di "`nir'/`nc' = " `nir'/`nc' " in `lvl'% Z-range"
	}
end