/*creer un fichier ipi avec la date et l'imprimer*/ data a; set chron9.ipi; t=intnx('qtr','1jan63'd,_n_-1); run; proc print data=a (obs=5); format t yyq.; run; /*Dessiner la chronique*/ title1 'Indice Production Industrielle trimestrielle'; proc sgplot data=a; series y=ipi x=t; format t yyq.; run; /* on verifie que le filtrage arithmetique a 4 termes supprime la saisonnalite */ proc expand data=a out=c; /* convert ipi = moymob / transform=(cmovave (0.5 1 1 1 0.5)); */ convert ipi = moymob / transform=(cmovave 4); run; /* calcul de la (pseudo) composante saisonnière */ data c; set c; saison=ipi-moymob; run; title "La chronique IPI et sa moyenne mobile centrée"; proc sgplot data=c; series y=ipi x=t; series y=moymob x=t; format t yyq.; run; /* on voit bien l'influence des chocs pétroliers 1973 et 1979 avec effet retard */ title "Composante pseudo-périodique"; proc sgplot data=c; series y=saison x=t; format t yyq.; run; title "Illustration de la décomposition de ipi"; proc timeseries data=a plots=(series decomp TC SC IC SA); id t interval=qtr; var ipi; run; /* idem: influence des chocs pétroliers + mai 68 que l'on voit dans les données CVS et dans les résidus */