Re: [R] Estimate Fisher Information by Hessian from OPTIM ChenYen wrote:
> Dear All, 
> I am trying to find MLE by using "OPTIM" function.
>
> Difficult in differentiating some parameter in my objective function, I
> would like to use the returned hessian matrix to yield an estimate of
> Fisher's Information matrix.
>
> My question: Since the hessian is calculated by numerical differentiate, is
> it a reliable estimate? Otherwise I would have to do a lot of  work to write
> a second derivative on my own.
>
>  
>
> Thank you very much in advance
>
>
> 	[[alternative HTML version deleted]]
>
>
>   
When the objective function is based on a smooth function (in 
particular, a mix of exponentials) then in my experience the Fisher 
information matrix is the same as estimated via the finite difference 
approximation in numericDeriv or via analytical derivatives -- e.g., for 
the results discussed in
Katharine M. Mullen, Mikas Vengris, and Ivo H. M. van Stokkum. 
Algorithms for separable nonlinear least squares with application to 
modelling time-resolved spectra. /Journal of Global Optimization/, vol 
38, n 2, 201-213, 2007 (at http://www.nat.vu.nl/~kate/jgo2005.ps)

______________________________________________
R-help@stat.math.ethz.ch 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.

