/*
	r_cnrml <varname>

	generates contaminated Normal
		.95	Normal(0,1) 
		.05	Normal(0,10)
*/

program define r_cnrml
version 2.1
	gen %_1 = invnorm(uniform())
	quietly replace %_1 = 10*%_1 if uniform()<=.05
end