[R]- Estimate LD50 with bnlr{gnlm}

Jim Lindsey james.lindsey at luc.ac.be
Mon Oct 30 08:45:13 CET 2000


> 
> Hi,
> 
> I'm not yet familar with GLM and still learning.
> How can I perform a BNL (to estimate LDp values) with matrixes like this
> (N indicates the observed objects):
> 
> data from:
> Kerr D, Meador J. 1996. "Modeling Dose Response Using Generalized Linear
> Models". Environ Toxicol Chem 15(3)
> 
> conc	respond	n
> 0	0	20
> 0	1	19
> 0	1	20
> 0	0	20
> 13.5	5	23
> 13.5	2	20
> 29	9	20
> 29	6	20
> 53	12	20
> 53	15	20
> 85	19	20
> 85	18	20
> 110	18	18
> 110	18	18
> 
> The bnlr{gnlm} documentation talks about a two column matrix as input.  

The two column matrix would be created by
cbind(n,n-respond)
However, if you are using standard linear logistic or probit models, I
would strongly recommend using glm in R rather than bnlr in my gnlm
library.
  Jim

> 
> Thanks for any help
> Jens
> 
> --
> Jens Schmidt                    TU Dresden -Inst. f. Hydrobiologie-
> Tel. +49-(0)351-463-6915        Mommsenstrasse 13
> Fax  +49-(0)351-463-7108        01062 Dresden
> 
> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> 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
> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
> 

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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