开始对数据做估计第一个数据被系统删除 再做估计就出现division by zero
程序:
@logl huigou
res=huigou-exp(-c(1))*huigou(-1)+c(2)*(exp(-c(1))-1)
var=(c(3)*c(3))*(exp(-2*c(1))-1)*huigou(-1)^(2*c(4))/(-2*c(1))
huigou=log(c(7)*exp(-(res-c(5))^2/(2*var+2*c(6)^2))/sqr(2*3.14*(var+c(6)^2))+(1-c(7))*exp(-(res)^2/(2*var))/sqr(2*3.14*var))