Thanks Rui, i already tried changin the reltol, but it didn't change the outcome
On Wed, May 15, 2013 at 1:32 PM, Rui Barradas <ruipbarra...@sapo.pt> wrote: > Hello, > > It's impossible to tell what is happening without the function and the > values for the other variables (including the initial parameters) but your > setting of reltol is too small. Try using the default. It's > sqrt(.Machine$double.eps), about 1e-8, you are using 1e-16. > > Hope this helps, > > Rui Barradas > > Em 15-05-2013 17:02, Luis Felipe Parra escreveu: > >> Hello to all, >> >> I have been using an optim with the following call: >> >> optim(param_ini,fun_errores2,**Precio_mercado=Precio,anos_** >> pagosE2=anos_pagos,control=**list(maxit=10000,reltol=1e-16)**) >> >> depending on the intial values I'm getting the same solution but once I >> get >> the convergence message=10 (no convergence) and for the others I get >> convergence message = 0 >> >> Solution1: >> $par >> beta1 beta2 beta3 lambda >> beta1 0.06537061 0.001474821 -0.0759236 0.5 >> >> $e2 >> [1] 74.84273 >> >> $conv >> [1] 0 >> >> $v >> [1] 74.84273 >> >> Solution2: >> $par >> beta1 beta2 beta3 lambda >> 1 0.06537061 0.001474822 -0.0759236 0.5 >> >> $e2 >> [1] 74.84273 >> >> $conv >> [1] 10 >> >> $v >> [1] 74.84273 >> >> My intuition tells me the correct solution is no convergence. Does anybody >> know why this might be happening. >> >> Thank you >> >> Felipe Parra >> >> [[alternative HTML version deleted]] >> >> ______________________________**________________ >> R-help@r-project.org mailing list >> https://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/mailman/listinfo/r-help> >> PLEASE do read the posting guide http://www.R-project.org/** >> posting-guide.html <http://www.R-project.org/posting-guide.html> >> and provide commented, minimal, self-contained, reproducible code. >> >> -- Este mensaje de correo electrónico es enviado por Quantil S.A.S y puede contener información confidencial o privilegiada. This e-mail is sent by Quantil S.A.S and may contain confidential or privileged information [[alternative HTML version deleted]]
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
