Re: [R] help in maximum likelihood estimation

2016-03-29 Thread Marc Girondot via R-help
Le 28/03/2016 22:19, heba eldeeb via R-help a écrit : Dear AllI'm trying to find the maximum likelihood estimator of a certain distribution using nlm command but I receive an error as: non-finite value supplied by 'nlm' can't figure out what is wrong in my function Any help? Thank you in ad

[R] help in maximum likelihood estimation

2016-03-28 Thread heba eldeeb via R-help
Dear AllI'm trying to find the maximum likelihood estimator  of a certain distribution using nlm command but I receive an error as:  non-finite value supplied by 'nlm' can't figure out what is wrong in my function Any help? Thank you in advance [[alternative HTML version deleted]] __

Re: [R] help in maximum likelihood

2016-03-04 Thread Ravi Varadhan
Standard error = sqrt(diag(solve(opt$hessian))) Ravi From: Alaa Sindi [mailto:alaasi...@icloud.com] Sent: Wednesday, March 02, 2016 3:22 PM To: Ravi Varadhan Cc: r-help@r-project.org Subject: Re: help in maximum likelihood Thank you very much prof. Ravi, That was very helpful. Is there a way t

Re: [R] help in maximum likelihood

2016-03-02 Thread Alaa Sindi
Thank you very much prof. Ravi, That was very helpful. Is there a way to get the t and p value for the coefficients? Thanks Alaa > On Mar 2, 2016, at 10:05 AM, Ravi Varadhan wrote: > > There is nothing wrong with the optimization. It is a warning message. > However, this is a good example

Re: [R] help in maximum likelihood

2016-03-02 Thread Ravi Varadhan
There is nothing wrong with the optimization. It is a warning message. However, this is a good example to show that one should not simply dismiss a warning before understanding what it means. The MLE parameters are also large, indicating that there is something funky about the model or the da

Re: [R] help in maximum likelihood

2016-03-01 Thread ProfJCNash
It's useful to add "print.level=2" inside your call to find that there's essentially nothing wrong. Rvmmin doesn't give the msg and numDeriv gives a similar (BUT NOT EXACTLY THE SAME!) hessian estimate. It's almost always worthwhile turning on the diagnostic printing when doing optimization, even

[R] help in maximum likelihood

2016-03-01 Thread Alaa Sindi
Hi all, what is wrong with this code? I am trying to estimate the model parameters by maximizing the likelihood function and I am getting this warning Warning message: In nlm(fn, p = c(-50, 20), hessian = TRUE) : NA/Inf replaced by maximum positive value x <- c(1.6907, 1.7242, 1.7552, 1