n _ 1000
mu.x_10
mu.y_10
mu_c(mu.x,mu.y)
sigma_c(2,2)
rho_.8

dat_rmvnorm(n,mean=mu,sd=sigma,rho=rho)
x_dat[,1]
y_dat[,2]

plot(x,y)
line.fit_lm(y~x)
abline(line.fit)
res_line.fit$residuals
qqnorm(res)
plot(x,res)
summary(line.fit)