全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Gauss专版
2012-2-10 21:18:29
zhangtao 发表于 2012-2-10 09:21
epoh老师,您好!
      有空帮我看看以下程序,谢谢!
zhangtao兄
请注意短信息
二维码

扫码加我 拉你入群

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

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

2012-2-14 23:54:40
tsp1.zip
大小:(7.94 KB)

 马上下载

本附件包括:

  • tsp.xls

epoh老师,您帮我试一下,我没有发现错误,可是觉得运行结果有点奇怪。
二维码

扫码加我 拉你入群

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

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

2012-2-15 12:06:08
epoh老师,我想问下跑Gauss程序时提示错误:G0562 Open /* */ comment
应该如何解决?
二维码

扫码加我 拉你入群

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

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

2012-2-15 13:18:57
ywh19860616 发表于 2012-2-15 12:06
epoh老师,我想问下跑Gauss程序时提示错误:G0562 Open /* */ comment
应该如何解决?
程序检查过没错.
ntrial=5000,我做出来的结果跟你接近,
可能跟数据有关.
二维码

扫码加我 拉你入群

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

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

2012-2-15 15:49:50
ywh19860616 发表于 2012-2-15 12:06
epoh老师,我想问下跑Gauss程序时提示错误:G0562 Open /* */ comment
应该如何解决?
请参考语法错误对照表
若排除有困难,请贴上程序
http://www.aptech.com/manuals/UserGuide11.pdf
chap 24 Error Messages
   G0002 File too large
   .....
   G0562 Open /* */ comment
   .....
   G0603 Failure to converge
二维码

扫码加我 拉你入群

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

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

2012-2-16 11:59:00
epoh 发表于 2012-2-15 15:49
请参考语法错误对照表
若排除有困难,请贴上程序
http://www.aptech.com/manuals/UserGuide11.pdf
好的,谢谢epoh老师,我先尝试看看
二维码

扫码加我 拉你入群

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

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

2012-2-21 10:37:08
ywh19860616 发表于 2012-2-16 11:59
好的,谢谢epoh老师,我先尝试看看
刚检查一下,
数据的dimension不同,需要修改
修改部分及结果如下:

fmlad.ox  请修改main()如下
原: // indexes: APX = 4, EEX = 5, EXAA = 6, Powernext = 9
改: // indexes: APX = 1, EEX = 2, EXAA = 3, Powernext = 6

main()
{
        decl asNames;
        decl mData = loadmat("MeanPrices.xls",&asNames)[168:][];
        // indexes: APX = 4, EEX = 5, EXAA = 6, Powernext = 9
        //decl iYindex = 9;
        //decl vXindex = 5;
//revised
        // indexes: APX = 1, EEX = 2, EXAA = 3, Powernext = 6
        decl iYindex = 6;
        decl vXindex = 2;
        decl vY = mData[][iYindex];
        println("\nY variable : ", asNames[iYindex]);
        decl mX = mData[][vXindex];
        println("X variables: ", asNames[vXindex]);
        fmlad(vY, mX, -1);
}

////////results:
Y variable : Powermean

LAD regression
       Estim.         Std.Err.      t-ratio         p-value
     -0.47472       1.7487     -0.27147      0.78603
      0.98907     0.053248    18.575        5.1496e-077

FM-LAD regression (Bartlett kernel, Bandwidth = 4)
       Estim.         Std.Err.      t-ratio          p-value
     -0.54928       1.7487      -0.31411      0.75343
      0.99079     0.014468       68.483      0.00000

Pseudo R2       = 0.721374
Regression MAD  = 3.4057
Regression RMSE = 6.3673
LAD-AR1 coeff   = 0.269091
OLS-AR1 coeff   = 0.347363
二维码

扫码加我 拉你入群

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

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

2012-2-21 11:29:16
epoh 发表于 2012-2-21 10:37
刚检查一下,
数据的dimension不同,需要修改
修改部分及结果如下:
非常感谢epoh老师
nh.ox问题我大致找到了,是作者把里面的程序用注释形式引起来了一部分。
二维码

扫码加我 拉你入群

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

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

2012-5-20 10:23:04
epoh 发表于 2012-2-7 14:02
bootstrap一开始产生的ystarm,不具代表性,必须去除.
这有点类似MCMC的burn in.
所以nall=69,就要去除前 ...
epoh老师,以下是我的data file,我增加了一项zs.(-2),
可是提示有错误,请您帮忙看下。

