maximum likelihood estimation程序问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位大侠
我在运行下面程序时，总是提示：   Label undefined G0090 : 'skip' ，错误指向：endp;
下面是程序：
library maxlik;
#include maxlik.ext;
fname="c:\gauss7.0\examples\crisp.w";
range="c5:c1331";
{bp,names}=import(fname,range,1);
var=vcx(bp);
u=(bp-meanc(bp))^2;
maxset;
_max_active={1,1,1};
_max_gradtol=1e-5;
_max_algorithm=5;
_max_linesearch=5;
_max_maxiters=1000;
_max_maxtry=5;
_max_randradious=1e-3;
_max_covpar=3;
__output=1;
iv=0.001|0.97|0.02;
{est,fmax,grad,cov,retcode}maxlik(u,1,&lr,iv);
proc lr(miv,u);
local a0,a1,a2,h,t,ll;
a0=miv[1]; a1=miv[2]; a2=miv[3];
if a0<=0;ll=1000000; goto skip;endif;
if a1<0;ll=-1000000; goto skip; endif;
if a2<0; ll=-1000000; goto skip; endif;
if (a1+a2)>1; ll=-1000000; goto skip; endif;
h=zeros(rows(u),1);
h[1]=a0+(a1+a2)*var;
t=2;
do until t>rows(u);
h[t]=a0+a1*h[t-1]+a2*u[t-1];
t=t+1;
endo;
ll=-.5*ln(2*pi)-.5*ln(h[2:rows(u)])-.5*(u[2:rows(u)]./h[2:rows(u)]);
skip:
retp (lr);
endp;
__title="garch(1,1) estimates";
print maxprt(est,fmax,grad,cov,retcode);

另外，我使用是GAUSS7.0

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Estimation Likelihood Maximum ATION like 程序 Estimation Likelihood Maximum

所编程序在下面文章中可见
GAUSS Programming for Econometricians
Soosung Hwang
Dept. of Banking and Finance,
City University Business School
September 25, 2001
这篇文章在论坛中可以下载