5 pièce(s) jointe(s)
Prévision de la consomamtion électrique du japon
Bonjour,
Je réalise un devoir sur SAS, cependant je suis ultra néophyte. De plus mon partenaire qui détient le logiciel chez lui, n'est pas la du weekend donc je me sens perdu dans mon analyse.
Je dois realiser la prévision de la consomamtion électrique du japon pour les années à venir. Mon camarade m'a pondu ceci, mais je ne suis vraiment pas sur de son modèle. Pouvez vous m'aider s'il vous plait ?
Code SAS
Code:
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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
libname U "D:\Sas";
/* importer des données excel*/
procimport datafile = 'D:\Sas\Japon.xlsx'
out = U.Japon
dbms = xlsx
replace;
run;
data U.Japon;
set U.Japon;
if electricity= "." thendelete;
run;
goptionreset=global;
procgplot data=U.Japon;
plot electricity*date;
symbol1interpol=join;
title"Electricity Supplied in Japon";
run;
quit;
procarima data=U.Japon;
identifyvar=electricity;
run;
quit;
procarima data=U.Japon;
identifyvar=electricity(1,14);
estimateq = (1)(14) p = 1noint method = cls;
run;
quit;
procarima data= U.Japon;
identifyvar = electricity(1,14);
estimatep=1q=(1)(14);
forecastid=date lead=14interval=month out=u.elec_Japon_sarima;
run;
procprint data =u.elec_Japon_sarima;
run;
data u.Japon_forecast;
merge u.Japon u.elec_Japon_sarima;
by date;
if date <mdy(6,1,2019) then forecast2=.;
if date >= mdy(6,1,2019) then forecast2=forecast;
run;
goptionreset=global;
procgplot data=u.Japon_forecast;
plot electricity*date forecast2*date / overlaylegend=legend1 href=11354;
symbol1color=black interpol=join;
symbol2color=black interpol=join line=2;
legend1label=none value=('electricity_prev');
title"Prevision";
run;
quit;[ATTACH=CONFIG]524450[/ATTACH] |
