[求助]GAUSS code 调试

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

lngdpq.rar (9.56 KB) 本附件包括：

• lngdpq.txt
• uc_uni.dta
• uc_uni.out
• uc_uni
• tvp.opt
• uc_uni.opt

2010-11-11 19:07:17 上传

本人在http://www.pinggu.org/bbs/viewthread.php?tid=775317&extra=&page=2下载了版主的一个程序，做了些修改见附件，在GAUSS7下运行，结果却是D:\uc_uni.opt(144) : error G0048 : Matrix singular ，奇怪的是，退出以后再次运行，结果就是c:\gauss7.0\src\optutil.src(1321) : error G0043 : Not implemented for complex matrices ，如果直接运行版主给的原程序uc_uni.opt，结果就是D:\gauss\uc_un\uc_uni.opt(26) : error G0014 : '' : File not found ，我的文件的路径设置的和程序要求是一样的！不知道为什么！求赐教！！！！！！

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：GAUSS code Aus COD ODE 求助 GAUSS code 调试

Matrix singular
本文来自: 人大经济论坛 Gauss专版 版，详细出处参考：http://www.pinggu.org/bbs/viewth ... =1&from^^uid=9913
数据问题吧？你没上传数据。

版主！我重新上传了！请帮忙看一下！！非常感谢！！！！

已修改

版主威武！！！！！！

4# xuehe

版主太强了！另外版主可知道下面这条语句的意思
     p = zeros(2,2);
      p[1,1] = exp(-c0[1]);
      p[2,1] = c0[3];
      p[2,2] = exp(-c0[2]);
      varcov = p'p;
      c1[1]=sqrt(varcov[1,1]);
      c1[2]=sqrt(varcov[2,2]);
      c1[3]=varcov[2,1];
   
        aaa=c0[5]./(1+abs(c0[5]));
      ccc=(1-abs(aaa))*c0[6]./(1+abs(c0[6]))+abs(aaa)-aaa^2;
     c1[5]=2*aaa;
     c1[6]= -1* (aaa^2+ccc) ;
retp(c1);
endp;
其中aaa的赋值好像是将参数变为无约束形式，而ccc的赋值就不清楚了，为什么是这种形式？它是整么得到的？c1[5]和c1[6]为什么会是那种形式的赋值呢？感觉像是与平稳性有关，不知道版主可否解释一下？非常感谢！

这个请你看文献琢磨，不能把程序与程序有关的文献割裂。

我发觉Gauss板块里面的人都是高手。。。。

7# xuehe

谢谢版主帮助！最近学习Kim的书的第四章，遇到这两的程序，运行时都是死循环！版主可知道如何解决？非常感谢！！ hmt.rar (7.82 KB)

2010-11-13 14:23:01 上传

It had already been modified.