Hi:
The MASS package has a function dose.p() to produce a CI for ED50, ED90 or
EDp in general (0 < p < 100). It takes a model object (presumably from a
suitable logistic regression) as input. You could always take the code
already available and adapt it to your situation or you could investigate
one or more of the packages devoted to dose-response models.
A useful thing to know in R is the sos package, starting with its primary
function, findFn():
library(sos)
findFn('ED50')
scares up the three major packages related to dose-response models with 40
matches to ED50.
HTH,
Dennis
On Mon, Nov 22, 2010 at 11:24 AM, Kenney, Colleen T CTR USA AMC <
[email protected]> wrote:
> Classification: UNCLASSIFIED
> Caveats: NONE
>
> A similar question has been posted in the past but never answered. My
> question is this: for probit analysis, how do you program a 95%
> confidence interval for the LD50 (or LC50, ec50, etc.), including a
> heterogeneity factor as written about in "Probit Analysis" by
> Finney(1971)? The heterogeneity factor comes into play through the
> chi-squared test for homogeneity and is equal to h=chi^2/(k-2), where k
> is the number of doses and k-2 are the degrees of freedom.
>
> I have done a lot of research on this and really appreciate any help
> anyone can provide. Thank you!
>
>
> Classification: UNCLASSIFIED
> Caveats: NONE
>
> ______________________________________________
> [email protected] 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.
>
[[alternative HTML version deleted]]
______________________________________________
[email protected] 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.
