search for: optpar3

...y.sigmasq[i-1]>d3) } my.mean<-rep(0,n+1) for(j in 2:(n+1)) { my.mean[j]<-theta[13]*data[j-1] } 1/2*sum(log(my.sigmasq[2:(n+1)])) + n/2*log(2*pi) + 1/2*sum((data[2:(n+1)]-my.mean[2:(n+1)])^2/(my.sigmasq[2:(n+1)])) } constLow=c(rep(0,(length(par.start)-1)),-2) my.optpar3<- nlminb(par.start,my.loglike.normal,lower=constLow,control=list(eval.max=500,iter.max=100) ) Second, I estimate a similar function but with only 7 instead of 13 parameters, I fix theta[1]-theta[6] to some constant, but vary d3 in a loop. It seems like that the optimizer faces some NA/Inf is...

...the function, as it causes the following error message Fehler (error)* in while (nf < maxfeval & restarts < restarts.max & dist > ftol & : Fehlender Wert (missing value)* , wo (where)* TRUE/FALSE n?tig ist (is required)* *translation Do I need to adjust the control ? my.optpar3<-nmkb(par=par.start,fn=my.loglike.normal,lower=constLo,upper=999, control=list(restarts.max=9)) If you need more info, please let me know. I would appreciate some help, thank you very much. Best, Marcial -- View this message in context: http://r.789695.n4.nabble.com/Error-message-in-nmkb-t...