[R] constrained optimization using nlm/optim?
Roger D. Peng
rpeng at jhsph.edu
Tue Oct 12 20:15:34 CEST 2004
My guess is that approach will not work particularly well with nlm(),
which expects the objective function to be smooth. You might have
better luck using the Nelder-Mead method in optim().
-roger
Jeff D. Hamann wrote:
> I'm looking for an example of a simple R script that impliments a
> contrained nonlinear function using nlm or optim. I'm not exactly sure how
> to impliment the constraints within the objective function that is passed
> to nlm/optim.
>
> obj.func <- function( p ) {
>
> x(p) <- unconstrained obj function value
>
> if( constraint1 > something ) {
> obj.func <- x(p)
> } else {
> obj.func <- some super huge number
> }
>
> }
>
> p <- c(0.1,2.4, 5)
>
> nlm( obj.func, p, and a bunch of other arguments )
>
> Any suggestions would be very helpful...
>
> Jeff.
>
>
>
More information about the R-help
mailing list