options nodate nocenter ls=80 ps=2000;
libname a 'd:\corinne\cours\desscos\appli2';

proc contents data=a.base;
run;

proc sort data=a.base;
by effectif;
run;

proc print data=a.base(obs=20);
run;

proc freq data=a.base;
table secteur;
run;

proc means data=a.base;
var prod;
run;

proc univariate data=a.base;
var prod;
freq prod;
histogram prod; 
run;

proc gplot data=a.base;
symbol1 v=star i=none;
plot prod*idobs;run;
plot effectif*idobs;run;
plot pub*idobs;run;
plot prod*effectif;run;
plot prod*pub;run;
quit;

proc corr data=a.base;
var prod effectif pub ;
run;

proc reg data=a.base corr simple ;
Reg1: model prod=effectif ;
output out=res p=fit r=res ;
quit;

proc gplot data=res;
symbol1 v=star i=none c=green;
symbol2 v=none i=join c=red;
plot res*idobs=1;
plot prod*effectif=1 fit*effectif=2/overlay;
run;

/* en logarithme */

proc gplot data=a.base;
symbol1 v=star i=none;
plot lprod*idobs;run;
plot leffectif*idobs;run;
plot lprod*leffectif;run;

proc reg data=a.base;
model lprod=leffectif /*lpub */;
output out=ress p=fit r=res ;
run;

proc gplot data=ress;
symbol1 v=star i=none c=green;
symbol2 v=none i=join c=red;
plot lprod*leffectif=1 fit*leffectif=2/ overlay ;
run;
quit;

data a.base1;
set a.base;
if secteur='F' then secF=1;
else secF=0;
run;

proc reg data=a.base1;
model lprod=leffectif lpub secF;
output out=resss p=fit r=res ;
run;