Discussion :

[med79]

Bonjour,

Je souhaiterais désaisonnaliser la série trimestrielle ci-jointe et ainsi calculer les coefficients saisonniers.

J'ai utilisé la procédure X12 mais je n'ai pas tous compris des résultats ( une boite noire).

Je vous remercie d'avance de votre aide.

Voici mon programme :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
Data WORK.Donnees    ;
     %Let _EFIERR_ = 0; 
     Infile "C:\Users\Desktop\test2.csv"
            Delimiter = ';' 
            MISSOVER 
            Lrecl=32767 
            Firstobs=2 ;
	 Informat annee best32. ;
	 Informat trimestre best32. ;
     Informat Date yyQ7. ;
	 Informat Var numx5.1 ;
     Format Date yyQ7. ;
	 Input Annee Trimestre Date VAR;
     If _ERROR_ then call symputx('_EFIERR_',1);  
Run;
 
Proc means data = donnees;
Class annee;
Var Var;
Run;
 
Proc means data = donnees;
Class trimestre;
Var Var;
Run;
 
Proc x12 data=donnees date=date;
      var Var;
      x11;
      ods select d1 d10d d11 d12 d13;;
Run;
 
Proc timeseries data=donnees plot=series;
   id date interval=quarter;
   var Var;
Run;
/*
Proc timeseries data=donnees outseason=season
                plots=(series TC SC IC SA);
   id date interval=quarter;
   var Var;
   decomp ORIG TC SC IC SA/mode=MULTORADD;
Run;*/
 
title1 "Seasonal Statistics";
Proc sgplot data=season;
   series x=_season_ y=max  / lineattrs=(pattern=solid);
   series x=_season_ y=mean / lineattrs=(pattern=solid);
   series x=_season_ y=min  / lineattrs=(pattern=solid);
   yaxis display=(nolabel);
   xaxis values=(1 to 4 by 1);
Run;

Bonne journée