在gauss帮助里 搜索 proc 或者 procedure 就可以了

比如，简单回归的最大似然估计。
y = beta0 + beta1 * x +u
如果扰动项u服从正态分布，均值为0，方差sigma2.

下面y表示N*1数据，x表示N*1数据，u表示残差序列，似然函数可以这样：
proc (1)=lik_fcn(parm);
        local beta,sigma2,xx,u,lv;
        beta=parm[1:2];
        sigma2=parm[3];
        xx=x~ones(N,1);
        u=y-xx*beta;
        lv=-0.5*N*ln(2*pi*sigma2)-0.5*u'*u/sigma2;        
        retp(-lv);
endp;

l请问，lik-fcn有下载的没？

你还是没搞清楚！
你这个问题表明了根本不知道optmum 是干嘛用的

额，我要是知道了就不用着急向大家请教了么，或者你可以帮我推荐几本参考书，我现在只有林光平的教材和KIM的代码，不知道从哪里入手学习才好。