程序运行问题

zhangtao 发表于 2012-2-10 09:21
epoh老师，您好！
      有空帮我看看以下程序，谢谢！

zhangtao兄
请注意短信息

tsp1.zip (7.94 KB) 本附件包括：

• tsp.xls

2012-2-14 23:54:33 上传

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

epoh老师，我想问下跑Gauss程序时提示错误：G0562 Open /* */ comment
应该如何解决？

ywh19860616 发表于 2012-2-15 12:06
epoh老师，我想问下跑Gauss程序时提示错误：G0562 Open /* */ comment
应该如何解决？

程序检查过没错.
ntrial=5000,我做出来的结果跟你接近,
可能跟数据有关.

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

epoh 发表于 2012-2-15 15:49
请参考语法错误对照表
若排除有困难,请贴上程序
http://www.aptech.com/manuals/UserGuide11.pdf

好的，谢谢epoh老师，我先尝试看看

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

epoh 发表于 2012-2-21 10:37
刚检查一下,
数据的dimension不同,需要修改
修改部分及结果如下:

非常感谢epoh老师
nh.ox问题我大致找到了，是作者把里面的程序用注释形式引起来了一部分。

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 10:22:59 上传

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);