求牛人debug

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

proc selif(x,e);
    retp( submat(x,submat(packr(seqa(1,1,rows(x))~miss(e,0)),0,1),0) );
endp;


proc 1=binitial(initd);
local b1, b2, b3, sigma1, sigma2, selctval;

initd = packr(initd);

selctval = selif(initd[.,1:_n1+1],initd[.,1].>=0);
b1 = selctval[.,1] / selctval[.,2:_n1+1];
sigma1 = stdc( selctval[.,1] - selctval[.,2:_n1+1]*b1 );

selctval = selif(initd[.,1 _n1+2:_n1+_n2+1],initd[.,1].<0);
b2 = selctval[.,1] / selctval[.,2:_n2+1];
sigma2 = stdc( selctval[.,1] - selctval[.,2:1+_n2]*b2 );

selctval = initd[.,1] - meanc(initd[.,1]);
selctval = cdfn( selctval./stdc(selctval[.,1]) );
b3 = selctval[.,1] / initd[.,_n1+_n2+2:_n1+_n2+_n3+1];

retp(b1|b2|b3|sigma1|sigma2);
endp;

报错
c:\gauss9.0\pitfalls\swinit2.g(16) : error G0276 : 'selif' : Illegal use of reserved word
c:\gauss9.0\pitfalls\swinit2.g(18) : error G0055 : RETP outside of procedure
c:\gauss9.0\pitfalls\swinit2.g(19) : error G0158 : ENDP without PROC

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：debug bug Initial Sigma Local debug

把selif proc  删除就可以了。