capture program drop simalpha
program simalpha, rclass
	clear
	set obs 100
	forvalues i=1/5 {
		generate x`i' = rnormal()
		}
	alpha x1-x5, std asis
	return scalar alpha = r(alpha)
end

preserve
simulate alpha=r(alpha), reps(1000): simalpha
summarize alpha
histogram alpha
restore

alpha x1-x5, std asis
*display r(alpha)/(_N-(_N-1)*r(alpha))
*display r(rho)/(_N-(_N-1)*r(rho))
display r(alpha)/(5-(5-1)*r(alpha))
*display r(rho)/(5-(5-1)*r(rho))
preserve
generate n = _n
quietly reshape long x, i(n) j(test)
icc x n
*icc x test
restore