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] |