谢谢。
tspeast.rar
大小:(4.52 KB)

 马上下载

本附件包括:

  • tspeast.csv


二维码

扫码加我 拉你入群

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

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

2012-5-20 19:36:16
ywh19860616 发表于 2012-5-20 10:23
epoh老师,以下是我的data file,我增加了一项zs.(-2),
可是提示有错误,请您帮忙看下。
修正了两个地方:
SET nc=mlagy+1+1+1;  ?constant,z(-1) z(-2)

SET ystarm(u,v)=beta(v,1)+beta(v,2)*ystarm(u1,v)+beta(v,3)*ystarm(u2,v)+beta(v,4)*zstarm(u1,v)+beta(v,5)*zstarm(u2,v)+epsstarm(u,v);
二维码

扫码加我 拉你入群

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

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

2012-5-20 20:10:00
epoh 发表于 2012-5-20 19:36
修正了两个地方:
SET nc=mlagy+1+1+1;  ?constant,z(-1) z(-2)
epoh老师,我对zstarm(u1,v)以及zstarm(u2,v)还有一个疑问。
您之前提到过:

ystarm(u,v)=beta(v,1)+beta(v,2)*ystarm(u1,v)+beta(v,3)*ystarm(u2,v)+beta(v,4)*tzstarm(u1,v)+epsstarm(u,v);
请更改为
ystarm(u,v)=beta(v,1)+beta(v,2)*ystarm(u1,v)+beta(v,3)*ystarm(u2,v)+beta(v,4)*tzstarm(u,v)+epsstarm(u,v);

上面是只包含z(-1)时的程序,现在是加入z(-2),是这样修改吗?
二维码

扫码加我 拉你入群

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

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

2012-5-20 20:46:49
ywh19860616 发表于 2012-5-20 20:10
epoh老师,我对zstarm(u1,v)以及zstarm(u2,v)还有一个疑问。
您之前提到过:
u1=u-1;(lag 1)
u2=u-2;(lag 2)

ystarm,zstarm两者公式相同,只是系数不同

beta(v,2)*ystarm(u1,v)+beta(v,3)*ystarm(u2,v)+
beta(v,4)*zstarm(u1,v)+beta(v,5)*zstarm(u2,v)
二维码

扫码加我 拉你入群

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

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

2012-5-20 21:18:51
epoh 发表于 2012-5-20 20:46
u1=u-1;(lag 1)
u2=u-2;(lag 2)
epoh老师,从公式来看,的确是这样,
我先看看,谢谢您。
二维码

扫码加我 拉你入群

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

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

2012-5-22 08:55:20
epoh 发表于 2012-5-20 20:46
u1=u-1;(lag 1)
u2=u-2;(lag 2)
epoh老师,您好,
我感觉这程序还存在一些问题,比如:
我对x、y、z的滞后分别取为1,1,2时,按
原来程序应该要把:
set mlagy=2 改为set mlagy=1
setmlagx=2 改为 set mlagx=1
然后下面的方程也对应修改。
可是,我运行的时候既然不修改mlagy和mlagx滞后为1,也
能出结果。这个不理解,难道mlagy和mlagx这两个变量在新
修改的程序中没有作用吗?

请您有空帮我核对下程序,谢谢您。
二维码

扫码加我 拉你入群

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

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

2012-5-22 11:03:30
ywh19860616 发表于 2012-5-22 08:55
epoh老师,您好,
我感觉这程序还存在一些问题,比如:
我对x、y、z的滞后分别取为1,1,2时,按
呵呵,可以运行并不代表对
恰巧你是1,1,2
取样范围没受影响
如果是1,1,1
取样范围就会影响结果
当然不只这些,公式也要配合更改

SET mlagy=2;
SET mlagx=2;
SET maxlag=mlagy+POS(mlagx-mlagy);  ?2
SET cfyear=fyear+maxlag;                      ?1992
SET cnyear=lyear-cfyear;                        ?18
SET syear=lyear-nall+1;                          ?1969
SET firsty=maxlag+2;                             ?4
SET startu=nextra+maxlag+1;               ?33
SET byear=fyear-startu;                         ?1957

Current sample:  1992 to 2010

#################
SET mlagy=1;
SET mlagx=1;
SET maxlag=mlagy+POS(mlagx-mlagy);  ?1
SET cfyear=fyear+maxlag;                      ?1991
SET cnyear=lyear-cfyear;                        ?19
SET syear=lyear-nall+1;                          ?1969
SET firsty=maxlag+2;                             ?3
SET startu=nextra+maxlag+1;               ?32
SET byear=fyear-startu;                          ?1958

