2# epoh
谢啦,是说改成这样吗?试了一下还是不行啊
series d1 = 0
smpl @first @first
d1 = 1
smpl @all
@logl logl
@param c(1) 0.3 c(2) 0.5 c(3) 0.1 c(4) 0.05
res=rus-c(1)*rus(-1)
sig2=c(2)+c(3)*res(-1)^2+c(4)*sig2(-1)
z=1+res^2/sig2/(tdf(1)-2)
logl=@gammalog((tdf(1)+1)/2)-@gammalog(tdf(1)/2)-log(3.14159265359)/2-log(tdf(1)-2)/2-log(sig2)/2-(tdf(1)+1)*log(z)/2