libname jason "/home/jaosborn/classes/st512sum08/datafiles";
data cars;
   set jason.cars1993;
   weight2=weight*weight;
run;

proc print data=cars (obs=10);
run;


/*
proc print data=cars;
   var numcyl  engliters  length  width  wheelbase  weight;
run;
*/
proc reg data=cars;
   model citympg = numcyl  engliters  length  width  wheelbase  weight;
   mycomparison: test numcyl=engliters=length=width=0;
   model citympg = numcyl  engliters  length  width  wheelbase  weight/selection=cp;
   model citympg = engliters  length  width  wheelbase  weight weight2/selection=cp;
   model citympg = weight wheelbase;
   output out=myresidz p=p r=r;
   plot citympg*predicted.;
run;
proc glm;
   model citympg = weight wheelbase weight*weight/clparm;
   estimate "Osborne's junker" intercept 1 weight 2957 wheelbase 97.4 weight*weight 8743849;
run;