/*
	r_chi2 <#df> <varname>

	generates Chi2(#df) distributed random variable
*/

program define r_chi2
version 2.1
	quietly {
		gen %_2 = invnorm(uniform())^2
		mac def _i = 1 
		while (%_i<%_1) { 
			replace %_2 = %_2 + invnorm(uniform())^2
			mac def _i = %_i + 1 
		}
	}
end