用optim函数

n<-10;
LL<-function(params,data)
{
    theta<-params[1]
    a<-params[2]
    b<-params[3]
ll<-log(1-theta+theta/sum(prod(a+b/i)))+n*log(1-theta+theta/sum(prod(a+b/i))-1+theta)+sum((n-i+1)*(a+b/i))
return(-ll)
}
geyser.res<-nlminb(c(0.5,0.5,0.2),LL,data=rnbinom(n,2,0.5),lower=c(-2,-2,0),upper=c(2,2,1))
geyser.res$par

可以帮忙修改一下吗？

i是什么