[R] 4PL algorithm

From: Greg Tarpinian <sasprog474_at_yahoo.com>
Date: Thu 12 Oct 2006 - 23:36:42 GMT


WinXP, Splus7 and R2.3.1.

All,

I have been using the SSfpl and SSlogis self-starting functions in the nlme library to fit 4PL and restricted 4PL models. I need to adapt these routines to fit the alternative model

   f(x) = A + (B-A)/(1 + abs(x/EC50)^C)

My Question: How do I obtain good starting values for this alternative model?

(The pseudo-code found on pages 517 - 520 of "Mixed Effects Models in S and Splus" is not applicable to my problem because it deals with f(x)=A+(B-A)/(1+exp(EC50-x)/C)) which is not the same function.)

Many thanks,

      Greg



R-help@stat.math.ethz.ch 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. Received on Fri Oct 13 09:41:07 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Fri 13 Oct 2006 - 00:30:09 GMT.

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