Your function is evaluating to Inf. Is this expected behavior? I don't
know what starting values you are using, but it looks like the initial
values evaluate to Inf. If this is also true for values near your initial
conditions then I am not surprised the optimization returns initial
values.
Try changing parametres[,3]
--------------------------------------
Jonathan P. Daily
Technician - USGS Leetown Science Center
11649 Leetown Road
Kearneysville WV, 25430
(304) 724-4480
"Is the room still a room when its empty? Does the room,
the thing itself have purpose? Or do we, what's the word... imbue it."
- Jubal Early, Firefly
kamel gaanoun <kamel.gaan...@gmail.com> wrote on 12/15/2010 09:22:57 AM:
> [image removed]
>
> Re: [R] Help about nlminb function
>
> kamel gaanoun
>
> to:
>
> Jonathan P Daily
>
> 12/15/2010 09:23 AM
>
> Cc:
>
> r-help, r-help-bounces
>
> I appologize for the lack of information in my previous mail, but it
> is my first one.
>
> Yes I red -in part- the Port documentation to learn what does the
> message X-convergence (3) mean.
>
> so this is my nlminb and the results :
>
> res1<-nlminb(parametres[,3], V, lower=c(rep(0.01, 18), rep(0.01, 3),
> rep(-Inf, 6)), upper=c(rep(Inf, 18), rep(0.99, 3), rep(Inf, 6)),
> control = list(maxit=1000) )
>
>
> Results : > res1
> $par
> [1] 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000
1.0000000
> [8] 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000
1.0000000
> [15] 1.0000000 1.0000000 1.0000000 1.0000000 0.6899225 0.7993528
0.6480687
> [22] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
>
> $objective
> [1] Inf
>
> $convergence
> [1] 0
>
> $message
> [1] "X-convergence (3)"
>
> $iterations
> [1] 1
>
> $evaluations
> function gradient
> 2 27
>
>
> So as you see, it stops after just one iteration but it is converging.
>
>
> Thanks
> 2010/12/15 Jonathan P Daily <jda...@usgs.gov>
> How are we supposed to help you if we have no idea how you are using it?
> Did you read the PORT documentation? Did you look at the control
> parameters outlined in the function help page? Often times for me a
return
> of the initial conditions is an issue with either inappropriate
tolerance
> or having variables on very different scales.
>
> Since you did not follow the posting guide (important part *'d below),
my
> only suggestion is to ?nlminb
> --------------------------------------
> Jonathan P. Daily
> Technician - USGS Leetown Science Center
> 11649 Leetown Road
> Kearneysville WV, 25430
> (304) 724-4480
> "Is the room still a room when its empty? Does the room,
> the thing itself have purpose? Or do we, what's the word... imbue it."
> - Jubal Early, Firefly
>
> r-help-boun...@r-project.org wrote on 12/15/2010 08:53:39 AM:
>
> > [image removed]
> >
> > [R] Help about nlminb function
> >
> > kamel gaanoun
> >
> > to:
> >
> > r-help
> >
> > 12/15/2010 08:56 AM
> >
> > Sent by:
> >
> > r-help-boun...@r-project.org
> >
> > Hi Everyone,
> >
> > Can anyone help me resolve a problem that i'm having with nlminb.
The
> > problem is that it stops after just one iteration and returns the same
> > values as "start" ones.
> >
> >
> > Thank you very much for your help.
> >
> >
> > Sincerely.
> >
> >
> > --
> > Kamel Gaanoun
> > (+33) (0)6.76.04.65.77
> >
> > [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help@r-project.org 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.
>
******************************************************************************
>
>
>
> --
> Kamel Gaanoun
> (+33) (0)6.76.04.65.77
______________________________________________
R-help@r-project.org 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.
