cap log close
clear all
eststo clear
log using "PS1.txt", text replace

*PRE*

***********Directory Settings************
global path "/Users/levivanboekel/Desktop/ku/øko/ugesedler"
global outtex "/Users/levivanboekel/Desktop/ku/øko/ugesedler/a1.tex"

******************************************

***********Data************
cd "/Users/levivanboekel/Desktop/ku/øko/ugesedler"
******************************************

clear all

global numobs = 50

global beta_0 =1
global beta_1=2
global beta_2=0
global rho = 0.5


*DEFINE PROGRAM THAT SPECIFIES THE DGP
program olsdata, rclass
	drop _all
	
	*SET NUMBER OF CURRENT OBSERVATIONS 
	set obs $numobs
	
	*DATA GENERATING PROCESS
	generate rho = $rho
	generate beta_0 = $beta_0
	generate beta_1 = $beta_1
	generate beta_2 = $beta_2
	generate x1 = 25 + 5*rnormal()
	generate x2star = 10 + 20*runiform()
	generate x2 = $rho *x1 + x2star
	generate u = -50 + 100*runiform()
	generate y = $beta_0 + $beta_1 * x1 + $beta_2 * x2 + u
	
	*CALCULATE OLS ESTIMATES AND SAVE RESULTS 
	regress y x1
	return scalar betahat_slr = _b[x1]
	return scalar se_slr    = _se[x1]
	regress y x1 x2
	return scalar betahat_mlr = _b[x1]
	return scalar se_mlr   = _se[x1]
end

olsdata

simulate betahat_slr=r(betahat_slr) se_slr=r(se_slr) ///
         betahat_mlr=r(betahat_mlr) se_mlr=r(se_mlr), ///
         seed(1479) reps(1000) nodots saving("sim_results.dta", replace): olsdata


twoway (hist betahat_slr, color(navy)) (hist betahat_mlr, color(gs10)),



use "sim_results.dta", clear


estpost summarize betahat_slr se_slr betahat_mlr se_mlr

esttab using "${outtex}", replace ///
    cells("mean sd min max count") ///
    nonumber nomtitle noobs  ///
    title("Resultater af MC for n=$numobs")