<p>data zlh;<br/>  input air@@;<br/>  date=intnx('month','1jan49'd,_n_-1);<br/>  format date monyy.;<br/>  cards;<br/>  112 118 132 129 121 135 148 148 136 119 104 118 <br/>  115 126 141 135 125 149 170 170 158 133 114 140 <br/>  145 150 178 163 172 178 199 199 184 162 146 166<br/>  171 180 193 181 183 218 230 242 209 191 172 194<br/>  196 196 236 235 229 243 264 302 293 259 229 203<br/>  229 242 233 267 269 270 315 364 347 312 274 237 <br/>  278 284 277 317 313 318 374 413 405 355 306 271 <br/>  306 315 301 348 355 422 465 467 404 347 305 336 <br/>  315 301 356 348 355 422 465 467 404 347 305 336 <br/>  340 318 362 348 363 435 491 505 404 359 310 337 <br/>  360 342 406 396 420 472 548 559 463 407 362 405 <br/>  417 391 419 461 472 535 622 606 508 461 390 432<br/>  ;<br/>run;<br/>proc gplot data=Zlh;<br/> symbol1 i=join v=dot c=red;<br/> plot air*date=1;<br/>run;<br/>data lair;<br/> set Zlh;<br/> lair=log(air);<br/>run;<br/>proc gplot data=lair;<br/> symbol2 i=spline c=green;<br/> plot lair*date=2;<br/>run;<br/>proc arima data=lair;<br/> identify var=lair nlag=36;<br/>run;<br/>proc arima data=lair;<br/> identify var=lair(1) nlag=36;<br/>run;<br/>proc arima data=lair;<br/> identify var=lair(1,12) nlag=36;<br/>run;<br/>proc arima data=lair;<br/> identify var=lair(1,12) nlag=36;<br/> estimate q=(1)(12) noconstant method=uls plot;<br/>run;<br/>proc arima data=lair;<br/> identify var=lair(1,12) nlag=36;<br/> estimate q=(1)(12) noconstant method=uls plot;<br/> forecast lead=36 interval=month id=date out=b;<br/>run;<br/>proc print data=b;<br/>run;<br/>data lair;<br/>set b;<br/>air=exp(lair);<br/>forecast=exp(forecast+std*std/2);<br/>l95=exp(l95);<br/>u95=exp(u95);<br/>run;<br/>proc print data=lair;<br/>run;<br/>proc gplot data=lair<font color="#f73809" size="4">;{<strong>就是从这里开始做不下去了，做出来就是最下面这样的提示</strong>}<br/></font>symbol1 I=none v=star r=1 c=red;<br/>symbol2 I=join v=plus r=1 c=green;<br/>symbol3 I=join v=none l=3 r=1 c=blue;<br/>proc gplot data=lair;<br/>where data>=’1jan59’d;<br/>plot air*date=1 forecast*date=2 l95*date=3 u95*date=3/<br/>overlay haxis=’1jan59’d to ‘1jan62’d by year;<br/>run;</p><p><font size="3">{<font color="#4822dd">551  proc gplot data=lair;<br/>552  where data>=’1jan59’d;<br/>ERROR: Variable data is not on file WORK.LAIR.<br/>NOTE: SCL source line.<br/>553  plot air*date=1 forecast*date=2 l95*date=3 u95*date=3/overlay haxis=’1jan59’d to ‘1jan62’d by year;</font></font></p><p><font color="#4822dd" size="3">                                                                                     --<br/>                                                                                     22<br/>                                                                                        -----------<br/>                                                                                        202<br/>ERROR: "’1JAN59’D" is not a valid name.<br/>NOTE: The previous statement has been deleted.<br/>ERROR 22-322: Syntax error, expecting one of the following: ;, ANNOTATE, AREAS, AUTOHREF, AUTOVREF, C, CAXIS, CFRAME, CHREF,<br/>              COUTLINE, CTEXT, CVREF, DATAORDER, DESCRIPTION, FRAME, HAXIS, HMINOR, HREF, HTML, HTML_LEGEND, HZERO, IFRAME,<br/>              IMAGESTYLE, LEGEND, LHREF, LVREF, NAME, NOAXIS, NOFRAME, NOLEGEND, OVERLAY, REGEQN, SKIPMISS, VAXIS, VMINOR, VREF,<br/>              VREVERSE, VZERO.<br/>ERROR 202-322: The option or parameter is not recognized and will be ignored.<br/>554  run;</font></p><p><font color="#4822dd" size="3">ERROR: At least one PLOT or BUBBLE statement must be given.}</font></p><p><font color="#4822dd" size="3">请大家告诉是什么问题，我是按书上来的，可是就是做不出来最后的预测图，</font></p>