楼主: zodica2008
3661 9

[原创博文] 各位大侠我有个ARIMA的用SAS做到有一部做不下去了 [推广有奖]

  • 0关注
  • 0粉丝

高中生

80%

还不是VIP/贵宾

-

0%
威望
0 级
论坛币
2010 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
151 点
帖子
26
精华
0
在线时间
9 小时
注册时间
2008-6-23
最后登录
2021-3-24
楼主
zodica2008 发表于 2009-2-26 10:48:00 |只看作者 |坛友微信交流群|倒序 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币
<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>
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0
关键词：ARIMA 各位大侠 Rim ima Description

相关帖子

回帖推荐

buqinghuan 发表于2楼  查看完整内容

data zlh;  input air@@;  date=intnx('month','1jan49'd,_n_-1);  format date monyy.;  cards;  112 118 132 129 121 135 148 148 136 119 104 118   115 126 141 135 125 149 170 170 158 133 114 140   145 150 178 163 172 178 199 199 184 162 146 166  171 180 193 181 183 218 230 242 209 191 172 194  196 196 236 235 229 243 264 302 293 259 229 203  229 242 233 ...

本帖被以下文库推荐

使用道具

沙发
buqinghuan 发表于 2009-2-26 11:13:00 |只看作者 |坛友微信交流群
data zlh;
  input air@@;
  date=intnx('month','1jan49'd,_n_-1);
  format date monyy.;
  cards;
  112 118 132 129 121 135 148 148 136 119 104 118
  115 126 141 135 125 149 170 170 158 133 114 140
  145 150 178 163 172 178 199 199 184 162 146 166
  171 180 193 181 183 218 230 242 209 191 172 194
  196 196 236 235 229 243 264 302 293 259 229 203
  229 242 233 267 269 270 315 364 347 312 274 237
  278 284 277 317 313 318 374 413 405 355 306 271
  306 315 301 348 355 422 465 467 404 347 305 336
  315 301 356 348 355 422 465 467 404 347 305 336
  340 318 362 348 363 435 491 505 404 359 310 337
  360 342 406 396 420 472 548 559 463 407 362 405
  417 391 419 461 472 535 622 606 508 461 390 432
  ;
run;
proc gplot data=Zlh;
 symbol1 i=join v=dot c=red;
 plot air*date=1;
run;
data lair;
 set Zlh;
 lair=log(air);
run;
proc gplot data=lair;
 symbol2 i=spline c=green;
 plot lair*date=2;
run;
proc arima data=lair;
 identify var=lair nlag=36;
run;
proc arima data=lair;
 identify var=lair(1) nlag=36;
run;
proc arima data=lair;
 identify var=lair(1,12) nlag=36;
run;
proc arima data=lair;
 identify var=lair(1,12) nlag=36;
 estimate q=(1)(12) noconstant method=uls plot;
run;
proc arima data=lair;
 identify var=lair(1,12) nlag=36;
 estimate q=(1)(12) noconstant method=uls plot;
 forecast lead=36 interval=month id=date out=b;
run;
proc print data=b;
run;
data lair;
set b;
air=exp(lair);
forecast=exp(forecast+std*std/2);
l95=exp(l95);
u95=exp(u95);
run;
proc print data=lair;
run;
proc gplot data=lair;
symbol1 I=none v=star r=1 c=red;
symbol2 I=join v=plus r=1 c=green;
symbol3 I=join v=none l=3 r=1 c=blue;
proc gplot data=lair;
where date>='1jan59'd;
plot air*date=1 forecast*date=2 l95*date=3 u95*date=3/
overlay haxis='1jan59'd to '1jan62'd by year;
run;
已有 1 人评分经验 论坛币 收起 理由
bakoll + 3 + 3 精彩帖子

总评分: 经验 + 3  论坛币 + 3   查看全部评分

使用道具

藤椅
buqinghuan 发表于 2009-2-26 11:14:00 |只看作者 |坛友微信交流群
这回试试

使用道具

板凳
456852 发表于 2009-2-26 12:08:00 |只看作者 |坛友微信交流群
proc gplot data=lair;{就是从这里开始做不下去了，做出来就是最下面这样的提示}
symbol1 I=none v=star r=1 c=red;
symbol2 I=join v=plus r=1 c=green;
symbol3 I=join v=none l=3 r=1 c=blue;
proc gplot data=lair;

where data>=’1jan59’d;（童鞋。。变量名叫date，不是data!）

plot air*date=1 forecast*date=2 l95*date=3 u95*date=3/
overlay haxis=’1jan59’d to ‘1jan62’d by year;
run;

使用道具

报纸
rdzr 发表于 2009-2-26 13:07:00 |只看作者 |坛友微信交流群
lz,偶个人认为，用 sas 之 ARIMA 过程 作时间序列预测并不是一个好的选择，推荐用 SPSS ，这方面 SPSS 比 SAS 直观，操作简单，而预测结果与残差均可保存为单独的列变量，而 SAS 就没有这方面的优势。所以，建议 LZ 不妨试一试SPSS 之 ARIMA。

使用道具

地板
456852 发表于 2009-2-26 13:25:00 |只看作者 |坛友微信交流群
楼上，谁和你说SAS不能保存预测结果和残差？
我就经常做预测值与真实值的比较图和残差图。

使用道具

7楼
zodica2008 发表于 2009-2-26 17:08:00 |只看作者 |坛友微信交流群
谢谢各位，问题已经解决，就是那个date弄错了

[此贴子已经被作者于2009-2-26 17:13:34编辑过]

使用道具

8楼
rdzr 发表于 2009-2-27 15:11:00 |只看作者 |坛友微信交流群
6楼的，我说SAS不能保存 预测结果和残差 了吗？ 我也经常用 sas，我的意思是说 sas 与 spss 各有所长，明白不？！

使用道具

9楼
456852 发表于 2009-2-28 01:17:00 |只看作者 |坛友微信交流群
而预测结果与残差均可保存为单独的列变量，而 SAS 就没有这方面的优势。
===============================================
这是你的原话不?
各个软件都各有所长这不假，但是也得实事求是的分析。

使用道具

10楼
wyxdlzw 发表于 2011-5-1 20:07:18 |只看作者 |坛友微信交流群
您好，能不能麻烦您把每一步的注释标一下呢？非常感谢，我是刚入门sas。请求帮助！

使用道具

返回列表
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
  • 手机版  |
  • 意见反馈  |
  • 帮助  |
  • 新手入门  |
  • 用户手册  |
  • 友情链接  |

    • 如有投资本站、合作意向或投放广告，请联系：13661292478（刘老师）

    联系客服

    邮箱：service@pinggu.org 投诉或不良信息处理：（010-68466864）

    010-68466864

    京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问：王进律师 知识产权保护声明   免责及隐私声明

    GMT+8, 2024-9-17 14:07