Dear Andy and all
Thank you very much for your quick reply.
I misunderstood that SSfpl is for not only numerical response but also binomial
response.
You know that the below data includes binary data for four doses.
In fact, I would like to apply a three parameter logistic model with binary
data.
I believe that glm is for the fitting to linear models with multiple link
functions.
Again, is that my misunderstanding for glm?
Thanks in Advance
-Baba
Liaw, Andy wrote:>Wouldn't you just use glm(..., family=binomial) for binary data? SSfpl
_is_
>for numerical response, where the upper and lower asymptotes need to be
>estimated. Binary data won't need those as parameters--- they are 1 and
0
>automatically!
>
>V&R's MASS has an example, and a function `dose.p()' for
estimating EC50.
>
>Andy
>
>> -----Original Message-----
>> From: baba at muj.biglobe.ne.jp [mailto:baba at muj.biglobe.ne.jp]
>> Sent: Friday, September 27, 2002 12:20 PM
>> To: r-help at stat.math.ethz.ch
>> Subject: [R] How to apply SSfpl with binary data
>>
>>
>> Dear R-help subscribers
>>
>> Would you tell me how to apply SSfpl with binary data as below?
>> Unfortunately, there is not the EXAMPLE in help(SSfpl) for
>> binary data but for quantitative data(Chick).
>>
>> V1: dose
>> V2: log-transformed dose
>> V3: response (rate)
>>
>> V1 V2 V3
>> 1 0.775 -0.2548922 0.1666667
>> 2 5.000 1.6094379 0.8148148
>> 3 10.000 2.3025851 0.5000000
>> 4 20.000 2.9957323 0.6428571
>>
>> I also attach the error message in R 1.5.1 in following.
>>
>> > targan
>> V1 V2 V3
>> 1 0.775 -0.2548922 0.1666667
>> 2 5.000 1.6094379 0.8148148
>> 3 10.000 2.3025851 0.5000000
>> 4 20.000 2.9957323 0.6428571
>> > targan.1 <- nls(targan$V3 ~ SSlogis(targan$V2, Asym, xmid,
>> scal), data = targan)
>> Error in solve.qr(a, b, tol) : singular matrix `a' in solve
>>
>> Best and Kind Regards,
>> Baba
>>
>> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.
>> -.-.-.-.-.-.-.-.-
>> r-help mailing list -- Read
>> http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
>> Send "info", "help", or "[un]subscribe"
>> (in the "body", not the subject !) To:
>> r-help-request at stat.math.ethz.ch
>> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._.
>> _._._._._._._._._
>>
>
>------------------------------------------------------------------------------
>Notice: This e-mail message, together with any attachments, contains
information of Merck & Co., Inc. (Whitehouse Station,
New Jersey, USA) that may be confidential, proprietary copyrighted and/or
legally privileged, and is intended solely for
the use of the individual or entity named on this message. If you are not the
intended recipient, and have received this
message in error, please immediately return this by e-mail and then delete
it.>
>=============================================================================>
>
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._