请教 error G0058 : Index out of range

我运行的程序总是报错，拜请高手指点。附件里是程序和数据，可以直接运行的，但是运行时有时候会出现error G0058 : Index out of range ，但也有时候能正确运行出来，不知道问题在哪里？ 麻烦高手指点，谢谢
程序疑问.zip (6.28 KB)

2015-10-15 20:19:01 上传

最好先说出程序的功能，这个程序运行的最终结果应该没出来，感觉调用sorthat(yh,crit,kq)出了问题，crit不知道起什么作用？还有148行赋值好像有问题。

是你改写的程序还是源程序？

tianyuwuhe 发表于 2015-10-16 00:48
是你改写的程序还是源程序？

这个程序是用来进行俱乐部收敛检验的，我只是修改了原始程序开头的n和t（原始程序n=19,t=84），运行就总报错

tianyuwuhe 发表于 2015-10-16 00:47
最好先说出程序的功能，这个程序运行的最终结果应该没出来，感觉调用sorthat(yh,crit,kq)出了问题，crit不知 ...

我做出来结果指出下面这段有问题，请指教！！
dasi1:
  ww = yhh[.,ik1:ik1+1];
  {b1,t1} = logts(ww,kq);
if t1[1] < -1.65; ik1 = ik1 + 1; goto dasi1; endif;

nk1230 发表于 2015-10-16 09:36
我做出来结果指出下面这段有问题，请指教！！
dasi1:
  ww = yhh[.,ik1:ik1+1];

别客气，谈不上请教，互相学习。
kk = yh[t,.]'~seqa(1,1,n);
kk = rev(sortc(kk,1));
kin = kk[.,2]';
yhh = yh[.,kin];
我认为就是这几个赋值导致了下面问题的出现。kk应该是个矩阵，kin应该是个列向量，而 yhh = yh[.,kin]，那么yhh是什么呢？这块我没弄明白。

问题出现在150-160 行 ， 用debug 功能，把break point 设定在150行

前面几个循环没问题，但是yhh 的columns 一直在变小，最后停止的时候ik1 = 2， columns = 2， 这样yhh[.,ik1:ik1 +1] 就超出范围了
  
我加了两行 code 在里面
            print "ik1 = " ik1;
            print "cols of yhh = " cols(yhh);
你可以run 一下 看看问题

w915w@163.com 发表于 2015-11-12 06:36
问题出现在150-160 行 ， 用debug 功能，把break point 设定在150行

前面几个循环没问题，但是yhh 的col ...

太感谢了，我试试