set mem 300m

set obs 50000

gen z_1=invnorm(uniform())
gen z_2=invnorm(uniform())
gen z_3=invnorm(uniform())
gen z_4=invnorm(uniform())
gen z_5=invnorm(uniform())
gen z_6=invnorm(uniform())
gen z_7=invnorm(uniform())
gen z_8=invnorm(uniform())
gen z_9=invnorm(uniform())
gen z_10=invnorm(uniform())
gen z_11=invnorm(uniform())
gen z_12=invnorm(uniform())
gen z_13=invnorm(uniform())
gen z_14=invnorm(uniform())
gen z_15=invnorm(uniform())
gen z_16=invnorm(uniform())
gen z_17=invnorm(uniform())
gen z_18=invnorm(uniform())
gen z_19=invnorm(uniform())
gen z_20=invnorm(uniform())
gen z_21=invnorm(uniform())
gen z_22=invnorm(uniform())
gen z_23=invnorm(uniform())
gen z_24=invnorm(uniform())
gen z_25=invnorm(uniform())

gen F_1=(z_1^2)/(z_1^2+z_2^2+z_3^2+z_4^2+z_5^2+z_6^2+z_7^2+z_8^2+z_9^2+z_10^2+z_12^2+z_13^2+z_14^2+z_15^2+z_16^2+z_17^2+z_18^2+z_19^2+z_20^2+z_21^2+z_22^2+z_23^2+z_24^2+z_25^2)