libname bus 'C:\Documents and Settings\hraissi\Bureau\toto';
data bus;
 input NB CAR REV;
cards;
2.073	0.88	17.293	
2.136	1.03	17.768	
1.878	0.91	17.823	
0.937	0.91	15.163
7.343	0.97	17.48	 
0.837	0.88	15.329	
1.648	0.91	16.141	
0.739	0.89	15.326	
1.07	0.89	17.115	
0.274	0.89	17.117	
0.312	0.87	16.127	
1.879	0.94	17.242	
1.941	0.99	17.34	
2.317	0.87	15.108	
0.471	0.93	15.809	
0.594	0.79	16.321	
7.632	0.93	18.027	
0.51	0.93	18.023	
0.63	0.93	12.349	
1.65	1.03	17.886	
1.618	0.86	16.537	
2.009	0.96	13.019	
1.562	0.96	13.019	
1.139	0.88	13.13		
13.103	1.00	20.513   
run;
proc print data=bus;
run;

proc means data=bus;
run;
proc univariate data=bus normal;
var REV NB;
run;
proc chart data=bus;
vbar NB /midpoints=1 to 13 by 2;
vbar CAR /midpoints=0.85 to 1.05 by 0.02;
vbar REV /midpoints=14 to 20 by 0.5;
run;

proc reg data=bus;
model NB=CAR REV /I VIF TOL COLLIN COLLINOINT;
output out=bus.B R=res;
run;
proc univariate data=bus.B normal;
var res;
qqplot;
run;

DATA busun;
SET bus;
carun = CAR*100;
RUN;
proc reg data=busun;
model NB=carun REV /I VIF TOL COLLIN COLLINOINT;
run;

proc reg data=bus gout=bus.graph;
model NB= REV/ r influence;
plot NB*REV='*';
plot R.*P.;
output out=bus.res residual=ei;
run;

data busprime;
 input NBp CARp REVp;
cards;
2.073	0.88	17.293	
2.136	1.03	17.768	
1.878	0.91	17.823	
0.937	0.91	15.163  
7.343	0.97	17.48	
0.837	0.88	15.329	
1.648	0.91	16.141	
0.739	0.89	15.326	
1.07	0.89	17.115	
0.274	0.89	17.117	
0.312	0.87	16.127	
1.879	0.94	17.242	
1.941	0.99	17.34	
2.317	0.87	15.108	
0.471	0.93	15.809	
0.594	0.79	16.321	
7.632	0.93	18.027	
0.51	0.93	18.023	
0.63	0.93	12.349	
1.65	1.03	17.886	
1.618	0.86	16.537	
2.009	0.96	13.019	
1.562	0.96	13.019	
1.139	0.88	13.13	
run;

proc reg data=busprime gout=bus.graphprime;
model NBp= REVp/ r influence;
plot NBp*REVp='*';
plot R.*P.;
output out=bus.res residual=vi;
run;
data bustrois;
 input NBt CARt REVt;
cards;
2.073	0.88	17.293	
2.136	1.03	17.768	
1.878	0.91	17.823	
0.937	0.91	15.163
0.837	0.88	15.329	
1.648	0.91	16.141	
0.739	0.89	15.326	
1.07	0.89	17.115	
0.274	0.89	17.117	
0.312	0.87	16.127	
1.879	0.94	17.242	
1.941	0.99	17.34	
2.317	0.87	15.108	
0.471	0.93	15.809	
0.594	0.79	16.321	
0.51	0.93	18.023	
0.63	0.93	12.349	
1.65	1.03	17.886	
1.618	0.86	16.537	
2.009	0.96	13.019	
1.562	0.96	13.019	
1.139	0.88	13.13	
run;

proc reg data=bustrois gout=bus.graphtrois;
model NBt= REVt/ r influence;
plot NBt*REVt='*';
plot R.*P.;
output out=bus.res residual=wi;
run;

proc univariate normal  gout=bus.graphtr;
var wi;
qqplot wi;
run;

data busquat;
   set bus end=eof;
   output;
   if eof then
      do;
         REV=16.5;
         NB=.;
         output;
      end;
run;
proc reg data=busquat;
model NB=REV/ CLI CLM;
run;
quit;