Current sample:  1991 to 2010
二维码

扫码加我 拉你入群

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

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

2012-5-22 11:34:30
epoh 发表于 2012-5-22 11:03
呵呵,可以运行并不代表对
恰巧你是1,1,2
取样范围没受影响
epoh老师,谢谢您。
您的意思是不是说mlagy和mlagx也和公式一样
对应着修改,程序还是没有问题的?
二维码

扫码加我 拉你入群

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

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

2012-5-22 11:51:09
ywh19860616 发表于 2012-5-22 11:34
epoh老师,谢谢您。
您的意思是不是说mlagy和mlagx也和公式一样
对应着修改,程序还是没有问题的?
这个程序我对照过取样
没错的
二维码

扫码加我 拉你入群

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

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

2012-5-22 11:56:52
epoh 发表于 2012-5-22 11:51
这个程序我对照过取样
没错的
谢谢epoh老师,呵呵,这样我就没有疑问了
早上看的时候,感觉mlagy和mlagx没有相应
修改也能得到答案,所以困惑了。
现在我也对mlagx和mlagy做相应修改就可以了。
二维码

扫码加我 拉你入群

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

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

2012-5-22 12:22:20
epoh 发表于 2012-5-22 11:51
这个程序我对照过取样
没错的
epoh老师,在您那修改mlagy=1和mlagx=1后,z滞后取2.其他公式也对应改变,程序可以正确运行?
二维码

扫码加我 拉你入群

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

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

2012-5-22 15:27:16
ywh19860616 发表于 2012-5-22 12:22
epoh老师,在您那修改mlagy=1和mlagx=1后,z滞后取2.其他公式也对应改变,程序可以正确运行?
SET mlagy=1;
SET mlagx=1;
SET mlagz=2;
请注意短信息
二维码

扫码加我 拉你入群

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

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

2012-5-22 16:16:49
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中最大值?
二维码

扫码加我 拉你入群

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

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

2012-5-22 18:30:07
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);
二维码

扫码加我 拉你入群

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

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

2012-5-22 20:03:49
epoh 发表于 2012-5-22 18:30
maxlag定阶方式是取mlagx,mlagy, mlagz中最大值
对,要不然底下是子会产生错误
是的,谢谢epoh老师
您这样修改后,运行其他的组合就方便许多了。
谢谢您。
二维码

扫码加我 拉你入群

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

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

2012-5-23 14:50:03
epoh 发表于 2012-5-22 18:30
maxlag定阶方式是取mlagx,mlagy, mlagz中最大值
对,要不然底下是子会产生错误
epoh老师,在程序中设置SET nall=42;即nall=2*T。
现在我的T为21,我想请教下抽样时是舍弃前面的21次?
就是“burn-in”。这个在程序哪处体现?

  
二维码

扫码加我 拉你入群

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

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

2012-5-23 15:50:45
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
二维码

扫码加我 拉你入群

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

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

2012-5-23 16:24:00
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.
我想问下,这个警告要紧吗?
二维码

扫码加我 拉你入群

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

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

2012-5-23 19:08:58
ywh19860616 发表于 2012-5-23 16:24
恩,清楚了,谢谢您。
如果不要把
?print zstarm;
请注意短信息
tspeast_revised.xls
大小:(28.87 KB)

 马上下载

二维码

扫码加我 拉你入群

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

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

2012-5-23 19:45:08
epoh 发表于 2012-5-23 19:08
请注意短信息
非常感谢epoh老师
谢谢您,我先试一下,有不理解之处再请教您。
二维码

扫码加我 拉你入群

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

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

2012-5-27 15:01:16
epoh 发表于 2012-5-23 19:08
请注意短信息
epoh老师,我想请教一下,我最近跑了一个程序,矩阵有点大,程序中出现大概3000万*100和1000万*100的矩阵,然后gauss就显示 "the results are too large". 请问我该如何设置gauss能让我跑这个程序呢?谢谢
二维码

扫码加我 拉你入群

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

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

2012-5-28 20:00:17
theawang 发表于 2012-5-27 15:01
epoh老师,我想请教一下,我最近跑了一个程序,矩阵有点大,程序中出现大概3000万*100和1000万*100的矩阵, ...
哈哈,通常比较大的矩阵
我都用matlab
二维码

扫码加我 拉你入群

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

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

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群