Displaying 1 result from an estimated 1 matches for "d1j".
Did you mean:
d1
2012 Nov 04
1
Struggeling with nlminb...
...e.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 issues for some d3.
for(j in (my.d1j+1):7){
cat(j,"\n")
d3 <- emp.quant[j]
constLo=c(rep(0.00001, (length(par.start)-1)), -999999)
my.optpar3 <- nlminb(par.start, my.loglike.normal, lower=constLo,
control=list(eval.max=60,iter.max=30))
value <- valore.normal(my.optpar3$par)
}
Thank...