clear
log using awd01                            /* rozpoczecie loga */ 
set obs 1000 			           /* ustawienie liczby obserwacji */
gen x=invnorm(uniform())	           /* wylosowanie zmiennej z rozkladu normalnego*/
gen e=uniform() 			   /* wygenerowanie zmiennej z rozkladu U[0,1] */
gen x1=2*x+10*e+6
gen f=4*invnorm(uniform())                 /* wygenerowanie zmiennej z rozkladu U[0,4] */
gen x2=0.4*x-f
gen g=5000*uniform()
gen d1=0                 		   /* wygenerowanie zmiennej rownej 0 */
recode d1 0=1 if mod(round(100*g),2)==1	   /* wygenerowanie zmiennej 0-1 */
gen d2=mod(round(170*g),10)        	   /* zmienna o 10 kategoriach */				

label var x  "Zmienna X"                   /* nadanie nazw zmiennym */
label var x1 "Zmienna X_1"
label var x2 "Zmienna X_2"
label var d1 "Zmienna 0-1" 
label var d2 "Zmienna jakosciowa"

describe
summarize 
histogram x			 	   /* wykres slupkowy */
histogram x1, bin(50)
histogram d1
histogram d1, discrete
histogram d2
histogram d2, discrete

twoway (scatter x1 x) (scatter x2 x)
twoway (scatter x1 x) (lfit x1 x)
twoway (scatter x2 x) (lfit x2 x)

reg x x1 x2 d1 d2
matrix list e(b)			   /* wyswietlenie wektora "beta" */
matrix list e(V)     	  	           /* macierz wariancji-kowariancji estymatorow */
ereturn list

sort x1
twoway (line x x1, clcolor("0 128 255") clpat(solid) clwidth(medthick)) , ytitle(Stopa zwrotu) yscale(range(-5 10)) ylabel(-6(2)10) xtitle(Długość inwestycji) xlabel(0(1)25, valuelabel) xmlabel(, noticks nolabels valuelabel) title(Zwrot z inwestycji) legend(off)

predict res,resid                         /* generowanie reszt*/
avplots			  	          /* wartosci rzeczywiste vs dopasowane */
symplot res			  	  /* wykres symetrycznosci reszt*/
pnorm res			 	  /* wykres normalnosci reszt*/	
graph matrix x x1 x2 d1 d2
 
sktest res			  	  /* test normalnosc */
vif                       	          /* test wspoliniowosci */
ovtest, rhs                               /* test RESET */ 
hettest, rhs                              /* test homoscedastycznosci */ 
whitetst
net search whitetst                       /* przeszukiwanie internetu w celu znalezienia komendy */
whitetst


log close