Anαlisis: Generando números aleatorios de un modelo lineal en SAS

febrero 02, 2013

Simulando el siguiente modelo lineal: y = β₀ + β₁x + ε. En donde ε ~N(0,2²). Asumiendo x~N(0,1²), y β₁ = 2

/* Simulando números aleatorios de un modelo de regresión simple*/

data simu;

call streaminit(4); * semilla para numero aleatorio;

do n = 1 to 22;

x = rand('normal', 0, 1);

e = rand('normal', 0, 2); * e ~ U[0,1.5] ;

y = 0.5 + 2*x + e;

output;

end;

run;

ods html;

proc print data=simu;

var x y;

title 'Simulando numeros aleatorios';

run;

* Estadística descriptiva;

proc means;

var x y;

run;

* Trazando grafico;

symbol1 v=dot c=blue h=0.4;

proc gplot data = simu;

plot y * x ;

run; quit;

* Histograma ;

proc univariate normal plot;

var y ;

histogram / normal endpoints=-10 to 10 by 1;

title 'Histograma';

run;

ods html close;

SAS output: