Hi,
I encounter a problem in R. My program is as follows.
lower <- c(-Inf,-Inf,-Inf,-Inf,0,0,0,-1,-1,-1)
upper <- c(Inf,Inf,Inf,Inf,Inf,Inf,Inf,1,1,1)
out <-
optim(parm,logLik,method='L-BFGS-B',lower=lower,upper=upper,hessian=hessian)
As you can see, I have restricted parameter[5], parameter[6] and parameter[7] to
be greater than 0, but the program error is shown in the attached picture
('??' means 'produce'). The reason for this error seems to be
the limit in the function does not work. I don't know why.
Can you help me? Thanks very much!.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: picture.png
Type: image/png
Size: 35455 bytes
Desc: picture.png
URL:
<https://stat.ethz.ch/pipermail/r-help/attachments/20200726/0dc535e9/attachment.png>
This is not reproducible. [1] http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example [2] http://adv-r.had.co.nz/Reproducibility.html [3] https://cran.r-project.org/web/packages/reprex/index.html (read the vignette) On July 26, 2020 7:33:32 AM PDT, Zixuan Qi <zixuan.qi at duke.edu> wrote:>Hi, > >I encounter a problem in R. My program is as follows. >lower <- c(-Inf,-Inf,-Inf,-Inf,0,0,0,-1,-1,-1) >upper <- c(Inf,Inf,Inf,Inf,Inf,Inf,Inf,1,1,1) >out <- >optim(parm,logLik,method='L-BFGS-B',lower=lower,upper=upper,hessian=hessian) > >As you can see, I have restricted parameter[5], parameter[6] and >parameter[7] to be greater than 0, but the program error is shown in >the attached picture ('??' means 'produce'). The reason for this error >seems to be the limit in the function does not work. I don't know why. > >Can you help me? Thanks very much!.-- Sent from my phone. Please excuse my brevity.
For this and the nlminb posting, a reproducible example would be useful. The optimx package (I am maintainer) would make your life easier in that it wraps nlminb and optim() and other solvers, so you can use a consistent call. Also you can compare several methods with opm(), but do NOT use this for production work as it chews cpu cycles. It is intended to test and compare methods for a family of objective functions. JN On 2020-07-26 10:33 a.m., Zixuan Qi wrote:> Hi, > > I encounter a problem in R. My program is as follows. > lower <- c(-Inf,-Inf,-Inf,-Inf,0,0,0,-1,-1,-1) > upper <- c(Inf,Inf,Inf,Inf,Inf,Inf,Inf,1,1,1) > out <- optim(parm,logLik,method='L-BFGS-B',lower=lower,upper=upper,hessian=hessian) > > As you can see, I have restricted parameter[5], parameter[6] and parameter[7] to be greater than 0, but the program error is shown in the attached picture ('??' means 'produce'). The reason for this error seems to be the limit in the function does not work. I don't know why. > > Can you help me? Thanks very much!. > > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. >
Dear Zixuan, On 2020-07-26 07:36 -0700, Jeff Newmiller wrote:> On July 26, 2020 7:33:32 AM PDT, Zixuan Qi <zixuan.qi at duke.edu> wrote: > > Hi, > > > > I encounter a problem in R. My program is as follows. > > lower <- c(-Inf, -Inf, -Inf, -Inf, 0, 0, 0, -1, -1, -1) > > upper <- c(Inf, Inf, Inf, Inf, Inf, Inf, Inf, 1, 1, 1) > > out <- optim(parm, logLik, method='L-BFGS-B', lower=lower, upper=upper, hessian=hessian)Can you include parm?> > As you can see, I have restricted > > parameter[5], parameter[6] and > > parameter[7] to be greater than 0, > > but the program error is shown in > > the attached picture ('??' means > > 'produce'). The reason for this > > error seems to be the limit in the > > function does not work. I don't know > > why. > > > > Can you help me? Thanks very much!. > > This is not reproducible. > > [1] http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example > [2] http://adv-r.had.co.nz/Reproducibility.html > [3] https://cran.r-project.org/web/packages/reprex/index.html (read the vignette)Best, Rasmus -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 833 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20200726/46a617fc/attachment.sig>