[R] Confidence interval on parameters from optim function

Arun.stat arun.kumar.saha at gmail.com
Wed Aug 19 17:49:40 CEST 2009


Regrading your second question, I guess somehow you get undefined value like
logarithm of zero of your target function for some unfortunate parameter
values in the parameter space.



Devred, Emmanuel wrote:
> 
> Hi everyone,
> 
> I have two questions:
> 
> I would like to get confidence intervals on the coefficients derived
> from the optim() function.
> I apply optim() to a given function f
>> res <-
> optim(c(0.08,0.04,1.),f,NULL,method="L-BFGS-B",lower=c(0.,0.,0.))
> And I would like to get the p-value and confidence intervals associated
> with 
>> res$par
> 
> My second question deals with error message. I am doing a loop with the
> optim() function in it, when I get an error message like below, the loop
> is stopped, however I would like to change my initial values to avoid
> this error message and keep the loop going, so if it crashes when I am
> away the program can still run, any idea ?
> 
> Error in optim(c(0.08, 0.04, 1), f, NULL, method = "L-BFGS-B", lower =
> c(0,  : 
>   L-BFGS-B needs finite values of 'fn'
> 
> Thank you for any information on these two problems.
> 
>   Emmanuel
> 
> -------------------------------------------
> Dr. Emmanuel Devred
> Bedford Institute of Oceanography,
> 1 Challenger Drive,
> Dartmouth, Nova Scotia, B2Y 4A2,
> Canada
> 
> Ph:  (1) 902 426-4681
> Fax: (1) 902 426-9388
> 
> devrede at mar.dfo-mpo.gc.ca
> 
> http://myweb.dal.ca/edevred/
> -------------------------------------------
> 
> 
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at 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.
> 
> 

-- 
View this message in context: http://www.nabble.com/Confidence-interval-on-parameters-from-optim-function-tp25044388p25046772.html
Sent from the R help mailing list archive at Nabble.com.




More information about the R-help mailing list