[Rd] nmle: gnls freezes on difficult case

From: Nicholas Lewin-Koh <nikko_at_hailmail.net>
Date: Thu, 18 Oct 2007 10:17:50 -0400


Hi,
Following up on my own post, if in gnlsControl I specify opt='optim' gnls exits quite nicely, with an error, which is what I would expect. Is this a bug in nlminb?

Nicholas

tt<-gnls(response~SSllogis(conc,A,B,xmid,scal),tdat,start=start,weights=varPower(),verbose=TRUE,control=gnlsControl(opt='optim')) **Iteration 1
GLS step: Objective: 29.04254varStruct parameters:

    power
0.3373206

NLS step: RSS = 0
 model parameters:-0.799941 8.99983 -0.522623 212.314  iterations: 2

Convergence:

   params varStruct
 1.172208 1.000000
Error in optim(c(coef(gnlsSt)), function(gnlsPars) -logLik(gnlsSt, gnlsPars), :
  initial value in 'vmmin' is not finite In addition: Warning messages:

1: In log(xmid) : NaNs produced
2: In log(xmid) : NaNs produced
3: In log(xmid) : NaNs produced

______________________________________________
R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Thu 18 Oct 2007 - 14:25:46 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 25 Oct 2007 - 11:37:11 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.