Displaying 1 result from an estimated 1 matches for "nll21".
Did you mean:
nll23
2007 Jan 10
2
problems with optim, "for"-loops and machine precision
...te("populations/TE23simset_",s,".txt",sep="")),e
nv=MyEnv) #reading a data set
M23=optim(rep(0.1,6),NLL23,method="L-BFGS-B",lower=0,
upper=c(Inf,Inf,Inf,Inf,1,1),control=c(maxit=150))
if (M23$convergence==0)
{
M21=optim(rep(0.1,7),NLL21,method="L-BFGS-B",lower=0,
upper=c(Inf,Inf,Inf,Inf,Inf,1,1),control=c(maxit=150))
}
res[s,1]=M23$par[1]
res[s,2]=M23$par[2]
res[s,3]=M23$par[3]
res[s,4]=M23$par[4]
res[s,5]=M23$par[5]
res[s,6]=M23$par[6]
res[s,7]=M23$value...