s-news
[Top] [All Lists]

4PL model fitting

To: Splus Listserv <s-news@lists.biostat.wustl.edu>
Subject: 4PL model fitting
From: Greg Tarpinian <sasprog474@yahoo.com>
Date: Thu, 12 Oct 2006 17:01:35 -0700 (PDT)
Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=s1024; d=yahoo.com; h=Message-ID:Received:Date:From:Subject:To:MIME-Version:Content-Type:Content-Transfer-Encoding; b=5rTBWXwLRo3LPIyw1zsnoklGI1n/D11B9AJbLOC07DRfXrIdojcZ7h/gCQQ+dKCtLfBb2fSq94MDm2Yf/VBj9qDZ69acD59HvmRFhcKoKqziuy2dQoIloP7S/CbPvaEu+NsHmkC7OTa8JLmOEiJPCwelBK0jwthHi1/TJUnsMHw= ;
WinXP, Splus7.


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

__________________________________________________
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com 

<Prev in Thread] Current Thread [Next in Thread>
  • 4PL model fitting, Greg Tarpinian <=