/*
	r_t <#df> <varname>

	generates t(#df) distributed random variable
*/

program define r_t
version 2.1
	quietly {
		r_chi2 %_1 _work
		gen %_2 = invnorm(uniform())/sqrt(_work/%_1)
		drop _work
	}
end