[R] How to apply SSfpl with binary data
baba@muj.biglobe.ne.jp
baba at muj.biglobe.ne.jp
Fri Sep 27 20:25:37 CEST 2002
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list