程序运行问题

epoh 发表于 2012-5-22 15:27
SET mlagy=1;
SET mlagx=1;
SET mlagz=2;

epoh老师，您好
?SET maxlag=mlagy+POS(mlagx-mlagy);
SET maxlag=2;
我还想请教您，maxlag定阶方式是取mlagx，mlagy
mlagz中最大值？

ywh19860616 发表于 2012-5-22 16:16
epoh老师，您好
?SET maxlag=mlagy+POS(mlagx-mlagy);
SET maxlag=2;

maxlag定阶方式是取mlagx，mlagy, mlagz中最大值
对,要不然底下是子会产生错误

FRML eq. ys.=a.+b1.*ys.(-1)+c1.*xs.(-1)+d1.*zs.(-1)+d2.*zs.(-2);

epoh 发表于 2012-5-22 18:30
maxlag定阶方式是取mlagx，mlagy, mlagz中最大值
对,要不然底下是子会产生错误

是的，谢谢epoh老师
您这样修改后，运行其他的组合就方便许多了。
谢谢您。

epoh 发表于 2012-5-22 18:30
maxlag定阶方式是取mlagx，mlagy, mlagz中最大值
对,要不然底下是子会产生错误

epoh老师，在程序中设置SET nall=42;即nall=2*T。
现在我的T为21，我想请教下抽样时是舍弃前面的21次？
就是“burn-in”。这个在程序哪处体现？

  

ywh19860616 发表于 2012-5-23 14:50
epoh老师，在程序中设置SET nall=42;即nall=2*T。
现在我的T为21，我想请教下抽样时是舍弃前面的21次？
...

RESIDNULL beta eps;
OPTIONS NORESID;
DO bs=1 TO ntrial;
  BSTR eps epsstarm ystarm ;         
  SMPL syear lyear;          ?1969 2010
  UNMAKE ystarm ys1-ys12;
  print ys1;
  SMPL fyear lyear;          ?1990 2010
  print ys1;
  SELECT t>maxlag;
  SUR(SILENT) eq;
  
  DOT(VALUE=i) 1-12;
    IF choice(i)=1; THEN;
      DO;
        ANALYZ(SILENT) constr1. ;
        ANALYZ(SILENT) constr2. ;
        SET waldstatm(bs,i)=@WALD;
      ENDDO;
  ENDDOT;
ENDDO;

????????
Current sample:  1969 to 2010

                   YS1
1969           0.00000
1970           0.00000
1971           0.00000
1972           0.35828
1973           0.66968
1974           0.96253
1975           1.25189
1976           1.52328
1977           1.80099
1978           2.05117
1979           2.33767
1980           2.58059
1981           2.84942
1982           3.10118
1983           3.36972
1984           3.61968
1985           3.84803
1986           4.07669
1987           4.32239
1988           4.51712
1989           4.71783
1990           4.92725
1991           5.16040
1992           5.27736
1993           5.43503
1994           5.56883
1995           5.69835
1996           5.83428
1997           5.96448
1998           6.07569
1999           6.21769
2000           6.35735
2001           6.47893
2002           6.61587
2003           6.73905
2004           6.86939
2005           7.01200
2006           7.12818
2007           7.23465
2008           7.34393
2009           7.43982
2010           7.55028


Current sample:  1990 to 2010

                   YS1
1990           4.92725
1991           5.16040
1992           5.27736
1993           5.43503
1994           5.56883
1995           5.69835
1996           5.83428
1997           5.96448
1998           6.07569
1999           6.21769
2000           6.35735
2001           6.47893
2002           6.61587
2003           6.73905
2004           6.86939
2005           7.01200
2006           7.12818
2007           7.23465
2008           7.34393
2009           7.43982
2010           7.55028

epoh 发表于 2012-5-23 15:50
RESIDNULL beta eps;
OPTIONS NORESID;
DO bs=1 TO ntrial;

恩，清楚了，谢谢您。
如果不要把
?print zstarm;
?print tzstarm;
这两句命令如果没有注释了，然后运行，
会提示：
*** WARNING: frequency of series T        is    0.
我想问下，这个警告要紧吗？

ywh19860616 发表于 2012-5-23 16:24
恩，清楚了，谢谢您。
如果不要把
?print zstarm;

请注意短信息
tspeast_revised.xls (28.87 KB)

2012-5-23 19:09:14 上传

epoh 发表于 2012-5-23 19:08
请注意短信息

非常感谢epoh老师
谢谢您，我先试一下，有不理解之处再请教您。

epoh 发表于 2012-5-23 19:08
请注意短信息

epoh老师,我想请教一下，我最近跑了一个程序，矩阵有点大，程序中出现大概3000万*100和1000万*100的矩阵，然后gauss就显示 "the results are too large". 请问我该如何设置gauss能让我跑这个程序呢？谢谢

theawang 发表于 2012-5-27 15:01
epoh老师,我想请教一下，我最近跑了一个程序，矩阵有点大，程序中出现大概3000万*100和1000万*100的矩阵， ...

哈哈,通常比较大的矩阵
我都用matlab