ods pdf file="Ex11_1Gujarati.pdf";

options nocenter linesize=75 pagesize=66;
data gujarati;
   input y x;
datalines;
3396 9355 
3787 8584 
4013 7962 
4104 8275 
4146 8389 
4241 9418 
4388 9795 
4538 10281 
4843 11750 
;

/*  
Exemple du test de Park
*/

proc reg data = gujarati;
  model y =  x;
  output out=datares r=residus p=predites;
run;
 
proc print data=datares;
run;

data datares2;
  set datares;
  res2 = log(residus**2);
  logx = log(x);
  absres = abs(residus);
run;
  
/* Regression de log(residus au carre) sur log(x) */
proc reg data=datares2;
  model res2 =  logx;
  title 'Information pour le test de Park';
run;


/*  
Exemple du test de Glejser
*/
proc reg data=datares2;
  model absres =  x;
  title 'Information pour le test de Glejser';
run;




ods pdf close;