希望通过最大似然函数来做预测，
但是书上提供的，均值只是一个常数
dec vect b(n)
nonlin(parmset=meanparms) b
do i=1,n
   frml resid(i) = (y(&i)-b(&i))
end do i
nlsystem(parmset=meanparms,resids=u) gstart gend resid
compute rr=%sigma
怎么把均值设置成arma形式啊？
dec vect c(4) d(4)
nonlin(parmset=meanparms) c d
frml resid(1)=y(1)-c(1)-d(1)*y(1){-1}
frml resid(2)=y(2)-c(2)-d(2)*y(2){-1}
frml resid(3)=y(3)-c(3)-d(3)*y(3){-1}
frml resid(4)=y(4)-c(4)-d(4)*y(4){-1}
nlsystem(parmset=meanparms,resids=u) gstart gend resid
compute rr=%sigma
但写成这样,但并不可以求大神解答