sas如何实现脉冲响应函数图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

本人用sas做了脉冲响应分析，想得到脉冲响应函数图，请大侠指导

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：脉冲响应函数 脉冲响应 如何实现 响应函数 如何

可以自己编写

谢谢，我自己已经搞定了。

   /* set the graphics environment */
goptions reset=global gunit=pct border cback=white
         colors=(black blue green red)
         ftitle=swissb ftext=swiss htitle=6 htext=3;
  /* create the data set STOCKS */
data stocks;
   input year @7  hdate date9. @17 high
              @26 ldate date9. @36 low;
   format hdate ldate date9.;
   datalines;
1955  30DEC1955  488.40  17JAN1955  388.20
1956  06APR1956  521.05  23JAN1956  462.35
1957  12JUL1957  520.77  22OCT1957  419.79
1958  31DEC1958  583.65  25FEB1958  436.89
1959  31DEC1959  679.36  09FEB1959  574.46
1960  05JAN1960  685.47  25OCT1960  568.05
1961  13DEC1961  734.91  03JAN1961  610.25
1962  03JAN1962  726.01  26JUN1962  535.76
1963  18DEC1963  767.21  02JAN1963  646.79
1964  18NOV1964  891.71  02JAN1964  768.08
1965  31DEC1965  969.26  28JUN1965  840.59
1966  09FEB1966  995.15  07OCT1966  744.32
1967  25SEP1967  943.08  03JAN1967  786.41
1968  03DEC1968  985.21  21MAR1968  825.13
1969  14MAY1969  968.85  17DEC1969  769.93
1970  29DEC1970  842.00  06MAY1970  631.16
1971  28APR1971  950.82  23NOV1971  797.97
1972  11DEC1972 1036.27  26JAN1972  889.15
1973  11JAN1973 1051.70  05DEC1973  788.31
1974  13MAR1974  891.66  06DEC1974  577.60
1975  15JUL1975  881.81  02JAN1975  632.04
1976  21SEP1976 1014.79  02JAN1976  858.71
1977  03JAN1977  999.75  02NOV1977  800.85
1978  08SEP1978  907.74  28FEB1978  742.12
1979  05OCT1979  897.61  07NOV1979  796.67
1980  20NOV1980 1000.17  21APR1980  759.13
1981  27APR1981 1024.05  25SEP1981  824.01
1982  27DEC1982 1070.55  12AUG1982  776.92
1983  29NOV1983 1287.20  03JAN1983 1027.04
1984  06JAN1984 1286.64  24JUL1984 1086.57
1985  16DEC1985 1553.10  04JAN1985 1184.96
1986  02DEC1986 1955.57  22JAN1986 1502.29
1987  25AUG1987 2722.42  19OCT1987 1738.74
1988  21OCT1988 2183.50  20JAN1988 1879.14
1989  09OCT1989 2791.41  03JAN1989 2144.64
1990  16JUL1990 2999.75  11OCT1990 2365.10
1991  31DEC1991 3168.83  09JAN1991 2470.30
1992  01JUN1992 3413.21  09OCT1992 3136.58
1993  29DEC1993 3794.33  20JAN1993 3241.95
1994  31JAN1994 3978.36  04APR1994 3593.35
1995  13DEC1995 5216.47  30JAN1995 3832.08
;
   /* restructure the data so that      */
   /* it defines a closed area          */
   /* create the data sets HIGH and LOW */
data high(keep=year value)
     low(keep=year value);
   set stocks;
   value=high; output high;
   value=low; output low;
run;
   /* reverse order of the observations in LOW */
proc sort data=low;
   by descending year;
run;
   /* concatenate HIGH and LOW to create data set AREA */
data area;
   set high low;
run;
   /* define titles and footnote */
title1 'Dow Jones Industrial Average';
title2  h=4 'Highs and Lows From 1955 to 1995';
footnote j=l ' Source: 1997 World Almanac'
         j=r 'GSYFAPL1 ';
  /* define symbol characteristics */
symbol interpol=m3n90  /* map/plot pattern fill  */
       cv=red                   /* pattern fill color        */
       co=blue;                  /* outline color            */
   /* define axis characteristics */
axis1 order=(1955 to 1995 by 5)
      label=none
      major=(height=2)
      minor=(number=4 height=1)
      offset=(2,2)
      width=3;
axis2 order=(0 to 5500 by 500)
      label=none
      major=(height=1.5)  offset=(0,0)
      minor=(number=1 height=1);
  /* generate plot using data set AREA */
proc gplot data=area;
   plot value*year / haxis=axis1
                     vaxis=axis2
                     vref=(1000 3000 5000);
run;
quit;

liping7939 发表于 2012-11-16 21:56
谢谢，我自己已经搞定了。

你好，遇到跟您一样的问题，你是怎么搞定的？谢谢！

同问楼主是如何解决的？
我在varmax程序中加入了plots=(impulse forecast)，但是得到的脉冲响应图貌似和书上不太一样。。。

ods listing gpath='F:\SAS';
ods graphics on;
proc varmax data = ndgdp plot=impulse; /*脉冲响应函数与预测方差分解*/
model ddy1t ddy2t ddy3t ddy4t/p=3 lagmax=6
print = (IMPULSE = (simple) decompose(6)) printform = both;
run;quit;
ods graphics off;

liping7939 发表于 2012-11-16 21:56
谢谢，我自己已经搞定了。

厉害！！！！