*/analyse en composantes principales/*; data tempville; infile "C:/Users/christine/desktop/SAS/temp.dat" ; input ville $ janv fevr mars avri mai juin juil aout sept octo nove dece; run; proc print; run; proc princomp data=tempville cov plots(ncomp=2)=(pattern(circles=1.0)score); var janv--dece; id ville; run; data ocde; infile "C:/Users/christine/desktop/SAS/ocde.dat"; input pays $ natal chomage a_prim a_sec pib fbcf infl recc m_inf prot nrj ; annee = 1975+2*mod(_n_ - 1, 4); run; proc univariate data=ocde normal plot; var natal--nrj; run; proc sort data=ocde out=tocde; by year; proc univariate data=tocde plot; var natal--nrj; by annee; run; proc corr data=ocde plots=matrix; var natal--nrj; run; proc princomp data=ocde plots(ncomp=2)=(pattern(circles=1.0)score); var natal--nrj; id pays; run; %acp(ocde,pays,natal--nrj); %gacpbx;%gacpsx; %gacpvx; %gacpix; data p_ocde; set ocde; if pays eq "PO" then poids=0 else ; Poids=1; run; %acp(p_ocde,pays,natal--nrj,poids=poids); %gacpbx;%gacpsx; %gacpvx; %gacpix; %let x=1; %let y=2; data anno; retain xsys ysys ’2’; set coorindq nobs=nind; style=’swiss’; if mod(_n_, 4) ne 1 then delete ; y= prin&y; x= prin&x; text=ident; size=1; run; proc gplot data= coorindq; title; footnote ; axis1 length=14cm; / * attention taille * / axis2 length=8cm; symbol v=dot i=join r=13 height=.5; plot prin&y * prin&x=ident / annotate=anno frame href=0 vref=0 nolegend haxis=axis1 vaxis=axis2; run; goptions reset=all; quit; data ukcomp1 ; infile 'C:/Users/Christine/Desktop/SAS/ukcomp1.dat' dlm='09'x; input RETCAP GEARRAT CAPINT WCFTDT LOGSALE LOGASST CURRAT QUIKRAT NFATAST INVTAST FATTOT PAYOUT WCFTCL; poids=1; run; proc reg data=ukcomp1; model RETCAP=WCFTCL WCFTDT GEARRAT LOGSALE LOGASST NFATAST CAPINT FATTOT INVTAST PAYOUT QUIKRAT CURRAT; run; quit; data ukcomp2 ; infile 'C:/Users/Christine/Desktop/SAS/ukcomp1.dat' dlm='09'x; input RETCAP2 WCFTCL WCFTDT GEARRAT LOGSALE LOGASST NFATAST CAPINT FATTOT INVTAST PAYOUT QUIKRAT CURRAT; poids=0; run; data ukcomp; set ukcomp1 ukcomp2; run; proc princomp data=ukcomp out=comp; var WCFTCL WCFTDT GEARRAT LOGSALE LOGASST NFATAST CAPINT FATTOT INVTAST PAYOUT QUIKRAT CURRAT; weight poids; proc reg data=comp; model retcap=prin1--prin12/selection=rsquare cp best=1; run;quit; proc pls data=ukcomp1 cv=one nfac=10; model retcap=WCFTCL WCFTDT GEARRAT LOGSALE LOGASST NFATAST CAPINT FATTOT INVTAST PAYOUT QUIKRAT CURRAT; run; proc pls data=ukcomp nfac=2; model retcap=WCFTCL WCFTDT GEARRAT LOGSALE LOGASST NFATAST CAPINT FATTOT INVTAST PAYOUT QUIKRAT CURRAT; output out=plstat predicted=ychap; run; */ANOVA/*; data carotte; infile "C:/Users/christine/desktop/SAS/carotte.dat"; input Sol Variete Jour; run; proc glm data=carotte; class Variete Sol; model Jour=Variete Sol variete*Sol; run; quit;