Re: [R] question about optim

2011-04-14 Thread chirine wolley
is still there. I also made sure that the function doesnt return NAN or INF values... Does anyone have any idea where the problem could come from ? Many thanks. To: r-h...@stat.math.ethz.ch From: bbol...@gmail.com Date: Tue, 12 Apr 2011 17:24:05 + Subject: Re: [R] question about optim

Re: [R] question about optim

2011-04-12 Thread Ben Bolker
chirine wolley wolley.chirine at hotmail.com writes: Dear R-users, I would like to use optim( ) to minimize a function which depends on 4 parameters: 2 vectors, a scalar, and a matrix. And I have a hard to define the parameters at the beginning of the function, and then to call optim.

Re: [R] question on optim() fn.

2010-10-21 Thread Ravi Varadhan
of Geriatric Medicine and Gerontology School of Medicine Johns Hopkins University Ph. (410) 502-2619 email: rvarad...@jhmi.edu - Original Message - From: sharkoil whitewhale2...@gmail.com Date: Thursday, October 21, 2010 1:09 am Subject: Re: [R] question on optim() fn. To: r-help@r-project.org

Re: [R] question on optim() fn.

2010-10-20 Thread Ravi Varadhan
Hi, You are doing the right thing by being concerned about the warning messages. The warning message of `NaN' usually arises due to illegal arithmetic operations (in the real number field) such as logarithm or square root of a negative number. This could happen during the intermediate stages of

Re: [R] question on optim() fn.

2010-10-20 Thread sharkoil
Hi Ravi, Thanks a lot for your informations. I think you are right. Without the constraints of the parameters, it's easy to get the NaN Produced warning message. The pdfs I worked on are Burr, Lognormal and Inverse Gaussian distributions. Lognormal distribution is ok, but I always have warning

Re: [R] question on optim

2010-09-15 Thread Hey Sky
thanks. Ravi and Nash. I will read the new package and may use it after I am familiar with it. I may bother both of you when I have questions.thanks for that in advance. Nan from Montreal Hi Nan, You can take a look at the optimx package on CRAN. John Nash and I wrote this package to

Re: [R] question on optim

2010-09-08 Thread Ravi Varadhan
...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Hey Sky Sent: Tuesday, September 07, 2010 2:48 PM To: Ben Bolker; r-h...@stat.math.ethz.ch Subject: Re: [R] question on optim thanks. Ben after read your email, I realized the initial value of w[5]=0 is a stupid mistake. and I

Re: [R] question on optim

2010-09-08 Thread Paulo Barata
-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Hey Sky Sent: Tuesday, September 07, 2010 2:48 PM To: Ben Bolker; r-h...@stat.math.ethz.ch Subject: Re: [R] question on optim thanks. Ben after read your email, I realized the initial value of w[5]=0 is a stupid mistake. and I

Re: [R] question on optim

2010-09-08 Thread Ravi Varadhan
...@infolink.com.br] Sent: Wednesday, September 08, 2010 4:26 PM To: Ravi Varadhan Cc: 'Hey Sky'; 'Ben Bolker'; r-h...@stat.math.ethz.ch Subject: Re: [R] question on optim Dear R-list members, I am using R 2.11.1 on Windows XP. When I try to install package optimx through the GUI menu Packages / Install

Re: [R] question on optim

2010-09-07 Thread Hey Sky
sorry. there is a type in the following code. there is no w[5]*acwrk[,i] in the regw equation. the right one should be as following:    regw[,i]-w[1]+ w[2]*eta[m]+exp(w[3]+w[4]*eta[m])*actr[,i]   - Original Message From: Hey Sky heyskywal...@yahoo.com To: R r-help@r-project.org

Re: [R] question on optim

2010-09-07 Thread Ravi Varadhan
Hi, I do not see how `data' is used in your objective function. The objective function is not even evaluable at the initial guess. myfunc1(guess, mydata) [1] NaN I also think that some of the parameters may have to be constrained, for example, par[1] 0. At a minimum, make sure that the

Re: [R] question on optim

2010-09-07 Thread Hey Sky
thanks. Ben after read your email, I realized the initial value of w[5]=0 is a stupid mistake. and I have changed it. but I am sorry I cannot reproduce the result, convergence, as you get. the error message is non-finite finite difference value [12]. any suggestion about it? and could you plz

Re: [R] question on optim

2010-09-07 Thread Ben Bolker
On 10-09-07 02:48 PM, Hey Sky wrote: thanks. Ben after read your email, I realized the initial value of w[5]=0 is a stupid mistake. and I have changed it. but I am sorry I cannot reproduce the result, convergence, as you get. the error message is non-finite finite difference value [12]. any