allocate 100


compute beta1=0, beta2=1, sigmax=1,sigmau=2
compute ndraws = 10000




do draw=1,ndraws

compute x = %ran(sigmax)
compute u = %ran(sigmau)

set y = beta1+beta2*x+u


linreg(noprint) y
# constant x

end do draw