Dear list,

I need to make nls() or gnls() function to return something like "NA" when
it fails to fit. I found I can use "warnOnly", but I cannot make it work.

In the following program. I intentionally make it fail by making very wrong
initial values, but the "out" is empty. What is wrong with my approach?

Many thanks in advance!

library(nlme)
out<- gnls(conc~SSbiexp(time,A,lalpha,B,lbeta),
          start=list(A=-1,lalpha=2,B=20,lbeta=4),
          weights=varPower(fixed=1),control=list(warnOnly=TRUE),
          data=Indometh,subset=Subject==1)




-- 
Chaehyung Ahn, Ph.D.
4033 Remington Oaks Circle
Cary, NC 27519
Daum Cafe: http://cafe.daum.net/biometrika
Blog: http://cahn88.blogspot.com/

	[[alternative HTML version deleted]]


